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IMPROVED SONIC-BOX COMPUTER PROGRAM FOR CALCULATING 
TRANSONIC AERODYNAMIC LOADS ON OSC 1 LL'AT'i NG, W-I NGS WITH- 'TH ICKNESS 

By S. Y. Ruo 

Lockheed-Georgia Company 


. SUMMARY 


A computer program has been developed to account approximately for the 
effects of finite wing thickness in the transonic potential flow over an 
oscillating wing of finite span. The program is based on the orig.inal sonic- 
box program of Rodemich and Andrew which has previously been extended by Ruo 
and Theisen to account for the effects of wing thickness. The latter work is 
further extended herein to improve computational efficiency and accuracy and 
to account for swept trailing edges. Account for the non-uniform flow caused 
by finite thickness is made by application of the local linearization concept 
with appropriate coordinate transformation. A brief description of -each com- 
puter routine and the applications of cubic spline and the spline-surface 
data fitting techniques in the program are given } and the method of input is 
shown in detail. Sample calculations as well as -a complete listing of the 
computer program are presented.' 

INTRODUCTION 

In reference 1, the • son i c-box method computer program was developed for 
calculation of unsteady transonic flow aerodynamics for oscillating planar 
wings with unswept trailing edge by approximating the wing planform with a 
matrix of square boxes. Later, i.t was extended to include the swept trailing 
edge and control -surfaces in reference 2. Sonic-box method uses a velocity- 
potential doublet as the basic solution to satisfy the linearized transonic 
flow, unsteady small-perturbation velocity-potential equation with the associ- 
ated boundary conditions. '• 

In reference 3 , the wing thickness effect is partially recovered by the 
inclusion of local Mach number in the governing equation for the unsteady 



transonic flow. It uses the concept of local linearization to reduce the 
nonlinear .smal 1 -perturbat ion equation to a linear one wi th non-constant 
coefficients. This is further reduced to a linear equation with constant 
coefficients by the appropriate coordinate transformation. This final equa- 
tion and the associated boundary condition in the transformed space become 
identical to those treated in the physical space by Rodemich and Andrew in 
reference 1. The numerical results for the wing with thickness were obtained 
by adopting the sonic-box method in the transformed space. Because of the 
assumptions made in deducing the governing equation to a manageable form, this 
technique is applicable only to the thin wings. That is, the local mean Mach 
number on the wing surface must not be very different from unity. Further, 
it is assumed that there is no flow separation and no strong shock waves on 
the wing surface. 

IT 

The computer programs developed in references 1, 2, and 3 use the least- 
square method to fit some of the input data, such as wing deflection or steady 
Mach number distribution on the wing,, and to fit the computed velocity poten- 
tial with a form of predetermined polynomial surface for the subsequent 
calculation of the unsteady pressure and the generalized aerodynamic force 
coefficients. For wings with a rapidly changing Mach number distribution or 
wing motions in certain deflection mode shapes, the polynomial surface is not 
adequate because of its inaccuracy. 

The computer program described in this report adopts the natural cubic 
spline for fitting calculated velocity potential and the spline-surface for 
fitting input modal deflections and Mach number distribution instead of the 
polynomial -surface fitting used in references 1, 2, and 3. The present com- 
puter code allows the computation of generalized aerodynamic force coefficients 
for wings of zero and finite thickness; the swept trailing edges are allowed 
but not the control surfaces. It also uses a different scheme from that of 
reference 3 for calculation of the generalized aerodynamic force coefficients 
for wings with finite thickness. to improve the efficiency of the program and 
the accuracy of the numerical results. This is done by performing the force 
coefficient integration in the transformed space where the velocity potential 
is computed rather than map the ve loci ty potential from the transformed space 
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to the physical space before performing the integration as in reference 3- 
Additionally, the contribution from the partial boxes along the swept leading 
edge has been included to improve convergence and accuracy of the numerical 
results with a smaller total number of boxes in representing the wing 
planform. 


The wing is assumed to have zero mean angle of attack, and the same mean 
Mach number distribution on upper and lower surfaces, and to perform motions 
symmetric with respect to the root chord. These, however, are not limited 
by the basic concept. 

The regular output of the present computer code is the generalized 
aerodynamic force coefficients. In addition to. these, the values of down- 
wash, velocity potential, pressure coefficient and Mach number at each box 
center may be printed out. Also the portion of the dimensionless area 
representing the actual wing planform in each box and the information of the 
arrangement of the boxes used in the computation as well as the velocity 
influence coefficients may be printed out. Since no provision is made in 
the computer code to smooth any computed data, the pressure coefficient must 
be used with caution. 


SYMBOLS 


a 0> a l» a 2 



b 

C 

d 


P 

j 


e, exp 
f (x,y) 

Gi (x,y) ,G£ (y) 


constants 

area of ij-th box on wing (i-th spanwise column and j-th 
chordwise row) 

reference length (dimension = L) 
pressure coefficient 

x-di recti on distance between two adjacent points used in 
the spline-curve fit, equation (14) 
exponential function 
wing deflection mode shape 

leading edge adjustment terms, see equation (10) 
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H(x,y) 



t 

T 

Uoa 

w(x,y) 
x,y,z 
x £e (y) 

x te (y) 


Ymax 



x 


‘Potx.y) 

4»(x,y) 


n 

n 


dependent variable represented by a cubic spline or spline- 

surface, see equations (14) and (15) 

constant 

integrated value, see equation (13) 

/T 

reduced frequency, mb/U w 
unit of length 

generalized aerodynamic force coefficient 

local Mach number 

coefficients of cubic spline 

number of points used in data fitting 

dimensionless distance between two points, / (x-x; ) 2 +(y-y \ ) 2 , 
(reference length = b) 

dimensionless wing planform area of full wing (reference 
area = b 2 ) 

dimensionless time (reference time = b/Uoo) 
unit of time 

reference velocity (freestream) , (dimension = L/T) 
dimensionless downwash (reference velocity = Um) 
dimensionless Cartesian coordinates (reference length = b) 
dimensionless x-coordinate of wing leading edge (reference 
length = b) 

dimensionless x-coordinate of wing trailing edge (reference 
length = b) 

dimensionless maximum semi-span (reference length = b) 

phase angle of L;j 

maximum thickness to chord ratio 

dimensionless velocity potential of doublet 

downwash at point (x,y) due to doublet of unit strength at 

the origin, see equation (6) 

angular velocity (dimension = radian/T) 

superscript denotes the dimensional quantity of the correspond- 
ing dimensionless variable shown inside the parentheses 
superscript denotes the quantity in the transformed space 
of the corresponding variable shown inside the parentheses 
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MAJOR VARIABLES IN COMPUTER PROGRAM 


Some variables used in the computer program are described briefly in 
this section. The subroutine name in the parentheses following the descrip- 
tion of a variable indicates where its value is generated or defined. The 
number appearing in the parentheses indicates the value of the variable for 
the corresponding condition described. In the following, N=1 and -N=2 indicate, 
respectively, the real and the imaginary parts of a variable; NEW=1 -and NEW=2 
indicate, respectively, the zero thickness and nonzero thickness cases; i 
(i < MB) is related to the chordwise coordinate while j ( j < MB) is related 
to the spanwise coordinate. 


A ( N , j , i ) 


AMA(j,i,NEW) 


AR( j , i , NEW) 

AREA 

CK 

CTE(j) 


D 

DA ( k) 

DH 
D I 

EDG(NEWJ) 


influence coefficient, the integral in equation (6); the 
upwash at the center of a box caused by a unit doublet 
distribution over another box separated from the former 
by (j-1) number of boxes in spanwise direction and (i-l) 
number of boxes in chordwise d i rect ion, (P0T2) 

Mach number at center of ij-th box on physical wing (NEW=1), 
Mach number at a point on physical wing corresponding to 
the center of ij-th box on the transformed wing (NEW=2), 
(SHAPE) 

area of wing in ij-th. box, (PLNFM) 

area of the physical, .full wing planform, (PLNFM) 

reduced frequency, (MAIN) 

Mach number on the physical wing at a point corresponding 
to the trailing edge of the j-th chordwise row of the 
transformed wing, (SHAPE) 
length of box side, (MAIN) 
input data, (DATRD) 

one-half of the length of box side, (MAIN) 

maximum number (real) of boxes in streamwise direction 

(same as L integer), (MAIN) 

value of leading edge adjustment term in equation (10) 
at j-th chordwise row, (SHAPE) 
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I EDG 

IR 

IW 

JMAK(NEW) 
JMAX (NEW) 

KSFO(M) 

KSFM 

L 

M 

HB 

MO 

ML(NEW, i) 
MLC(NEW,k,i) 
MLT ( N EW , k , j ) 
MLW(NEW,-i) 

NB 

NEDGI 

NEW 


flag to identify whether the most outboard section of the 

leading edge is parallel (l) or is not parallel (0) to the 

freestream, (SHAPE) 

computer read-unit number, (MAIN) 

computer write-unit number, (MAIN) 

number of boxes on the wing along spanwise direction, (SHAPE) 
number of boxes on wing plus wake along spanwise direction, 
(SHAPE) 

number of points used in spline-surface fit of wing 
deflections in M-th mode, (DRED) 

number of points used in the spline-surface fit of 
Mach number, (MRED) 

maximum number (integer) of boxes in streamwise direction, 
(MAIN) 

deflection mode number, M < MD, (MAIN) 

maximum number of boxes in streamwise direction corresponding 
to the dimension of those subscripted variables related to 
the box distribution, (MAIN) 

maximum number of deflection modes allowed, corresponding to 
the dimension of those subscripted variables related to 
def feet ion mode, (MAIN) 

number of boxes, including partial boxes, in i-th spanwise 
column, (PLNFM) 

sequential number of first (k=l) and last (k=2) wing box 
in the i~th spanwise column, (PLNFM) 

sequential number of first (k=l) and last (k=2) wing box 
in the j-th chordwise row, (SHAPE) 

sequential number of first wake box next to the wing box 
in i-th spanwise column, (PLNFM) 

number of points to be used in the spline-surface fit, 
corresponding to the dimension of those subscripted 
variables related to the spl ine-surface , . (MAI N) 
maximum number of points allowed to describe the in- 
put leading edge of the wing, (MAIN) 

index for physical (1) or transformed (2) wings, (MAIN) 



NM 

NS (NEW) 


NSMAX 


NST(NEW) 


NTDGl 


S(N,j,i) 


SFDH (k, M,NEW) 


SFDX(k,M) 

SFDY (k,M,NEW) 

SFMH(k) 


SFMX(k) 

SFMY(k) 

T(k,i) 

XEDG(k) 

YEDG(k)' 

XEDG 1 (k) 
YEDG I (k) 

XLE(NEW,j) 
XTE (NEW, j) 

XTDG(k) 

YTDG(k) 


NB+3, (MAIM) 

number of leading edge points on the wing semispan, 
including point at apex only if the wing is pointed » (SHAPE ) 
maximum number of leading-edge points (including apex of 
pointed wing) that can be used to approximate leading edge 
of transformed wing corresponding to the dimension of those 
subscripted variables related to the leading-edge points 
used in computation, (MAIN) 

number of trai 1 ing-edge points including- the point at wing 
root, NST(l) < NTDGI-1, (SHAPE) 

maximum number of points allowed to describe the input 
trailing edge of the wing, (MAIN) 

downwash at the center of ij-th box, equation (8), (WVAL) , 
and velocity potential at the center of ij-th box, (BOXPO) 
known wing deflection at those points used in spline- 
surface fit of wing deflections in M-th mode; or, the 
coefficients in wing-deflection, spline-surface expression 
for M-th mode, k < NM, (DRED) 

x and y coordinates of points used in spline-surface fit of 
wing deflections in M-th mode, k < NB, (DRED) 
known Mach number at those points used in spline-surface fit' 
of Mach number; or, the coefficients in Mach number spline- 
surface expression, k < NM, (MRED) 

x and y coordinates of points used in spline-surface fit of 
Mach number, k < NB, (MRED) 

dimensioned temporary storage, k <_ NM and £ .£ (NM+l) 
x and y coordinates of those points used to describe the 
wing leading edge in the computations, k < NSMAX, (SHAPE) 
x and y coordinates of those input points used to describe 
the wing leading edge in physical space, k < NEDGI, (SHAPE) 
x-coordinate of leading and trailing edge along j-th 
chordwise row, (PLNFM) 

x and y coordinates of those points used to describe the 
wing trailing edge in the computations, (SHAPE) 



XTDGl(k) 
YTDGI (k) 

XX(i) 

YMAX (NEW) 
YY(j) 


x and y coordinates of those input points used to describe 
the trailing edge in physical space, k < NTDG I , (SHAPE) 
x-coordinate of i-th spanwise column, (SHAPE) 
maximum y dimension of the wing, (SHAPE) 
y-coordinate of j-th chordwise row, (SHAPE) 


EQUATIONS USED IN THE COMPUTER PROGRAM 


In this section, some equations and expressions used in the present 
version of sonic-box method computer program are listed. Additional in- 
formation on the derivation of these equations can be found in the standard 
text books or numerous other publications (e.g., refs. 1, 4, and 6). 

Small Perturbation Potential Equation 

The governing equation of the unsteady part of the harmonically 
oscillating thin wing in a sonic flow field can be written as 

% y + Vz ' M2(2ikCp O x * = 0 0) 


and 


^ (x ,y , z , t ) 


<c 0 (x,y,z) 


e ikt 


When M = 1 everywhere, equation (l) reduces to the familiar linearized 
equation. However, if the variation of local Mach number is taken into 
account, local linearization may be applied to obtain the following equation 
valid in a small region 


where 


^Oyy + - (2ikT> Q ~ - k*P Q ) = 0 

•(•2) 

X = X 

(3) 

y = My 


z = Mz 



$ * M9 .. 






o? 
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Equation ("2) is a linear equation with constant coefficients and is the same 
as the familiar linearized equation with M=1 (eq.(l)). Therefore, the prob- 
lems of w'ing with finite thickness can be solved in the transformed space 
with a method suitable for the problems of zero thickness wing in the physical 
space. The sonic-box method (ref. 1) has been chosen for the present 
implementation. A solution to the linearized form of equation ( 1 ) , when M=1 , 
can be written as (ref. 1) 

fo, x < 0 


c P 0 (x,y,z) = 


— 

|2tt 


z r i k / 

—*■ exp L — (x + 


y 2 -f 7 ?- 


•)] , x>0 


( 4 ) 


and in the wake 


^oCxjY.O) = cp 0te exp[-ik(x-x te ) ] , x>x te and z = 0, (5) 

which represents a doublet at the origin of coordinates pulsating with reduced 
frequency k. An expression for downwash at any point on a wi ng 1 y i ng approximate- 
ly in the xy-pl ane due to a un i t doublet at the coordinate origin can be.obta ined 
by taki ng part 5 al d i fferent i at ion of equat ion (4) wi th respect to z and letting 
z go to zero. The approximate solution to equation (2) with a known down- 
wash distribution on a lifting surface, which is approximated by a matrix of 
square boxes and the velocity potential is considered to be constant within 
each box and equal to that at the box center, can be expressed as: 


l 




^(xj -£,yj - n)d£dr| 


w(x; ,yj) 


( 6 ) 


where i * j 1 - magnitude of velocity potential at the center of box B j i j i 


^(xj-^yj-n) = downwash at (x;,yj) due to doublet of unit 
strength at (£,n) . 


Us 1 

2 Tx 


FU? exp {- iik [ (x,- 5 ) + 


( 7 ) 


w(xj ,y i ) = (-g— ) 

J dz at (x],yj) 


= downwash at the center of box B: 


‘J 
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Downwash Is related to the deflection mode shape in physical space as 


9q> 0 

dz 


3f 

3x 


+ ikf 


(8a) 


and in the transformed space as 



+ ik «- 


(8b) 


Since the left-hand sides of equations (8a) and (8b) are identical, the down- 
wash at (x,9) on the transformed wing is unchanged from that at the correspond- 
ing point, (x,y), on the physical wing. 


It is not possible to arrange the center of the box to follow the wing 
leading edge with the sonic-box method. Before the pressure coefficient and 
the generalized aerodynamic force coefficients are calculated, it is recommend- 
ed to adjust the computed velocity potential near the leading edge to conform 
to the theoretically known leading edge condition (e.g., = 0 for a swept 

(subsonic) leading edge). This adjustment is expressed as 


<P’o = G i * G 2 ♦ <P C 


( 9 ) 


where Gi and G 2 , depending solely on the shape of the wing leading edge, are: 


x - x £e sAx: Gi = 1.0 


x - X£ e <Ax: 


! x ~ x Ze f I x+x &e ~ \ m /m=0 for at x=0 
N Ax \ tl 2xjgf-Ax / ’ \m=l for y=0 at x-0 


Ymax-Y^^y : 


Ymax~Y <A Y : 


G 2 = 1.0 

r = \ f " y2 

2 { ^Y(2y max -Ay) 


f n=0 for ^ 0 at x=l 
. dy 

L n= l for dx" 0 at x=l 


( 10 ) 


where 

X£,e (y) = x coordinate of wing leading edge at y = constant 
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<- y coordinate of the wing tip 

= range of x for which leading-edge adjustment is applied , 

=• range of y for which 1 ead i ng-edge .adj ustment is applied. 

The selection of the range for leading-edge adjustment is arbitrary depending , 
on how extensively one wants to impose the leading-edge property on the cal- 
culated velocity potential. For example, if one choses Ax = x^ e - X£ e and 
Ay = y max , the resulting adjustment becomes identical to that used in earlier 
versions of •sonic-box method computer program. Presently, both Ax and Ay are 
set to the length of the box side used .in computation. These relationships 
are depicted in figure 2. 


Tmax 

Ax 

Ay 


Pressure Coefficient 

The unsteady pressure coefficient for the harmonically oscillating thin 
wing is obtained from Bernoulli’s equation as 

Cp 0 = "2( c P 0x + ik<P 0 ), on the wing - 

(11) 

Cp Q =0, in the wake • 

and Cp(x,y,z,t) = Cp Q (x,y,z)e lkt . 

Generalized Aerodynamic Force Coefficient 

The generalized aerodynamic force coefficient as defined in reference 
6 is written as 

L ij = f j( (V + ikCp o^i f j dxd Y (12a) 

S ' 

= | jj (T g ~ + rtf 0 ) , fj dxdy. • (12b) 

S 


ORIGINAL' PMB IS 
OF POOR QUALITY 
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An integration by parts is performed to cast the integral in equation (12) 
into the following form: 


■I 


i = U(^) v= 


0 ' x=x 


te 


' (fx+ikfj^dxldy, x,y on, S, 


(13a) 


1 = j [ ^ V)x=x te " j mT^x +_ Mcf)9 0 dx]dy, x,y, on S. ( 1 3b) 

y * 

Equation (12b) is integrated in the transformed space using equation (13b). 


DATA FITTING METHODS 


The polynomial surface method of data fitting was used in the earlier 
versions of sonic-box computer program (refs. 1, 2, and 3) to fit the wing 
deflection mode, Mach number, velocity potential and pressure coefficient. 

In the present program, two types of data fitting are used. One is a 
two-dimensional cubic spline fitting and the other is a three-dimensional 
spline-surface fitting. The former is based on the concept of thin' elastic 
beam and the latter on the thin elastic plate. Since they use different forms 
of expression, the two-dimensional fitting can not be obtained directly from 
the three-dimensional fitting by setting the coefficient of the terms involv- 
ing the extra independent variable to zero as it is usually done when the 
polynomial fitting is used. 


Cubic Spl ine 


Method . - It has been demonstrated over the years that the draftsman's 

spline has distinct advantage over other methods of fairing curve through a 

set of predetermined points (e.g.‘, reference 7). The basic equation can be 

written as 
\ 


H = -f 1 ' [(xj-x 3 ) 4- d? (x-Xj.j) - dj] + i [(x-xj.,) -dj-(x-Xj-i)] 
+ Tj (x-Xj_-|), for xj_ -j <x<Xj- and j = 2,N 


(14) 
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+ Tj (x~ Xj - i ) , 
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where dj = Xj - xj_^ 

H = value of dependent variable at point x 
Hj = value of dependent variable at point xj 
Mj = coefficients of cubic spline 
N = total number of predetermined points 
Tj = (Hj - Hj _ i ) / d j 

xj = coordinate of the predetermined points where H is known 

The values of M j 1 s are computed with the N known values of H at xj by requiring 
the slope at each of these points to have a value computed from both sides of 
the point. This provides only (N-2) equations for N unknowns, thus the follow- 
ing two additional equations are assumed: 

M x = M N = 0. 

Use . - Cubic spline fitting is used in the computation of pressure co- 
efficient and generalized aerodynamic force coefficients according to equations 
(11) and (13), respectively. 

In the computation of pressure coefficient d i str ibut ion_on both physical 
and transformed wings, the known values of the velocity potential at the center 
of wing boxes, in the respective spaces, along each chordwise row are fitted 
with a cubic spline to calculate the chordwise derivative of the velocity 
potential at the box centers. The velocity potential distribution on the 
transformed wing is also fitted along each spanwise column to facilitate the 
interpolation of the velocity potential at the center of wing boxes along the 
same column on the physical wing. 

In the computation of aerodynamic force coefficients, the integration is 
performed in the space where the velocity potential distribution is calculated 
and in the sequence of chordwise and then spanwise directions. The cubic 
spline is used to linearly extrapolate both the velocity potential at the 
trailing edge of each chordwise row for the chordwise integration and the 
integrated value of the chordwise integration at the root chord for the 
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spanwise Integration, Both chordwise and spanwise integrations are performed 
with the implied cubic-spline fitting as shown in equation (16). 


Spl i ne-Surface 


Method . - In reference 8, it is shown that a spline-surface can fit a 
three-dimensional distribution of data better than the polynomial surface. 
The expression of the spline-surface can be written as: 


H 


N 

a 0 + aix + a 2 y + l Hir* anr* 
j = 1 J J 


0 


where 

r j = ^ x_x j) 2 + (y-yj) 2 

ao,a 1} a 2 ,hj = coefficients of spline-surface 

H = value of dependent variable at (x,y) 

N = total number of predetermined points 
(xj ,'y j ) = coordinates of the predetermined points where H is known. 


The valus of ao,a.i,a 2 , and hj's are computed with the N, known values of H at 
(xj,yj) plus the foTlowing three additional relations: 

N 

I hj =0 
j=1 J 

N 

l hj-xj = 0 

j = 1 


N 

l h j y : = o. 
J-i 


This yields (N+3) equations with (N+3). unknowns. 

The partial derivatives, of equation ( 1 5) with respect .to x, y can -be 
easily oblained as follows: 

lit 


VJT 



■^ H (x,y) 


g^HU.y) 



Use . - There are five sets of discrete points on the lifting surface 
which in general do not coincide: 


1. (x,y) on physical wing where Mach number is known 

2. (x,y) on physical wing where modal deflection is known 

3- (x,y) of box centers on physical wing, or 

(x,y) of box centers on transformed wing 
h. (x,y) on transformed wing corresponding to (x,y) on physical wing 

where modal deflection is known 

5. (x,y) on physical wing corresponding to box centers of transformed wing. 

In the process of solving equation (6) for the velocity potential, spline-, 
surface fitting is used to facilitate the interpolation of quantities and 
calculation of their gradients among these sets of points. The known local' 

Mach number distribution on physical wing is surface-fitted for the inter- 
polation of local Mach number on the physical wing at the box centers (for 
print-out purpose only),. at those points where the modal deflection is 
specified, and at those points corresponding to the box centers and the 
t ra i 1 i ng-edge points of transformed wing (for use in the evaluation of the 
integrals in eq. (13b)). The modal deflection points in the physical space 
are then transformed into transformed space according to equation (3) for 
nonzero-thickness wing. Note that the modal deflection f(x,y) is unchanged 
in the transformed space, i.e.. f(x,y) = f(x,y). The modal deflection, in 
the physical space for zero-th i ckness wing and in the transformed space for 
nonzero-thickness wing, is surface-fitted to facilitate the evaluation of the 
modal deflection and its chordwise derivative at box centers and the modal 
deflection at trailing edge as required in equation (8) for downwash and in 
equation (13) 'for force coefficient in the respective spaces. 
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When the modal deflection of the physical wing is expressed in a function- 
al form, .a number of points (i.e., box centers) on the transformed wing are 
selected for specifying the modal deflection and the corresponding points on 
the physical wing are then located with the aid of the surface-fitted Mach 
number. The modal deflect-ion at each of these points on the physical wing is 
then calculated with the known expression and this, in turn, is the value at 
each of the corresponding points on the transformed wing. 

PROGRAM FLOW 


The flow chart of the program SBOXR is presented in figure 1. The 
function of those controlling variables appeared in the flow chart is as 
f o 1 I ows : 


Variable 

Value 

Function 

NEW 

1 

indicates case without thickness effect 

2 

indicates case with thickness effect 

DA(26) 

0 

Indicates first frequency for a wing 

1 

indicates additional frequency for the 
same wing 



0 

Calculates cases with and without thickness 



effects 

DA (50) 

1 

Calculates case without thickness effects 



on 1 y 


2 

Calculates case with thickness effects 



on 1 y 

K 

M 

N 

<DA(28) 

Counter of the modes of deflection 

DA (28) 

5 MD 

Total number of modes to be considered 


16 




jFigure 1. - Flow chart of program SBOXR. 
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FUNCTION OF SUBROUTINES 


A brief description of the function of each subroutine is given in this 
section. Some equations are noted here for the purpose of identifying the 
steps in the sequence of obtaining the final results, i.e., the generalized 
aerodynamic force coefficients. It is suggested to refer to reference 1 for 
the "deri vat ions of these equations. The freestream flow is in the positive 
x direction and the apex of wing is at the origin. The spanwise direction 
is designated by y. Since the plane passing through the centerline chord 
and perpendicular to the wing planform is the plane of symmetry of wing 
geometry and motion, the input of wing geometry, deflection, and steady-flow 
Mach number distribution needs to be made only for the portion of the wing 
where x>0 and y>0 (see figure 3). It is assumed that the camber, twist, and 
mean angle of attack of the wing are all small; and the difference in local 
Mach number at the corresponding points of the upper and lower wing surfaces 
can be neglected. As the computer program is presently formulated, the motion 
is limited to symmetric modes. This limitation is imposed by the expression 
used in the present program to evaluate the Influence coefficient of upwash 
due to a unit doublet distribution. 

In the following, the variable name used in the computer program to 
represent a certain quantity mentioned in the description of a subroutine 
is indicated in the parentheses following such quantity, also, whenever a 
variable NEW appears as a subscript to another variable it is implied that 
the operation is applicable to either physical (NEW=1 , without thickness 
effect) or transformed (NEW=2, with thickness effect) wings. 

SBOXR - a controlling routine; also performs the final calculation of 
the generalized aerodynamic force coefficient for the full wing, as given in 
equation (12), and prints out the results. 

DATRD - data input routine; for a new wing, the data array (DA(k)) is 
cleared every time but for the same wing, the old value is used if a new 
value is not entered (blank datum will not affect the old value). 
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SHAPE - wing geometry routine; the main objectives of this routine are 

as follows: 

(1) approximates the wing planform of the physical or tranformed. wing with 
a grid of square boxes and calculates the number of boxes (ML(N£W,i,)) 
along each spanwise column ( XX ( i ) ) 

(2) finds the number sequence of the leading (MLT(NEW, 1 ,j) ) and trailing 
edge (MlT(NEW,2,j)) wing boxes along each chordwise row (YY(j)) 

(3) finds the number sequence of the first (MLC(NEW, 1 , i ) ) and the last 
(MLC('NEW,2, i) j wing boxes along each spanwise column (XX C i ) ) 

(4) computes the wing tip adjustment term { EDG (NEW , i ) ) along each spanwise 
column ( XX ( I ) ) 

(5) finds the Mach number at each wing box center on the physical wing 
(AMA(j,i,l)) and at those points on the physical wing corresponding 
to the box centers (AMA(j,i,2)) and the trailing edge (XTE(NEW,j), 
CTE(j)) along the chordwise now { YY ( j ) ) of the transformed wing 

(6) redistributes the number of leading edge (XEDG(k) , YEDG(k)) and trail- 
ing edge (XTDG(k) , YTDG(k)) segments by subdividing the corresponding 
input data (XEDGl(k), YEDGl(k) and XTDGl(k), YTDGl.(k)) for the physical 
wing so that the transformed wing can be better approximated 

(7) performs the coordinate, transformat ion 'of the wing with thickness effect 
with the aid of the known Mach number distribution (SFMX(k), $FMY(k), 
SFMH(k) , KSFM) on the physical wing, using the relations shown 1 2 3 in 
equation (3) j and 1 al so- checks the possible creation of an artificial 
wake in the transformed wing resulting from an uneven distortion of the 
wing tip 

-PLNFM - wing geometry routine called from subroutine SHAPE; the main 

functions' of this routine are as follows: 

(1) computes actual wing area in each box (AR( j , i ,NEW) ) 

(2) computes the x-coord i nates of the wing leading edge (XLE(NEW.j)) and. 
trailing edge (XTE (NEW, j ) ) along each chordwise row (YY(jj) 

(3) calculates the full wing planform area (AREA) of the physical wing 
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(4) finds the total number of boxes (ML(NEW,i)). including the wing as well 
as the wake along each spanwise column (XX(i)) 

(5) finds the number sequence of the first (MLC (NEW, 1 , i ) } and' the last 
(MLC(NEW,2, ?) ) wing, boxes along each spanwise column ( XX (i ) ) in which 
the actual wing planform occupies at least 50% of the box area. 

(6) identifies, along each spanwise column ( XX ( i ) ) , the number sequence of 
the first wake box (MLW(NEW,i)) of the wing-tip wake' (MLW(NEW, i } <0) and 
the last wake box of the wing-root wake (MLW(NEW, i ) >0) 

P0T2 - evaluates the integral in the expression relating the downwash 
and the velocity potent ial (equation (6)). The integral in equation (6) is 
integrated over the surface of each box B j i j i centered at ( x j i , y j i ) once for 
every frequency, and the result (A(N,j,i)) is used in the calculation of the 
velocity potential (S(N,j,i)) in both physical and transformed spaces for all 
modes under consideration at the same frequency. 

DRED - fits the wing deflection for each mode (M) of motion in both 
physical and transformed space (SFDX(fc,M), SFDY (£,M,NEW) , SFDH(£,M,NEW) , 
KSFD(M)) with a spline-surface of the form shown in equation (15) expressed 
in terms of the dimensionless coordinates. For the physical wing, the wing 
deflection in each mode can be given as input either in the form of the co- 
efficients (SFDH(£,M, 1)) of an expression H = ag-Fa^x (KSFD(M)=0) or deflection 
(SFDH (£ ,M , 1 ) ) at a number of predetermined points ,(SFDX(£,M), SFDY(£,M,1), 
KSFD(m)> 0). For the transformed wing, the wing deflection in a particular 
mode is obtained either directly from the input data points for the physical 
wing (KSFD(M)>0) or indirectly from a number of selected points (KSFD(M)>0) 
on the physical wing of known deflection of form h = ag+a^x. The deflection 
(SFDH (£,M, NEW) ) and the chordwise coordinate (SFDX( £. ,,M) ) of a point on the 
transformed wing are unaltered from those of the correspond i ng .po i nt on the 
physical wing, but its spanwise coordinate (SFDY (£, M, NEW) ) in the transformed 
space is modified by a factor of the Mach number (SFMX(k) , SFMY(k) , SFMH(k), 
KSFM) from that of the physical wing (see equations (3) and (8)). 

SLCT - called from subroutine DRED; selects a maximum of NB number of 
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points from the input data on the physical (NEW=1). or transformed (NEW=2) 
wing for the spline-surface fitting of the wing deflection. 

WVAL - calculates the complex (real N=1 , imaginary N=2) • downwash (S(N,j,i,)) 
at the center of each box (ij-th) for a particular mode (M) and reduced fre-. 
quency (CK.) using the known wing deflection spline-surface (SFDX(x,,M), SFDY(£, 
M,NEW) , SFDH (X, , M,NEW) , KSFD(M)) for the planar wing in physical (NEW=l) or 
transformed (NEW=2) space. 

BOXP - incorporates the leading edge adjustment, equation (9), to the 
computed ve loci ty potential (S(N,j,i)); computes pressure coefficient, equation 
( 1 1 ) , when a print out of its value is desired. To compute the pressure co- 
efficient for the ,wing with thickness, the velocity potential on the physical 
wing is found by transforming the known velocity potential on the transformed 
wing with the known coordinates and Mach number on the physical wing. 

BOXPO - called from subroutine BOXP, solves a set of simultaneous equa- 
tions, equation (6), relating the complex velocity potential (S(N,j,i)) and 
the known complex downwash of a particular mode (M) at the center of each box 
(ij-th) to obtain the velocity potential of the zero thickness wing in the 
physical and the transformed -spaces. Since the disturbance is considered not 
to travel upstream, -the calculation is carried out for one spanwise column at 
a time starting from the .most forward column. The number of upstream columns 
influencing a part i cul ar. col umn is controlled through the input (0A(49)). If 
a wake is encountered, the velocity potential in the wake along a chordwise 
row (YY(j)) is computed using equation (5) for various chordwise stations; the 
velocity potential at the trailing edge (XTE(NEW,j)) is interpolated from the 
velocity potential computed for the last wing box and the first wake box along 
the same chordwise row by considering the first wake box to be a wing box. 

FORCi - performs the integration of the generalized aerodynamic force co- 
efficient, equation (12a) or (12b), according to whether the wing is in the 
physical or transformed space. At each box center on the wing, the values of 
f and f x as required in equation (13) are obtained from the fitted spline- 
surface, equation (15) and T 0 is known. For wing with thickness, the Mach 
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number at each box center on the transformed wing (eq. ( 13 b)) is interpolated 
from the spline-surface representing Mach, number- distribution on the physical 
wing. The integration in x-di rection along each chordwise row (Y-Y(j.)).is 
carried out by computing the value of the integrand at a number of points 
along a particular row and then completing the integration with the quadrature 
expression shown in equation (16). The integration in -y-di rect ion is accomp- 
lished in the same manner. For the wing with thickness effect, the procedure 
of evaluating the integral is the same except the involvement of the Mach 
number at those- points on the physical wing corresponding to the box centers 
( AMA ( j , i , 2) ) and trailing edge ( CTE ( j ) ) on the transformed wing. 


MRED - fits the mean steady-state Mach number distribution (SFMX(k) , 
SFMY(k), SFMH(k), KSFM) over the phys ical wing with a spline-surface, equation 
(15), expressed in terms of the dimensionless coordinates. The input for the. 
Mach number can be either in the form of the coefficients of an expression 
Mach number = ao+a^x (KSFM=0) , or as tabulation (SFMH(k)) of the Mach number* 
or the pressure coefficient at a number of predetermined points (SFMX(k) , 
SFriY(k), KSFMs-O) on the physical wing. 


1NTGL - called from subroutine F0RC1; carries out the evaluation of the 
integral in equation ( 13 ) by, first, fitting the values of the integrand, at 
a set of predetermined points with a cubic spline, equation (1*0, and then 
by using the following quadrature expression: 


1 -il d j * (Hj-1-Hj)]. 

J ^ 


( 16 ) 


Equation (l6) sums the sectional ly integrated area under a' cubic spline 
(equation (1*0); which fits N values of the integrand, between two adjacent 
points, j-1 and j. 

TR1DI - called from subroutine INTGL ; solves for the roots o.f .a set of 
equations whose coefficient matrix is. a tri-diagonal matri-x.. 


SPLN1 - provides a cubic-spline fit for a. variable -at a given set of 
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points, from.which the interpolated values of ' the dependent variable and 
its derivative may be obtained using subroutine SPLN2; also inserts an addi- 
tional point bet .een every successive pai.r of input points in the original set 
-to increase. the number of points for cubic-spline f i tting. and- performs the 
specified number of smoothing iterations using subroutine SMOOTH. 

SPLN2 - called from subroutine SPLN1; interpolates the value of the 
dependent variable and calculates its derivative at the set of desired 
locations. 

SP1SET - called from subroutine SPLN1; computes the slope of the dependent 
variable using the cubic spline method at those points where the values of the 
dependent variable are known. 

SMOOTH - called from subroutine SPLN1; replaces the value of the de- 
pendent variable at each of the predetermined points with the new value 
computed by passing through a five-point least-squares cubic before the slope 
at these points is computed. 

CHLSKY - called from subroutine SMOOTH; solves for the roots of a set of 
simultaneous equations whose coefficient matrix is a symmetric matrix. 

SURF1 - spline-surface fitting routine; with the known values of the 
dependent variable at a number of points of known coordinates, this routine 
finds the coefficients of the spline-surface associated with those known 
points of the form shown in equation (15)* 

SURF2 - interpolating routine; finds the value of the dependent variable 
and its partial derivatives at a specified point based on the 'fitted spline- 
surface. 

C 1 N - calculates sine and cosine integrals. 

MSIMER - finds solutions of the simultaneous real number equations. 

MSIMEC - finds solutions of the simultaneous complex number equations. - 
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PROGRAM DIMENSIONS 


The present computer program is dimensioned to handle a. maximum of 30 (MB) 
boxes either in chordwise or spanwise direction in approximating one-half of 
the wing planform. The maximum numbers of leading and trailing edge segments 
are, respectively 7 (NEDGI-1 ) and 2 (NTDGI-1). It can handle up to 3 (MD) wing 
deflection mode shapes. The maximum number of points used in the spline- 
surface fitting is 100 (NB) . These limitations can easily be increased by 
-changing the dimensions of the corresponding variables in the computer program. 

The maximum number of input points to describe the wing deflection and 
the thickness effect expressed in terms of the steady-state, either the Mach 
number distribution or pressure coefficient distribution is presently dimen- 
sioned 100 ( N B ) . This number can be increased to 128 by activating the unused 
portion of the input data array, DA, from locations 501 to 700 and by changing 
the value of KP, the starting location for the thickness effect input points, 
from 701 to 6 1 3 in subroutine MRED and the value of NB from 100 to 128 in 
program SB0XR. The dimensions of the variables used in the spline-surface 
fitting need to be changed according to the new value of NB. 

RESULTS 


Sample calculations were made using the present program as well as other 
versions of sonic-box computer program. It was found necessary to use a unit- 
box increment instead of a multi-box Increment In' the convergence test due to 
the inherent property of the box method. Since no smoothing is applied to 
any computed values in the present program, the pressure coefficient calcu- 
lated should be used with caution. In order to use it, one needs to put the 
pressure coefficient distribution through a smoothing process such as the 
built-in smoothing routine which was not utilized, however, for the results 
presented in this section. The pressure coefficient is obtained by differen- 
tiation of a set of numerical values whereas the generalized aerodynamic force 
coefficient is obtained by integration. Since integration itself is a smooth- 
ing process, the resulting generalized aerodynamic force coefficient is 
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considered to be acceptable within the bounds. of the accuracy of the. numerical, 
techniques and the adequacy of the sonic-box method. , . . 

Wings considered in the sample ca 1 cul at ions- .are .del ta , cropped- del ta , 
arrow and -rectangular wings and are given separately below. 

Aspect Ratio 1.5 Delta Wing 

Zero Thickness . - Shown in figure k -is the convergence with respect to 
the -number of boxes along the root chord of the generalized force coefficients 
due to plunge (mode 1) and pitch about its apex (mode 2) at a reduced fre- 
quency 0.2 for an aspect ratio 1.5 delta wing. The results obtained from the 
present and the other sonic-box computer programs given in references 1 and 2 
are plotted in the -same figure for comparison. The program of reference 3 
generated the same results as that of reference 1. The results of the present 
program appear to show more clearly the converging trend at a fewer number 
of boxes used in the computation than the other sonic-box programs. 

The computed and the smoothed velocity potentials along a chordwise row, 
y = 0.0166667 from root chord, of the aspect ratio 1.5 delta wing due 
to plunge at a reduced frequency 0.2 from reference 3 using 30 boxes along 
the root chord are plotted in figure 5. The results from the programs 
of references 1 and 2 show the similar pattern of velocity potential 
distribution. In these programs, a smoothing of the computed velocity po- 
tential distribution is applied to obtain the distribution for the subsequent 
computation of the pressure and force coefficients. The velocity potential 
distribution over the same wing at the same conditions using the present 
program is shown in figure. 6. in the present illustration, no smoothing is 
used; thus, the velocity potentials used in the computation di ffer from the 
computed values only near the leading edge where the leading edge adjustment 
- i s appl i ed . 

Nonzero Thickness . - Figure 7 shows the variation of, the force coeffi- 
cients due to plunge (mode 1) and pitch about its apex (mode 2) as function 
of the -reduced frequency for the aspect-ratio 1.5 delta wing computed with 
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the present program using 30 boxes along the root chord. The results with 
thickness were calculated with the assumption that the wing has an elliptic 
cross-section in a plane perpendicular to the chordwise axis and the steady- 
state pressure coefficient distribution used as input was obtained from the 
method given in reference 3 . The distribution of Mach number at the box 
centers interpolated from the fitted spline-surface is plotted in figure 8. 

The Mach number distribution in spanwise direction for this wing is supposed 
to be constant but the interpolated value has deteriorated near the trailing 
edge. Since the error in magnitude is not large, the effect on the integrated 
results shown in figure 7 may not be very serious. As it was observed in 
reference 4, the present results also show that the thickness effect on the 
force coefficients of a delta wing is not very large. However, the effect on 
flutter speed can be significant (ref. 10). 

Aspect Ratio 2.0 and 4.0 Delta Wings with Zero Thickness 

Convergence of the force coefficient due to pitch, f = -0.00617 + 

0.017446 x (b = 0.5 ft.), of aspect ratio 2.0 delta wing at a reduced fre- 
quency 0.106 is shown in figure 9 and that due to plunge (mode 1) and pitch 
about its apex (mode 2) of aspect ratio 4.0 delta wing at a reduced frequency 
0.1 is shown in figure 10. In figure 9, the present results are compared 
with that given in reference 5. The results of references 1 and 3 are in- 
cluded in figure 10 for comparison. According to the calculations performed 
for the delta wings, the results of the present program show a better and 
faster converging trend than that of the other versions of the sonic-box 
method computer program. 

The computed velocity potential distribution without smoothing, due to 
plunge (f = 1.0) at a reduced frequency 0.1 of aspect ratio 4.0 delta wing 
from the program of reference 1 is plotted in figure 11, and that from the 
present program is plotted in figure 12. The numerical values of these two 
sets of results are not very different. The fluctuation of the values of 
velocity potential along the chordwise row as observed in figure 5 is not 
observed in figure 11. This is due to the difference in the box arrangement 
along the swept leading edge of the delta wings. 
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The convergence plots of the force coefficients of a delta wing as func- 
tion of the number of boxes along the root chord always show fluctuations; it 
does not increas - or decrease monoton i ca 1 1 y toward the converged values. This 
is an inherent property of the box-method. 

As shown in figure 13, the box arrangement along the swept leading edges 
can only approximate that of the actual wing. The number of boxes in a span- 
wise column does not increase uniformly from the nose to the tail of the delta 
wing. For- example, the aspect ratio 1.5 delta wing has 3 spanwise columns 
with the same number of boxes, followed by another 3 columns with the same 
number of boxes and followed by 2 columns with the same number of boxes and 
this sequence, 3“3"2, is repeated toward the trailing- edge of the wing. .For 
the aspect ratio 2.0 and 4.0 delta wings, the sequences are, respectively, 

2-2 and 1-1 . 

From the plots shown in figure 5, one notices that the fluctuations of 
the computed velocity potential along a chordwise row is caused by the se- 
quence of increase of the number of boxes in the spanwise columns from the 
nose toward the tail of the wing. The plots in figure 5 clearly show that 
the values of the computed velocity potential along a chordwise row increase 
whenever the number of boxes in a spanwise column increases and remains the 
same or decreases whenever the number of boxes remains unchanged. This is 
the way the fluctuation is formed. 

The fluctuation of the computed velocity potential along the chordwise 
row of the aspect ratio 4.0 delta wing is not observed. This is due to the 
fact that the number of boxes in the spanwise column increases by one un i form] y 
from the leading edge to -the trailing edge, and the magni-tude of the computed 
velocity potential only increases in value. 

The fluctuation of the force coefficient as function of the number of 
boxes along the root chord shown in figure 9 is probably caused by the num- 
ber of boxes in the spanwise column at the trailing edge of the aspect ratio 
2.0 delta wing. The box number sequence for this wing is 2-2 and the last 
spanwise column may end at the middle or end of the sequence depending on 
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whether an even or an odd number of boxes along the root chord is chosen, 
the example shown in figure 5 ends at the middle of the sequence. 

The similar data from the present program are shown in figures 6 and 12. 
Except for the leading edge adjustment,, the .results presented in these figures 
are the computed values unsmoothed. The fluctuation of the computed velocity 
potential along the chordwise row is greatly reduced, this is achieved by 
including the contributions from the partial boxes along the leading edges. 

Cropped Delta Wing with Zero Thickness 

Shown in figure 14 is the convergence of the -force coefficients due to 
plunge (mode i), f = 1.0, and pitch (mode 2), f = -0.4034 + 1.832 x, (b = 

0.5833 ft.), at' reduced frequency 0.095 with respect to the number of boxes 
along the root chord of the aspect ratio 1.5 cropped delta wing with taper 
ratio 0.143. The results computed from the programs of references 1 and 2 
are also included for comparison. It is again observed that the results of 
the present program shbw a better converging trend than the other versions 
of the sonic-box computer program. 

Arrow Wing with Zero Thickness 

The convergence Of the force coefficients due to plunge (mode 1) and 
pitch about its apex (mode 2) at a reduced frequency 0.5 of a pointed tip 
wing with aspect ratio 4.0 and leading edge sweep angle 56.33 degrees is 
shown in figure 1-5. The results computed from the present program and that 
of reference 2 are included in the figure. The convergence rate as well as 
the converged values of these two programs are about the same. 

> 

Distributions of the velocity potential due to plunge on the wing and 
in the wake calculated from the present program and -that of reference 2 are, 
respectively plotted in f.igures 16 and 17- T^e reduced frequency is 0.5 and 
the number of boxes used in the chordwise direction is 30- .The leading edge 
adjustment is included in the resul ts . shown in fiugre 16, but no smoothing-. 

The results computed from the program in reference 2 shown in figure 17 are 
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the smoothed values and the computed values along 1st, 5th, 10th and 15th 
chordwise rows. These computed velocity potentials exhibit the similar 
fluctuation in t'ij chordwise direction as observed in' the case of delta wing. 

_ The velocity potential at the trailing edge of the arrow wing is inter- 
polated from the velocity potentials computed at the centers of the last wing 
box and the first wake box along the same chordwise row. This wake box is 
assumed to be a wing box in the preliminary computation to facilitate the 
evaluation of the velocity potential at the trailing edge; this box is treat- 
ed as a wake box thereafter. In reference 2, the velocity potential ?n the 
wake along a chordwise row is computed from the velocity potential at the 
trailing edge of the wing on the same row before the -smoothing is applied. 
After the smoothing, the velocity potential at the trailing edge may not re- 
main the same as the one originally used in the computation of the velocity 
potential in the wake. This, in turn, raises the question of the correctness 
of the velocity potential in the wake which was used in the computation of 
the velocity potential on the portion of the wing downstream of the wake in 
question. The problem area of the computer program in reference 2 described 
here is depicted in figure 17* In the present program, this does not become 
a problem because the computed velocity potentials in the chordwise direction 
are fairly smooth and no additional smoothing is applied. 

Rectangular Wing 

Zero Thickness . - The results for a rectangular wing of aspect ratio 2.0 
performing plunge (mode l) and pitch about its leading edge (mode 2) at re- 
duced frequency 0.6 are presented in figures 1 8 , 19, and 20. In figure 18, 
the convergence of the different sonic-box programs is shown. The convergence 
of the present program again showed a better trend than the other versions of 
sonic-box method computer program. Since the calculations were made only up 
to 20 boxes along the root chord, the converged Values of the other programs 
were not yet attained. The difference in the values of force coefficients 
of the present method and that of reference 3 might have been due to the 
different data fitting techniques used in the programs. The different ways 
of treating the wing leading edge in these two programs do not have effect on 
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the results for this wing because there is no partial box involved in the 
computation. 

Nonzero Thickness . - The convergence of the present program and that of 
reference 3 for both zero' and nonzero thickness is shown in figure 19.. .This 
rectangular wing has a biconvex airfoil section with a thickness ratio 0.0521 
and the steady-state Mach number distribution on the wing was' calculated with 
the method of reference 9. The thickness has greater effect on the rectangular 
wing than on the delta wing (fig. 7). This may be caused by the fact that the 
Mach number on the rectangular wing (fig. 21) deviates from unity more than 
that on the delta wing. For the nonzero thickness-, in addition to. the differ- 
ences in data fitting technique and handling of wing leading edge, the method 
of computing the force coefficients is also different. In the present method, 
the force-coefficient integration is performed in the transformed space, 
whereas the method in reference 3 integrates in the physical space. 

Variation of the force coefficients with- respect to reduced frequency is 
plotted in figure 20. Both zero and nonzero thickness using the present 
program and that of reference 3 with 20 boxes along the root chord are shown. 
The results of reference 6 for the zero thickness wing are included in. figures 
20(c) and '(d) for comparison. Agreement between the present results and that 
of reference 6 is better in the medium frequency range. This may be attri- 
buted to the different assumptions imposed on the parameters involving reduced 
frequency. The results of reference 6 are based on a theory for large values 
(greater than 0.8) of the product of aspect ratio and the square root of re- 
duced- frequency, whereas, the present program is based on the assumption of 
small values of the product of box-side length and reduced frequency. With 
the sonic-box method, one needs to use more boxes to represent the wing as 
the frequency increases. 

ORIGIN AJj 

•REMARKS AND RECOMMENDATIONS OF 

A sonic-box method computer program using the concept of local linear- * 
ization to account for the effects of wing thickness, in unsteady sonic flow... 





is presented. This program uses a coordinate transformation to reduce the 
small perturbation, unsteady transonic veloc i ty. potent ial- equation, which has 
variable coefficients, to a linear equation with constant coefficients and 
solves this boundary-value problem in. the transformed space by the sonic-box 
method by regarding the mean local Mach number as a parameter. 

Accuracy of this approach deteriorates when the variation of Mach number 
on the physical wing is large. in fact, the method fails when the lateral 
variation of Mach number becomes large enough to cause multivalued transforma- 
tion which implies wing-surface fold-over in the transformed space. Moreover, 
difficulties can also occur when artifical "wakes" created by the transforma- 
tion impinge upon downstream portion of the lifting surface. Specifically, 
these "wakes" exist if yM along leading edge (or tip) at a particular 
chordwise location is greater than that at another location further downstream 
Existence of these artificial "wakes" causes no difficulty, however, unless 
a portion of the transformed wing further downstream protrudes into these 
"wakes". Therefore, the user of this program should examine the geometry 
of the transformed wing to determine if these conditions exist. 

This program includes the contribution of the partial boxes along the 
wing leading edge in the calculation of velocity potential, thus, allows' 
one to use a smaller number of boxes to represent the wing than other versions 
of sonic-box computer program do. The earlier polynomial -surface data-fitting 
technique has been replaced by the cubic spline and the spline-surface to 
improve the accuracy. 

Since the box method itself is numerical in nature,' the distribution of 
the calculated values is not always smooth, that is, some degree of data 
oscillation is unavoidable. This means that some smoothing of the computed 
data will be desirable not only because of the oscillation of the computed 
data but also because of the very nature of cubic spline and the spline- 
surface techniques. Even though a smoothing scheme has been included in 
the program, it was not utilized in the samples presented.. The- adequacy of 
this smoothing scheme has not been extensively tested and the differentiation 
of the computed values has been avoided in the calculation of generalized 
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aerodynamic force coefficients. It is quite possible to- lose the desired 
feature of the interpolating surface after the smoothing unless a proper 
control of the -smoothing process is employed. 

In the program, the vibration mode shape and Mach number are represented 
by spline-surface. The accuracy of the interpolating surface can be affected 
by the choice of the data points used to generate it. To facilitate an 
assessment of the accuracy of the fitting technique, nonzero thickness and 
flexible modes were not emphasized. Instead, zero thickness and rigid body 
modes were studied more extensively in the samples presented. Emphasfs was 
placed on the consistency of the converging trend with the present program. 
According to the results presented, the thickness effect on the rectangular 
wing is quite significant, as much as 30 %, but smaller for the delta wing. 
This is attributed to the d i f ference in the Mach number variation for the two 
wings . 


Since the accuracy of the present program depends on the data-fltting 
techniques used in the program, it would be most desirable to modify or 
change the fitting scheme so that the flexible mode shape and the Mach number 
distribution can be better represented and interpolated during-the computation. 

One method to accomplish this is to generate accurate spline-surface for 
Mach number and each mode' shape from a separate program and use these as part 
of the input. Another method is to input Mach number and deflection of every 
mode at each box center on the physical wing without fitting the data with 
surfaces and use the cubic spline to facilitate the interpolation required 
during the computation. Ei.ther of these two methods undoubtly will increase 
the computation time, but the additional cost may be compensated wfth the use 
of a smaller number of boxes to represent the wing. 
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Figure 5- - Velocity potential along chordwise row at y = 0.01 66667 due to 
plunge for delta wing of aspect ratio 1.5 at k = 0.2 with 30 
boxes along root chord. Zero thickness. 
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Force coefficients due to plunge and pitch for delta wing of 
aspect ratio 1.5 with 30 boxes along root chord. 
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Figure 7. —Continued. 
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Figure 8. - Interpolated Mach number for t = 0.10 
delta wing of aspect ratio 1-5. 









-Q.64Q«- 


Lll 


- 0.630 

- 0.620 

-0.610 

- 0.600 

-0.590 

-0.580 

-0.570 

- 0.560 



Present method 

* Reference 3 . 

* • Reference 1 


-0.080 p 

-0.070 - 


-0.060 L 



_1 t 1 1 L 1 1 1 I 1 1 1 L_j 1 — 1 — 1 

14 16 18 20 '22 24 26 28 30 

Number of boxes along root chord 


(a) Lift due to plunge. 

Figure 10. - Convergence of force coefficients due to plunge and pitch 

for delta winq of aspect ratio 4.0 at k = 0.1. Zero thickness. 
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Figure 10. - -Concluded. 
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Figure 11. - Velocity potential along chordwise rows due to p.lunge for 
delta wing of aspect ratio 4.0 at k = 0.1 with 20 boxes 
along root chord — reference 1. Zero thickness. 



SSSp-tteK 

op PO0K 


Imaginary 









ORIGINAL' 

OF POOR QUAT-tf^ 


~N— /b 2 i, deg 


Present method 
Reference 1 
Reference 2 


/ V ' \ / . \ / 

- 

f/\V • w ; 


^ A 

T - pr~ 

\ i \ 

'■ ' / ; •. ¥ 
\ v/ v 


Number of boxes 
along root chord 


Number of boxes 
along root chord 


Figure 1-4 . - Convergence of force coefficients due to plunge (mode 1) and 
pitch (mode 2) for cropped delta wing of aspect ratio 1.5 and 
taper ratio 0.143 at k= 0.095- Zero thickness. 
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Figure 15. - Convergence -of force coefficients due to plunge (mode 1) and pitch 
(mode 2) for arrow wing of aspect ratio A.-O and' leading edge sweep 
angle 56 • 33 degrees at k = 0.5. Zero thickness.' 
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APPENDIX 


COMPUTER PROGRAM 
Input Guide 

Data are input through the subroutine DATRD us ing the one d imens ionai array 
DA wi th a- s i ze of 1005. The allowable maximum number for some of the input data 
as indicated below may be changed if the d imens ion of the corresponding storage 
array and computational operat ions are al so changed accordingly. Subroutine 
DATRD initializes DA ( 1 ) through DA (22) to bl ank, the weighting factors in 2-<. 1 04) , 

DA(108), ,DA(500) to 1.0, and the remaining portion of the DA array to 2.0. 

Consequently, these are the default values. The layout of the srrav C->.k) as 
it is presently used is as follows: 

1-7: Title 

8-12: Not used 

13“ 19: Mode title 

20-22: 'Not used 

23: Frequency, (cycle/sec) 

24: Overall length of wing in streamwise direction, (ft or meter) 

25: Speed of sound of the freestream, (ft/sec or meter/sec) 

26: (0) - -indicates the frequency is the first one for a new v.ing 
(l) - indicates the frequency is the add i.t ionai one for t-he 
sane wing 

27: Number of boxes in streamwise direction (maximum 30) 

28: Numbe r of deflection modes (maximum 3) 

29: Number (n) of segments of leading edge per semi span to be given, 

excluding segment from origin to yo (m max = 7) 

30-44:. Coordinates of points on the leading edge, (ft or meter) 

(in sequence of V ; ,xi ,y ! , x 2 ,Y2 > — >*n> Ym) > n rr.ax_"„ J 
45: Number (n) of segments of trailing edge per semispan to be given, 

(default: unswept trailing edge), n max = 2 

46-48: Coordinates of points on the trailing edge, (ft or meter) 

(in sequence of xj,y 1 ,x 1 for n=2, 
or Xr> (only) for n=l , 
no input for n=0; 



49 : 


50: 


51: 


52-53: 

54-70: 

71-72: 

73-95: 

96: 


97: 

98 : 

99 - 100 : 

101 - 500 : 

501-700: 

701 - 1000 : 


-ORIGINAL* PAGB 1$ 
OP POOR quality: 

last trailing edge point conincides with the last leading edge 
point and is set internally) 

Numbe- of boxes' allowed for upstream influence (if .this location 
is left blank or assigned a zero, it'will assume DA(49)=DA(27) 
and in no case DA(4g) >DA(27) is allowed) 

(0) - indicates to calculate cases with and without thickness 

effect 

(1) - indicates to calculate case without thickness effect only 

(2) - indicates to calculate case with thickness effect only 
Indicator to suppress calculation of potential for a mode 

(0) - no suppression 
(1-) - suppression 

Coefficients of the deflection polynomial (in the sequence of 
ag and a^) • • 

Not used* 

Coefficients of the Mach number distribution polynomial (in 
the sequence of ag and a^) 

Not used* 

Indicator of the type of wing- thickness effect input 

(1) - pressure coefficient 

(2) - Mach number 

Number of points at which pressure coefficient or Mach number 
to be given 

Number of points on which deflections to be given 
Not used* 

Deflection data for a maximum of 100 points (in the sequence 
of x, y, def 1 ect ion’ and weighting factor) 

Not used** 

Pressure coefficient or Mach number data for a maximum of 100 
points (in the sequence of x, y and pressure coefficient or 
Mach number) 


The remaining part of DA array is used for -the. control of -intermediate 
results print out. When the latter is desired, a non-zero positive integer 



number should be entered at locations in the DA array corresponding to the 
information from one particular subroutine is needed. 

1001: CBA: for wing deflection (DRED) 

A=l, for NEW=1 \ • 

f spl me-surface fitted results 
B=1 , for NEW=2 J 

C=1 - coordinates of the selected points used in 

spline-surface fitting 

1002: BA: for wing upwash (WVAL) 

A=1, for NEW=1 'l 

> upwash 

B=1 , for NEW=2 J 

1003: FEDCBA for velocity potential (B0XP and B0XP0) 

A=1, for NEW=1 'i 

?■ velocity potential 

B=1 , for NEW=2 J 

C=1, for NEW=1 ) . . 

f influence coefficient and solution matrices 
D=1 , for NEW=2 J 

E=1 , for NEW=1 1 

> pressure coefficient 
F=1 , for NEW=2 J 

1004: A: for Mach number (MRED) 

A=1 - spline-surface fitted results 

1005: DCBA: for wing shape (SHAPE and PLNFM) 

A=1 , for NEW=1 "1 distributions of box, box area , ' 1 ead ing and 

B=1 , for NEW=2 / trailing edges 

0=1, - Mach number at box centers 

D-1, for NEW=2 - redistributed leading and trailing edge segments 

The format of the input data card is (A1 , A5, i6, 6A10, A8) . The first 
field i.s for the control of clearing the data array, DA, for a new wing '(+) 
and the control to indicate the end of the set of data (-). The second field 
is the indicator for the type of data, either numeric (blank) or alphameric 
(ALPHA). The third field is the designator for the relative location in the 
data array of the first number to follow in the fourth field. if this field, 
is left blank, or a zero is entered, the execution will be terminated. The 
fourth and fifth fields are for five consecutive input data each occupying 12 
columns plus 8 blank columns at the end. All the fixed point numbers are 
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right-adjusted and the decimal point for the floating point number must be 
included. If an input datum is left blank, no change at the storage location 
for that particu'ar datum in the data array will occur .unless the set of the . 
input data is for a new wing. 

Those storages currently not used in array DA marked with * are reserved 
for future improvement in the method used for the functional form of data 
input. Those marked with ** are reserved for the case where a large number 
of data input points for deflection or Mach number is required. 


Sample Case 

A typical data deck set-up and output for an aspect ratio 1.5 delta wing 
having an elliptic lateral cross-section with 10% thickness ratio performing 
plunge and pitch about its apex are given below. 


Input 

The input format is (A1 , A5, 16, 6A10, A8) . 


Card 1 : 
.Card 2: 
Card 3‘ 

Card b: 

Card 5: 

. Card 6: 


Card 7: 
Card 8: 
Card 9: 


title of the case under consideration, 
title of the first mode of deflection. 

first frequency (cycle/sec), centerline chord length (ft), 
reference velocity (ft/sec). 

number of boxes along the centerline chord, number of deflection 
modes, number of total leading, edge segments of the wing, 
spanwise coordinate (ft) of the first section of' the leading edge, 
chordwise and spanwise coordinates (ft) of the next section (the 
sequence is y 0 , xj, yi — e.g., see figure 3) • 
first mode of deflection f = 1.0 

the 11 - 11 sign indicates the end of the' group of data "cards 
to be read at this stage, 
rt 1 1 1 e of the second mode, 
second mode of deflection f = O.lx. 

identification of the type of input regarding the wing thickness 
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effect (Mach number for this case), number of points on the wing 
this information to be given. 

Cards 10 to 69: 

chordwise and spanwise coordinates (ft) of a point on the wing, and 
the Mach number at this point. 

the " - " sign on the last card indicates the end of the group of 
data cards to be read at this stage. 

Cards 70 and 71 : 

additional frequencies for the same wing, one card is read in at 
one time. 

Card 72: blank card to make an exit from the computer. 
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IS I 2 2NAI; s&es 


The card images le as follows': 


V 2 3 4 S 6 ’ ~ “»--- 

123o55789ol23456789Ol234S67A9()l?34567a9ol234S6789Ol234567a9Ol234567895i23-*.567890 


♦ALPHA 1 ASPECT PATIO 1,5 DELTA WING <TAU*0,10) 
ALPHA 13PLUNGE 


23 

O ,159154941 

10.0 

looo. o 

27 

30 

2 

1 

30 

0.0 

10.0 

3.75 

■52 

1,.0 




ALPHA 13PITCH A90UT ROOT LEADING EDGE X*0.o 


52 

0,0 

0.1 


96 

1 

60 


701 

0.1 


0.140232 

704 

0.3 


9. 135125 

707 

0.5 


0.132175 

710 

0.7 


0.129943 

713 

0.9 


0.126094 

716 

1.3 


0.124980 

719 

1.7 


0.122340 

722 

2.1 


0.119963 

725 

2.5 


0.117746 

728 

2.9 


0.115627 

731 

3.3 


0.113564 

73<* 

3.9 


O.llOSio 

737 

4.5 


0.107428 

74C 

5.1 


0.104243 

743 

5.7 


0.100872 

746 

6.3 


0.097215 

749 

6.9 


0,093129 

752 

7.5 


0,080395 

755 

7.9 


0.084697 

758 

8.3 


0.080342 

761 

0,7 


0 , 07497c 

764 

9.1 


0.067325 

767 

9.3 


0.063062 

770 

9.S 


0.056823 

773 

9.7 


0.047661 

776 

9.9 


0. 029651 

7,79 

1.7 

0,5 

0.122340 

752 

3.3 

1.0 

0.113564 

755 

4.5 

1.0 

0,107428 

7RO 

5.7 

1.0 

0.1 00372 

791 

6.9 

1.0 

0.093129 

794* 

8.3 

l.o 

0.080342 

797 

9.5 

1." 

0.056823 

800 

4.5 

' 1.5 

C. 107*28 

303'. 

6.3 

2.0 

0.097215 

806 

7.5 

2.0 

0.088395 

809 

8.7 

2.0 

J. 0749 70 

812 

9.7 

2.0 

0.047661 

815 

o.3 

2.5 

0.060342 

B18 

9.7 

3.0 

0.047661 

921 

2.5 

0.5 

0.117746 

524 

3.9 

0.5 

o.nosio 

827 

5.1 

0.5 

0.104243 

830 

6.3 

0.5 

C. 097215 

633 

7.5 

0.5 

0.088395 

836 

9.1 

0.5 

0.067625 

B39 

9.9 

0.5 

: .029651 

84? 

5.1 

1 .5 

0.104243 

34$ 

6.3 

1,5 

C.C9721S 

348 

7.9 

1.5 

:.:£4697 

651 

9.1 

1.5 

:,C67S2S 

954 

9.9 

1.5 

0.029651 

657 

•5.7 

2.0 

C .100372 

56C 

6.9 

2.5 

0.093129 

863 

9.5 

2,5 

* • 

866 

9.9 

2.5 

<}. 02*651 

869 

6.3 

3.0 

0. C$02462 

872 

9.1 

3.0 

0.067325 

875 

9.5 

3.5 

0.056823 

37fl 

23 

23 

9.9 

1 .591549407 
6.366197626 

3.5 

0.029651 


1 

2 

3 

A 

5 

6 

7 

8 
9 

10 

11 

12 

13 

l* 

15 

16 
17 
l? 

19 

20 
21 
22 

23 

24 

25 

26 
27 
26 

29 

30 

31 

32 
‘ 33 

36 
35 
36- 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 
40 
4 ? 

50 

51 

52 

53 

54 

55 

56 

57 
56 

' 59 
60 

- -61 
62 

63 

64 

65 

66 

67 

68 
69 
73 
71 
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Output 



ASPECT PATIO 1,5 DELTA *1*0 

(TAU=0ol0) 




30 ROES ALONG ROOT CHORD 

ROOT 

CHORD LENGTH * 

10,00 ET 


REDUCED FREQUENCY = .01 G 

EPEE 

STREAM VELOCITY = 

1000,00 ET/SEC 


FREQUENCY = T.592E-01 CYCLE/SEC 




"•■ODE No. 1 PLUNGE 





>'OOE NO. 2 DTTCH APOU t ROOT LEADING 

EDGE X=o.O ' 



GENERALIZED FORCES (NO THTC<NESS FEEFcT) 



vnpr $ 

PPE?. 0€FU # ' RFAL ^APT 

ImaG part 

ASS. VALUE 

PHASE ANGLE 

i 

1 2.1Sll?E-0f 

-2.4278.3f-0? 

2.427A3E-02 

-89.9949 

1 

2. -1.S9 1 36f -07 

-1.611 59F-02 

1 .61159E-02 

-90.0006 


generalized pqpces (n.'' thickness effect) 



-odes 

. PRES. DE r L. »EaL "APT 

IMA-3 RAPT 

A3S, VALUE 

phase angle 

2 

1 -2.4?7S8E*00 

-2.44596E-02 

2.42200E+00 

-179.4226 

2 • 

2 -l.M162E*oo 

-l.ei749F-02 

1.61173E*00 

-179.3553 


NODE No . 1 




- 

M ODF NO. .2 





GENERALIZED fopCES fWTT-N thickness Effect) 



-odes 

PRES. OEPL, . ’ Rf AL p A=T 

I«A3 PART 

ASS. VALUE 

PHASE ANGLE 

I 

1 -1.31S30E-OS 

-2.427 Z9F-.12 

2.42779E-C2 

-90.0 3 I'D 

1 

? -1 . 1 18 14E-0S 

-1.61155E-02 

1.61156F-02 ' 

-90.0393 


GENERALIZED FORCES (rfl-TH tficknESR EFFECT) 



moT; 

PRES. 

ES 

OE r L. REAL PART 

T«Ao PART 

ABS. VALUE 

phase angle 

? 

1 -2.A?95C£ + n f) 

-2.293 7 8E-02 

2,427R6E*00 

-179.4587 

_2 _ 

2 -I.81l60E*00 

-1 .70333F-02 

1 .61 169E*0A 

-179.3945 
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oaasaHtf 

0F PO0B 


tfjAiarv 


ASPECT PATIO 1,5 Of |_T A WI*G (TAU=0.10) 

30 BOXES ALONG PO' , T CHORD ROOT CHORD i.t>15TK a 10,00 FT 

REDUCED FREQUENCY = .100 FPff STREAM i^U_QCITY a 1000.00 FT/SEC 

FREQUENCY = 1 .59EE+0O CYCLE/SEC 

MODE NO. 1 
MODE NT. 2 

GENERALIZED FOpCES (NO THICKNESS FFFfCT) 


MOOES 

PRES. DEFL. 

REAL FART 

ImaG part 

AS 5. aALUE 

PHASE ANGLE 

1 1 

3.63985E-0^ 

-2.42141E-OI 

2.A2l4tE-flI 

-89.9139 

_i._ 2 .. 

1.20C57F-O4 

-l.60604f.--oi 

1 .6Ofc-7*fcE-01 

-89.9572 

GENERALIZED FORCES {NO 

thickness ffefcti 



“ODES 

PRES. DEFL. 

RFaL PACT 

I«a-3 PAPT 

A6S. VALUE 

phase angle 

2 1 

“2.4?s56E+00 

-2.45150E-01 

2.43TC2C*00 

-174.2287 

2 _ J 

-l,ft094pE*00 

-I.P1B26F-01 

1.61«*2S*00 

-173.5544 


mOOE NT. 1 


— 

- 

m_OOE_ No, ? 





GENERALIZED FORCES {WITH 

THICKNESS EFFECT) 



“ODES 

PRES*. DEFL*. 

_ PEAL PART 

Ifc-AG PAPT 

A8s, value 

PHASE ANGLE 

1 

1 

-1.05402E-03 

-2.41737E-01 

2.M739E-01 

-90.2498 

„ i 

2 

•8.9254IF-Q4 

-1.60300E-01 

1.60302E-0I 

-90.3190 


GENERALIZED FORCES (WITH 

THICKNESS EFFECT) 



MODES 

PRES. DEFL. 

RFAL PART 

IHfto PART 

A6$, VALUE 

phase ANGLE 

2 

1 

-2 ,42?83t ♦CO 

-2.305^1E-Ol 

2.43377E*00 

-174.5642 

2 

2 

-1.607455 *00 

-1.713^0e-01 

1.61656E*00 

-173.9147 


8i 



ASPECT RATIO 1.5 CELTA KING (T AU=0 •'! 0 ) 


ROOT CHnCQ LENGTH = 10. C* FT 

FRFP STREAM VELOCITY = moo, 00 '"T/SEC 


30 R'VES AL^VG Po n T CHORD 
REDUCED FREOUEnCY = .400 

FREQUENCY = f’.366E*00 CYOLE/SFC 

MODE NO, 1 
UQHE NO. 2 

9 

GENERALIZED FORCES (NO THICKNESS EFFECT) 


MODES 

PRES. DEFL. 

PEAL PART 

IMAG part 

48S, VALUE 

PHASE ANGLE 

1 

1 

2,ftC220f-02 

-P.S1330E-01 

9.51 736E-C 1 

-PR. 4332 

1 

7 

1,937125-92 

-6,284^2E-01 

6.2S740E-01 

-PR. 2345 


GENE 

PALI-7PQ FORCES (NO 

THICKNESS EFFrcT) 



MODES 

PRES. -PEPL. 

ARAL PART 

TMAG PART 

AoS. VALUE 

PHASE ANGLE 

2 

1 

-2.M5ce : ».if 

-9.9S9/2E-01 

2.6l373E*C0 

-157,5298 

2 

? 

-1.60030 r .*nn 

-7.44CT5E-01 

1 .T6<.63F*00 

-155.0635 


'•COE -NO. ! 

'■OCE NO. 2 


GENERAL 

I ZED FORCES fKITH 

THICKNESS FffECT) 



MO0FS 

PRES, DE r L, 

C EAL PART 

IMiG PAS'* 

ASS. VALUE 

phase angle 

1 1 

1 ,543Wt-0? 

-9.3970HP-AI 

9.39A93E-P1 

-39,0591 

1 ? 

1 ,PP?PPE-02 

-6.?O?95f-01 

O.2U416E-01 

-PR, 8706 

5f* E p.‘ i. 

TIED FCPCPS c ^ I T“ 

T<JD<NFSS E crc o T ) 



"MOOES 

PRES. DEP^. 

P^A_ PART 

Tw/,,1 =A = T 

A8s. VALUE 

PHASE ANGLE 

2 1 

.«f| 

-3. 61025* -0 ) 

2.58p?ie*C0 

-158.1326 

2 7 

-!.5c-696: *0 9 

-7.1?5~2>:-91 

l .T4.-'92F.00 

-155,6578 


82 



ORIGINAL' -BA&J5 IS 

Computer Program Listing OF PG0R QUALITY 



PROGRAM SBQXRT INPUT, OUTPUT ,T AP [ 5 = I NPU T, T AP£6»0UTPU T > 

SBOXR 

2 

c 

RAIN PROGRAM 

seoxa 

3 

c 

Ho = HAXI MUM NUPfcEK OF a'JXES ALONG THE CE NTE RL INE ' CHORD ALLOWED 

sboxr 

6 

c 

IN DIMENSION STATE tt NT 

sboxc 

5 

c 

ND=H AXI HUH HUPEfck OF MODES ALLOWED IN OlHENSION STATEMENT 

S50XR 

6 

c 

NB = HAX I HUH NUPock UF djXES ALLOWED FOR SPLINE FITTING 

SBCXR 

7 

c 

NK=N6+3 ALLCWEG IN ^InENSION STATEMENT 

SBOXR 

8 


COMMON xxl 2G),YY(3u> jXEDGHo) ,YED6H8)«XEDGI 32),YEOG(32) 

SBOXR 

9 


1 ,NS<2),nST(2I ,XTdGI<3»,YTDGK3J,XTDG(12I,YTDC<12) 

SBOXR 

10 


2 ,XLE <2,301 ,XTE<2,301,YnaX<2) , E DG ( 2 , 30 ) , AR I 30 , 30 , 2 1 

SSOXR 

11 


3 ,SFNX I iOOj ,s c iT HOOT, SFrml0 3J ,KSF N , T l i 03 , 10 A ) » 0 A ( 100 5 1 

SBCXR 

12 


A ,MLl2,3ol ,rlLT(2,30>,flLTU,2,301,MLC(2,2,30) 

SBOXR 

13 


5 » JHAX (2> » JNAK (2 ) ,CTE I 3 0 ) , AHA I 30,30,2 > 

SBOXR 

16 


& * A R £. A * CK » 0 > J ^ 1 i Iu!j t Ii' YLj^B>NtW»piH^NSM A"X y N E 0 G 1 1 NT 0 G I 

SBOXR 

15 


OlHENSION AT 2 , tiO,jOI , jl2,30»3G) 

SBCXR 

16 


1 ,SFOX < 100 ,3) .SF0YI100, 3 , 2 ) , SFOH ( 10’3 , 3 , 2 ) ,KSF C C 3 ) , G 1 3 ) 

SBOXR 

■ 17 


Z. tXcOGT(o) ,Y cOGTia»,/TOGTt3J,YTJGT(3J 

SBOXR 

18 


DATA Z/1H / 

SBOXR 

19 


IR=: 

SBOXR 

20 


IW=6 

SBOXR 

21 


NEOGI=0 

ssoxr 

22 


N TDG 1 -3 

SBCXR 

23 


HE = 30 

SBOXR 

26 


ND=3 

s SCX«. 

25 


NSH AX = 32 

SBOXR 

26 


N c - i 0 0 

SBCXR 

27 

c 

HATCh AtiOVc M.r'iL'O T; 1 l r I n£l JtNcNSIOl** IN STORAGE ARRAYS 

SECXR 

28 


HE=»S+1 

5BCX=. 

29 


nC*2«f u 

SBCXR 

30 


NN=Ad*3 

SBOXR 

31 


write ii-,??) 

SBCXR 

32 

r 

^tA f ) GAT A f-C* fur Al t Ja u *1*40 

SBOXR 

33 


100 Call 2-IPi. 16A-! 

5EQXR 

36 


Nt N=. 

SBCxa 

35 


T$5TI = _/M i JGI j **w*A( MCC J ) ♦»-A{ 11*1.4 ) +i,A f ICO ^ ) 

SBCXR 

36 


LK=j*<23)6CAi261 / jaU, J^S. 2?31obH 

SBOXR 

37 


D 1 =D A { 2 ? 1 

S b GxS 

38 


L =D1 

SBOXR 

39 


C= 1 ,0'/G I 

SEOXR 

60 


EN=n. b-0 

SBOXR 

61 


-SIT- <IS,,60) l „A4 X i , i = a, ? )- 

SBOXR 

62 


1 10 IF {l> cOt*cGC«I2G • 

SBCXR 

• 63 


120 IT lrie-L) ov'C , 130 , * 30 

SBOXR 

A A 


130 HklTE (IW,b-j) L*CA(2**)«L>,,0A(25I * A ( 2 3 ) 

SBCXR 

65 


IF (DAtZun 160,130.160 

.SBCXR 

66 


150 C AlL SHAPE 

SBCXR 

67 


IFINet.cO.2) Go TO 130 

SBOXR 

68 


A C = o * 0/ Ar t A 

sack’s 

69 


160 L 1 H=ML INEh.l J 

SBC-xi 

SC 


If ( L IP— MB 1 1?C, 170,630 

533*k 

51 


170 .LIN2 = 2»NB 

SoQXR 

52 


■NFLNS = 0At69> 

SBCXR 

53 


IF (NFLNS.EQ.03 HFLnS = L . 

SBOXR 

56 


LPCT = NINOIL ,NFLNS1 

seox'R 

55 


L1H1«2**8. 

SBOXR 

56 


CALL -.P0T2(LI«1,LIH2,LP0T,CK»0,A) 

S30XR 

57 


180 CONTINUE 

SBOXR 

58 


M=0 

SBOXR 

59 


K =DX<’<T»7 

'TB’JTO' - 

' 60 


GO TO 230 

SBOXR 

61 


83 





ORIGINAL' PA&S'ia 
OP POOR QUAIiTF2 


S1-AC»T<1.1> 


5 BOTH 


S2»ftC*Tl2,l> 

S3* S0RT<S1*»2.S29*2) 


SA»5?.295?8*ATAN2IS2.SU 
ASO WRITEtlHUOl H1.M2 »S 1.S2 «S3»S9 
30 FORMAT (1*1021 fc , lP3fc 19,5. OPiFlb. 4» 
»70 CONTINUE 


S88X*. "p33' 

SSQi&t a- - SHr 


ssmst-— -S32f 
SBO & ■_ -- -a33 


ssesa. 

osaSad- 


480 WRITE 1IH.55I 
■MOCaNTimjE 


IFTBAi-JU.GT.G.) GO TO 510 
DO 500 1-13,22 


560 BAtI1=2 


__ 2«i 

s goxir^ 

mgr** '■*** 




CHECK Uf ?RArrSFnRPI£0-TtS&& 



c 


IF SOat^-GO'-IO 

* Joj 


$H4?£ 


star re&tEaaigg&ssi 









550 

570 


3S0 


XT CG ( t J=XTOC! ( I ) 

YTDG I 1 J = T T Ov> I m " 

CONTI HUt 

CONTINUE 

IF INtX.Ea.2 ) GG T J 1J0 
NEW = 2 

IFI0A{2t>U 18C, 580,133 
IPAlNfOUlCOA)"* 


SBOXR 

‘SBOxR 

seoxR 

seoxR 

SBOXR 

SBOXR 

5B0XR 

SBOXR 



c 



CALL N*LC < Ji., r,Nn,N=, 

>.SFH,SFHX,SFKY,SFMH,iW, I PRINT) 

SBOXR 

181 


C-C To oG 


SBOXR ‘ 

182 


. 


SBOXR 

•183 


•t^RCx ; < i T b 


SBOXR * 

184 




SBOXR 

165 

5 4 n 

I p ' - .> c 


S 8 OXR 

186 


or i . t.i-j 


SBOXR 

187 

r ' C 

:»*=:? 


SBOXR ‘ 

IBS 

e- i 0 

^ ! i * (!>»•.£ ‘JJ-Ih*'- 


S&QXR 

189 

35 

f- I - v ~ t ( 1 o at tnr-i-j 

' 

SBOXR ■ 

190 


0 »T T * 7 J '•£ 


SBOXR 

191 

c ^ ' 

''-it: ! ! -, , *,0 


SBOXR 

192 

- - 

^ > * : 1 l 1' L.LMn’ L-I'-'-'l L 1 

w 1 I _'n r*U^ ::■< Gf- BGXt $ * cXC t£ D£0« ) 

SEOXR 

153 


, * .■* 


SBOXR 

194 

- - 

f i - " - 1 1 ; - *. i 


SBOXR 

195 

~ 1 

h i * w ► 1 1 *"i.i i’,^ i /-i. i 


SBOXR 

L96 

“ 

‘ M . < * l 7 | 1 » 1 . 5 * * o ' I j 

-L -*tvj * O'j T CriL^l>*r 

SBOXR 

19 7 


1. . . *.• 

. ,_-..,T- =,F3.2i 3N Pi/ 

seoxR 

198 



* • -Jt N Ct z if 3 i - 

SBOXR 

199 


1 --• ■. 1-- 

-*- ■ - Lt C i T Y =,Eb.2» 7h ET/SeC/ 

SBOXR 

200 


* .* ■ _ » ♦ , *• r • . - 

=..--Ii.3»IGn CYCLc/StC) 

SBOXR 

201 


. 


SEOXR 

202 
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DATA «ESt/5HAtPHA/*bTESt/lH /ijcTeSX/ltW, 5TES7/1H-»/. 


£Ss- 





< 4 .. -->gfr'* 5 S 


lift *21 • SttlN »AL*U I MS >»CRBU*1> » j^ii7l 




■ 

- s i ,, ‘ - 


- **&~ - r r 


DATRD 

OATRO 

DATRO 

a»Taa 

0 &B% 


~rsv n 


7 

3 


->9 ‘ 

SO 


Xf^TKMeOitfi Wtg zo 

IF tEttlN.ttf.SIEST) GO TO IPS 


C NEW W 1 HG- IF GO LORN 1 CONT A I#S A # 5 - VS S IGM 
C INITIALIZATION Of OAT A ARRAY 


DAI 
DATRO 


12 


DATRD 

0ATB8 


13 

14 


101 


00 101 - 1 - 23*1005 
0 AU(U* 0.0 


DATRD 

BATR 0 


102 


00 102 I»l »22 
DATAID-Z 


DATRO 

DATED 

'DATRO 

DATRO 

DATRO 

DATRO 

DATRO 

DATRO 

DATRO 

DATRO 


15 

16 


17 

10 


103 


DO 103 1 = 104 , 700,4 
DATA! 1 ) * 1.0 


19 

20 

XT" 

22 


TOT 


CONTINUE 
IF lALP.EC. 


A TEST), GO TO 9 


IF IALP.NE.DTEST) GO TO 8 

NUMERIC CARD 

DO" 3 1 = 1, fc 

DDR6U ( I )=0RbU(l) 


CONTINUE 

DECODE <60*9 90 ,CJRB U) ( 3< BU (I),I=li 

cc 6 i=r,v 

IK DR3U l_1 1) 4 »_fc ,4 

TEST' FOR BLANK FIELD 

IFISIGN < 1. C , DP.DU ( I ) ) IS, 5,* 


10 

II 


9 -< 1 
9 *3 
i -*1 


Si 


DATRO 

DATRO 

DATRO 

DATRO 


J A TA ( I NO l = CR3L(I ) 
I N D = I NO *■ 1 
OCTi) 11 
ALRr-A C A«v 


9 0 


Gi. JL. 


A - AT * f 


0AJ_AI lNJl=Dkc HI) 
'IND=INJ *1 
Ir Hi Nt . t TlS T ) 
hi. TURN I) CCLi-KN 1 
f. l TUi\ N 

t. *>D ?F 3 AT A CAkLS 
t. k I T < 1 6 ) 

tALL LaiT 
aTEH 

-a.. ;j-._ 

CONTI.Nk^ 


DATRD 

OATRD 

DATRO 

OATRD 

DATRD' 

DATRD 

DATRD' 

DATRO 


uG T_‘ 1 
0 JN T A 1 h 6 


A NiNUS SIGN 


^ * I N , Ai. P , l*' / i l j) J U < I ) ,1 = 1,7 ) 


< It. ,991) 


1 I c 

< in, 9 93) 

„i IT’ 

jl%. J . . 

. G* "A I t !•.' *0) 

t- T < 3dn BAG GAIA J‘, ThIS GAt-.Cj. JOB TE*MJ<ATEj) 
f « = **A U 12HCC ARC INaGc = a i ,A6, 16,7*10) 

f- Of "A_T t_/ / / ir , IPX , I 7BN J f L-R E OAT* CA=C/) 

• NC 


DATRD 

OATRD 

DATRO 

DATRO_ 

’ daTrd 

DATRO 

TXTrIT ' 
OATRO 
OATSIT' 
DATRO 
DATRD 
DATRO 
OATRD 

_DATRD 

DATRO 

OATRD 

DATRO" 

OATRO 


23 

24 

25 

26 
TT 

20 


29 

30 


31 

32 

IT 

34 

35 

36 

37 

38 


39 

40 

4 ~r 

42 

43 

44 
'T 5 

46 

4 T 

48 

"49 

50 


51 

52 


53 

_54 

55 

56 




SO 


OATRO 





ORIGINAL BAGD IS 
OF POOR QUALITY 



sr 1 

i f < H&r&mi. larrs 

I F 2 1 "GO SO 200 

TECGiS; ' " 


SHAPE 

SHAPE 

“SHAPE 

SHAPE 

“SHAPE 

SHAPE 


“15“ 

20 

"zr 

22 

“zy 

24 

25 

26 

TT 

28 


H5m=SXf2 SGT 

IF tttStl*} 550,850*100 
100 IE(MS«].9^E0&I*1I 105*105, 850 
105 NSP=HS«1»*1 


HSUT-BSP 
IF IOAt24I I 855*855*110 


nooirn 

XEDG(ll-OA<2*.l*271/OA(241 


SHAPE 

SHAPE 

Shape 

shape 


115 YESSTrySA 1 2* IT73JV0A<24» 

XEDG«1)«0,0 ' 

Nsrm=T5inR5i 


"29 

30 

3T" 

32 

SJ- 

34 

“35~ 

36 


SHAPE 

SHAPE 

“SHAPE 

SHAPE 

SHAPE 

SHAPE 

"SHAPE 

SHAPE 


IF (NST 111 J 871,125,130 


425 XT0SI11*>1.1S 

NSTI1I=1 

gu~nrT4i 

130 IFIH5TUJ-H10GI*!) 135,135,871 


135 

DO 140 I=»l*tiS?P 

mmsssn^PsiToimi 

140 XTDGtSt«0Af241«44i/0A{24) 


“37“ 

38 

“39“ 

40 

41 

42 


SHAPE 

SHAPE 

SHAPE 

SHAPE 

SHAPE - 

SHAPE 


141 HSTPi-HSmiil 

HSTIII-pHSTP 

HTSSIKSYP fSlECCTISn 

YTDGIHSTPI'TEBGIHSPJ 


“43“ 

44 

“sy 

46 

"47 

48 

49 

50 


YiBG<u-a.fl r 

if JMSTIH„£0.NTCGI1 00 TO 150 

ji-NSTUm 

DO 142'I=JL,HTB6I 

xTD6iU=xTG£II-lJ 

142 mrosn^vmii-i] 

“rW“irr*35r5Tfff53HSWT-XEDG<W5?l7YlfUG<HSIP]-YE0GTHSPTr.Trr;ERKRr 


SHAPE 

SHAPE 

SHAPE 

SHAPE 

SHAPE 

SHAPE 

“SHAPE - 

SHAPE 

SHAPE 

SHAPE 

SHAPE 

SHAPE 


63 TO 871 


~W" t6& “I* IVftl 061 

xE06iin>XEsein 


“5T“ 

52 

53 

54 


“yMUg i^vtusi j* 
IrCI.GToHIDEIJ 00 TO 160 

^bm'cir^msrn . ■ 

VTOSUli-vtOOIll 


SHAPE 

SHAPE 

“SHAPE 

SHAPE 


55 

56 
“5T 

58 


160 CGHTIKOE 

GO TO 550 

c ! 

C FOR KCOIFIED MHG 

C COMPUTE HfltH KW8E1 OISTRIBUTIOH OH PHYSICAl U1H6 

200 DH*0. . 

DO 210 K«l»2 . * 

DO 210 

DO 210 1-1*88 - 

'210 AHA! J,I»KW.C 

00 215 I - 1 »L 


SHAPE 59 

SHAPE 65“ 

SHAPE 61 

SHAPE 62 

SHAPE 63 

SHAPE 64 

5HAPE 65 

SHAPE 66 

SHAPE 67 

SHAPE 6i“ 

SHAPE 69 


87 













ill.l-tH.EQ-O) GO TO 215 




m*£h 


>» 



«-ML#?¥i*iU : ~ 7':T~~: 

CALL SURF 2 4 XX 41 » tT» , J»X » It AMAH . i s 1» *0M,DM <SF8X» SF(4Y tSFHHiKSFr til 

— 

SNARE * 

73 

215 

CONTINUE 

Stipi 

rx^ 


iMjejTE.EQ.CJ GO *0 235 

SHAPE 

7» 

C 

PRINT OUT MACH DISTRIBUTION »* ' 

SHAPE 

7* 


WRITE 1 IB *60 > 

SHAPE 

77 


00 230 I«i.L 

SHAPE 

— nr 


IF(HLC41,1,X i.EO.O) GO TO 230 

SHAVE 

73 


JL*NLC«1,2,I )-NLC4 ltltl 1*1 

SHAPE 

«« 


IFIJL.EO.O) GO TO 230 

SHAPE 

*1 


WR ITE (IW t 70) I 

SHAPE 

82 


JLP= JL/6 

SHAPE 

83 


IFt JL~fcPJLP.NE.OI JLP«JLP*1 

shape 

85 


4U«HLC41,1,I 1-1 

SHAPE 

85 


JL*HlClit2tI> 

SHAPE 

86 


DO 225 J-I.JLP 

SHAPE 

87 


K-Kl* J 

SHAPE 

88 

225 

WRITE t IK ,801 4 !J1,AMA< Jltl.DIt J1«K,JL,JLP» 

SHAPE 

89 

230 

CONTINUE 

SHAPE 

■TO 

C 


SHAPE 

91 

235 

CONTINUE 

SHAPE 

97 


K*1 

SHAPE 

93 


HSP«l.*DA427>/2. ♦ ERRR 

SHAPE 

95 


02»2.*B 

SHAPE 

95 

C 

DEFINE EDGES GF MODIFIED HIM* 

SHAPE 

96 


DO 250 I=Z*NSP 

SHAPE 

97 


XEOGI I)<FLGATtI-lJPDZ 

SHAPE 

98 


IF 4XEDG4 D.GT.XED&KKI) JC«K*1 

SHAPE 

99 


TEOGI I J*YEDGI CK— 1 1*4 TEOGI IK1— YE06HK— 11 1P4XEDG4 11— XEOGI IK— 11 1/ 

SHAPE 

100 


i 4XE0GHKJ— XE0G1CK— 111 

SHAPE 

101 

250 

CONTINUE 

SHAPE 

102“ 


IF ( IF 1X4 DA4 2? > >-*2*CNSP-l ) ) 810,255.250 

SHAPE 

103 

. 250 

HSP-NSP+1 

SHAPE 

105 


NSl«NS(l) 

SHAPE 

105 


XED&INSP1— XEOGI C NS 1 1 

SHAPE 

106 


TEOGINSP 1- YEQG It NS 1 > 

SHAPE 

107 

255 CfflSTIHUE 

SHAPE 

Ida 

c 


SHAPE 

109 


KW 

SHAPE 

110 


A2*«$P 

SHAPE 

111 


NST-NSltl 

SHAPE 

112 


TO 300 I«2.*S1 

SHAPE 

113 


SO 260 K«IL1»k2 ■ 

SHAPE' 

nr 


K3<«- % 

SHAPE 

115 


.* t . 1 * /v i { I jl ' J ) ,lT j } 0<, TO 2S0 

SHAPE 

116 


if(> . j { 5 } . •, v *Ac.^^l(i)} 3 ~ f ^ j7- 

SHAPE 

nrr 


L-'^n --C!. 

SHAPE 

ne 


.. * ^ . 1 1 

SHAPE 


_ 7 ) 

f G 5* ^ tA * 

SHAPE 

120 



Shi FiT " 

121 


S*** — J 

SHAPE 

122 


Mcn-^r^xLccir. » - . . 

SHAPE 

T2T 

* r C' 

n * *i ) =t •_ -c u » 

SHAPE 

125 


■ " ‘* * 

shape 

I2F 


>E9G t a3 J = *fc iiC 1 1 1 1 

SHAPE 

126 


YfcCG(Aj)=YEuC*[ (1 ) 

SHAPE 127 


K2=>:2*1 

SHAPE 

128 

T 

Ki*< in 

SHAPE 

129 

300 

COHT i hUE 

SHAPE 

130 


IMkHITE) 310,320,310 ‘ “ ' "“SHAPE ITT 

c 

PRINT OUT Llli INC tOGtS TO 3E TRANSFORMED 

SHAPE 

132 

313 

Kl*l 

"SHAPE 

133 


WRITE! IWtfcSJ 

SHAPE 

135 

NStKfc 

SHAPE 

135 


WRITE 1 1 H , 501 K 1 ,NEN <NS ( hEM) , 1EDG 

WRITE! n,.5oT _ TK«'Xt OCI K i . yffiEliI .a^i 1'tTCVT — 

SHAPE 

136 

nr 


320 CONTINUE 


i Tf2T> 


SHAPE 

SHAPE - 


138 

W 



ORIGINAL'' BAGl IS 
OF POOR QUALITY 


c 


SHAPE 

1*0 

c 


SHAPE 


c 

TRANSFORM Y-CCCRDINATE 

SHAPE 

142 


TTO 350 I»1*NSP * " 

SHAfE 

143 


YMP-O. 

SHAPE 

144 



SHAPE 145 


3»1_ 

SHAPE 

146 



SHAPE 

147 

330 

CAlfc SURF2C XECGt I } * Y* 1* l*l*£tt»Dtt*DMDY*5FHX*SFMY* SFMH«KSFH*3I 

SHAPE 

148 




1W 


TeK»rS-uMpv 


SHAPE 

150 

y«b=-ymd - 

l^tVko.eE.9.1 60 TO 335 


“SHAPE 

SHAPE 

I5T 

152 

i if- 


SHAPE 

SHAPE 

153" 

154 



SHAPE 

SHAPE 

1ST 

156 



SHAPE 

SHAPE 

157" 

3.58 

~C '^6^ 
- ^M«£St$^»,«,CI.6.4»DH» 60 Td si 5 


SHAPE " 
SHAPE 

' 155 
160 

' fiaRQR * 1 T I J* U4J-1 * 1E*1 


SHAPE 

SHAPE 

1ST" 

162 



SHAPE 

SHAPE 

1ST 

164 

'- «r- 


SHAPE 

SHAPE 

165" 

166 

Y-^EBGII* 


Shape 

SHAPE 

168 



SHAPE 

SHAPE 

160 

170 



SHAPE 

SHAPE 

1737 

172 


IF(KRITE) 360*370* 360 


SHAPE 


173 


C 

360 

PRINT OUT LEACIMG EDGES OF TRANSFORMED WING 
Kl-2 

SHAPE 

SHAPE 

174 

175 

176 

177 



WRITE 118*401 K1*NEU*NS{NEU1«1EDG 
WRITE (18,501 (K.XEBGdO* YEOGtKl »K=1»NSP) 


SHAPE 

SHAPE 


370 

CONTINUE 


SHAPE 

178 

C 




SHAPE 

174 

c 


READJUST THE PGIHT DISTRIBUTION TO DEFINE 

NS 1 2 1 

SHAPE 

180 



Kl-0 


SHAPE 

181 



iCHECK-0 


SHAPE 

182 



YKP=0. 


Shape 

183 



DO 450 1=2* NS F 


SHAPE 

184 



YNP = AHAX 1 1 YMP * YEDG ( I )> 


SHAPE 

185 



1FIYE0GI I 1-YfcCGC 1-1 i ) * 00 * -.2 J* **0 


SHAPE 

186 

c 


ACCEPT Small H/«>E otHiNO TRANSFORMED "IMG 

LEADING EDGE AS PART 

SHAPE 

187 

c 


OF rflNG* IF NOT MORE THAN 1/5 OF A ?0X 


SHAPE 

188 


<•00 

T ( 1 * 1 1 =X£CC II) 


SHAPE 

189 



T 12*1 )=YEDg< i 1 


SHAPE 

190 



T (3 *1 ) = AtCG t i~ll 


SHAPE 

191 



T(4,l »=Y£DC< I - 1 > 


SHAPE 

192 



TI5,l»=YnP 


SHAPE 

193 



I ER=5 


SHAPE 

194 



IF lYMP-YEDOl I 1 .GT.O.*,*3rO 62 TC iC5 


SHAPE 

195 



I£RP0R=40C 


SHAPE 

196 



«RITclI»i,£.;j ILRiv.<jtt{T(J«i)*J=l*l£?) 

. . 

SHAPE 

197 



YECGl 11 = rrf 


SHAPE 

198 


4?6 

IF ( I •C'-**ICrc.v_* , '-*ll JU ^ 


SHAPE 

199 



ICHECr- 1 


SHAPE 

200 



L C TC *t C 


Shape 

201 


*. 30 

ICh£C*.= t 


Shape 

202 



KX = *v 1 ■*■ i 


SHAPE 

203 


4 AO 

K= i-K 1 


SHA»E 

204 



XtDGt X 


SHAPE 

205 



y E r g i * ) *r* . ; t ; i 


SHAPE 

206 


tSC 

COfTI-’-', 


Shape 

207 



IF is- 1 1 C (• )) .7 .:»U ) o : T„ »02 


SHAPE 

208 



a<) 3w^i*ic 1 i 


Shape 

209 


39 






c 

n60 

4 7.0 
480 

NS I 21 *K — 1 
60 TO 480 

NSI2MK ’ 

N5P*N<»<2_) : • ... 

IF(KiUTt) 450 , 495 a 490 

* , - ♦ 

SHAPE 

SHAPE 

"Shape 

SHAPE 

■ shape 

SHAPE 

210 

211 

■ 212 
213 
”214 
215 

£ 




SHAPE" 

" 216 


490 

Kl=*3 


SHAPE 

217 



WRITE ( iw"»40) Kl,'NEn ,NS<NEWT,IEDG 


"SHAPE 

2TTT 

* 


NR ITt ( I H 1 5 0 ) (KtXEOG(K) ,YE0GtK) »K*1,NSP) 


SHAPE 

219 


495 

CONTINUE 


5HWIT 

220 

c 




SHAPE 

221 



TRAILING EDGE TRANSF ORNA T 1UN 


SHAFT" 

7Z7r 

c 



* 

: SHAPE 

223 



NST tNEWMNST (1) 


SHAPE 

224” 



NS TP =NS T ( NEW ) 


SHAPE 

225 



CALL SURF2<'xTCG,YTuG,l,NSTP,0,TIl,l>,DM,DM*SFHX,SFMY»SFHH,tCSFM,lI 

SHAPE 

226 



DC 510 I = 1 » N STP 


SHAPE 

227 


510 

YTDGI I ) =YTDG III) *T.( I . 1) 


SHAPE 

. 228 



IF(XRITE) 520,530,520 

" 

SHAPE 

229 

c 



* 

SHAPE 

236 




SHAPE 

&£■££! m 



MRIT€fIW*50J <e„3?0G*R»VYTDG«O ,K«1,HSTP» 

' 

SHAPE 

232 


TTiTT TH3Bfc? : y < 5 ^ ■■*, i- : * - . - ■ - kj 



233 


550 

Y1-YEOGI1) 


SHAPE 

234 



YMAXfN£UI-YEB6«K$?> 






IFfYia.T.0.01 GO TO 060 


SHAPE 

236 





■fi ; hrjm 

237" 



AREA- 0.0 


SHAPE 

238 



XXI1MDH 

_ 

B1FJ24 

239 



VYUJ-OH 


SHAPE 

246 



DO 560 1*2 ,ilB 



h M 1 | m 



VVUMYYII-1M0 


SHAPE 

242 


5&0 

XX(1I*XX(I-11«0 




c 


, 


SHAPE 

244 

c 




SHAPE - 

iKiVl 

c 


COMPUTE AND IDENTIFY SOX DISTRIBUTION OH THE NIHG 

PLAMFORM 

SHAPE 

246 


565 

CALL PLMFMt XEB6, YED6,XTDG , YTDG,XX, YY»XLE ,XTE »AR, AR£A,D 

M.TI7.TJ4I 



V ,IM,LtHB,HL,HLC*HLH*MEM,HS,HST, IRITE) ' 

SHAPE 

248 

T 




"SHAPE 

2A9” 



IFIHLOtEUvLl.GT.HB) GO TO 825 


SHAPE 

250 

c 




SHAPE - 

251 



JMAXI H£Mi*0 


SHAPE 

252 





SHAFT" 

253 



□0 570 I»1,L 


SHAPE 

254 



JNAMHEWMMAXfl <JHAK<NEM> ( HLC(H£H,2,1 ) ) 


SHAPE 



570 

JHAXI HE M > =*MAXC 1 JNAXCNEW » , ML < HEM , I )> 


SHAPE 

256 





SHAPE 

Z5"7 

c 




SHAPE 

258 

L 


»»«- 1111 ■ I U— ■ 

"SHAFT" 


c 




SHAPE 

260 



JHA-JHAXIHEU 1 “ 


SHAPE” 

261 



00 620 J-1,JNA 


SHAPE 

262 



eugYHen,j»«o.c 


"SHAFT" 

263 



00 620 K-1,2 


-SHAPE 

• 264 


625 

MLTtK£H,K, JMO ' “ ~ 

- 1 

5MAPE - 

26T 



11-1 


SHAPE ‘ 

' 266 




* 

SHAPE - 

267” 



Xl-0 


SHAPE 

268 



00 655 J»l, JMX “ ' - 


"SHAPE" 

269 


630 

IF(HLC(M£M,2,I1I.GE.J) GO TO 632 


SHAPE 

270 





■MILM4 

JVfl 



GO TO 630 - 


, .SHAPE 

272 


6J2 

HZMI 1 


SHAPE 

rrr 

c 


FIND ORDER OF TRAILING EDGE BOX AT 3-TH CHORD H1SE 

ROM 

SHAPE 

274 


634 

IF< I2.GT.L) GC TO 638 


SHAPE 

275” 



IFIKLCIKEM* 2*121. 6T.0» Kl-1 


SHAPE 

276 



IF (Kl.EQ.O) GC TO 635 


'SHAPE 

277 



IF(MLCiK£H,2,I2>.EO.0.AN0.NZl.NE.0} GO TO 646 


SHAPE 

278 




-SHffPT" 

279” 



IF(KLtI(KEH,I2>..nE.O> GO TO 636 


SHAPE 

280 


or 










ORIGINAL BAGE Ife 
OP POOR. QUALITY 




12-12*1 " 
K2-1 




~ 5H5T>£ 
SHAPE 

2*2 


636 

GO -TO 634 

IFCHLS-J) 644*646*638 




SHAPE 

284 


6 36 

IFIK2.LT. Of GO TO 648 



* 

SHaWE" 

2*> 



IF4HLT(TKEW,I21 J 642*642**46 




SHAPE 

2«6 


642 

12-I2-*K2 




SHAPE 

_ 2*7 



IFU2.LE.L) GO TO 634 
H22-L 




SHAPE 

SHAPE 

311“ 

2*9 



12-4. 

GO TO 652' 



- 

SHAPE 

SHAPE 

299 

2*1 


644 

K2-HLWtK£H,I2>/m.S 
GO TO 642 




SHAPE 

SHAPE 

292 

293 


646 

NZ2-12-1 

GO TO 652 




SNAP# 

SHAPE 

294 

295 


648 

652 

N22-I 2 
CONTINUE 




SHAPE 

SHAPE 

296 

297 



NLT4NEW* 1* J1-NZ1 
MLT(NEW,2*J)-MZ2 




SHAPE 

SHAPE 

2^8" 

299 

c 

655 

CONTINUE 




SHAPE 

SHAPE 

"399“ 

301 

c 


FIND ORDER OF RING BOXES AT EACH I-TH SPANWJSE 
00 670 I-l.L 

COLUMN 


SHAPE 

SHAPE 

302 

303 


662 

IFiHLMNEH,!) ) 664*662*662 
NLC(NEH*1*I )«PLW4N£H, I ) *1 




SHAPE 

SHAPE 

304 

.305 



IF(HLC(NEW»2 » I I .EQ.O) HLC4 NE H * 1 »I 1-0 
GO TO 670 




SHAPE" 

SHAPE 

306 

307 


664 

MLC(NEN,2,I>«IABS(MLW(NEW,I) >-l 
HLCIN£W,1*I J-l 




SHAPE 

SHAPE 

308 

309 

c 

670 

CONTINUE 




SHAPE 

SHAPE 

310 

311 . 



IFURITE.EO.O) GO TO 678 
WRITE i IW.65) 




SHAPE 

SHAPE 

312 

313 



IF (NEH.EQ.l ) hRITtUH,771 





314 



IF4NEW.E3.2) HR1T6<IW,78T 




SHAPE 

315 



WRITE4IW,82) 

DO 672 J=1,JHX 




SHAPE 

SHAPE 

316 

317 


672 

WRITE <IH, 85) J»HLT(N£W»1»J),HLT (NEW,2,J> 
WRITE (18*83) 




SHAPE 

SHAPE 

318 

319 



DO 674 J-1»L 

WRITE! 18*85 1 J ,HLC (NEMil* J)*HLC (NEH*2* J1 




SHAPE 

320 

c 

674 




SHAPE 

321 

“322 

323 


674 

b78 

WRITE 418*85) J *NtC( NEW, 1 * J) * HLC (NEW, 2 * J) ,HLH(NEW, 
CONTINUE 

J) * NL (HEW, 

J) 

SHAPE 

SHAPE 

T 






SHAPE 

324 



IF(XEDG(NSP).GT.l.) GO TO 865 




SHAPE 

325 



1F(N£W.EQ.2> GC TO 680 
NSN-NSP-1 




SHAPE 

SHAPE 

SHAPE 

SHAPE 

326 

327 
328' 
329 


680 

DT=TEDG(NSP)— YEDGCNSrt) 
CCMTINUE 






I F (ABS(DY) » LE »ERRR ) IEDG-1 
IFIDY.LT.O.) GO TO 827 




Shape - 

SHAPE 

330 

331 

c 


COMPUTE VALUES FOR LEADING EDGE CORRECTION 




SHAPE 

332 



YPAX2— YHAX {NEW )*YHAX ( NE W ) 




SHAPE 

333 



DO 720 J-l* JHX 

IF ( IEDG ) 710,715,710 




SHAPE 

SHAPE 

334 

335 

c 

710 

D IN FOLLOWING EXPRESSION IS ARBITRARY 

FDG4KEW*J)-S«3t< (YHAXZ-YYl J)*YY t 41 >/(0* (2.4YHAX (NEWI-O) II 


SHAPE 

SHAPE 

336 

337 



IF t£DG( NEW* J) .GT. 1. 0 ) GO TO 715 
GO TO 720 




SHAPE 

SHAPE 

338 

339 


715 EDC(KEV*J>-1.0 
720 CONTINUE ■- 




SHAPE 

SHAPE 

340 

341 

T 






SHIPS - 

342 



IF(4tEN.EQ»U 90-70 750 - . 




SHAPE 

343 

~c 


COMPUTE Y-COORDIHATE AND MACH NUMBER OF POINTS 

ON 

PHYSICAL 

WING 

SHAPE 


c 

c 


AT TRAILING EDGE OF TRANSFORMED WING 
AT 60X CENTERS OF TRANSFORMED KING 




SHAPE 

SHAPE 

34* 

346 



DO 745 J-l, JHX 
M-NLT1NEH*1»J1 




SHAPE 

SHAPE 

34^. 

348 



K2-NLT (NEW *2 * J )— K 1*2 
I —XI 



v. * 

SHAPE 

SHAPE 

349 

-no. 


91 






00 740 K«1*K2 
K3-0 


SHAPE 

SHAPE 

351 

3M 

DEL-0. 

X-XX<1> 


SHAPE 

SHAPE 

3** 

354 

IF(K.eo.K2» X-XTEINEtf.Jl 
Y-YYIJ) 


SHAPE 

SHAPE 

355“ 

356 

730 Y-Y*DEL 


SHAPE 

35T 

K3-K3+1 


SHAPE 

358 

C FIND HACH MUHfiER *1 X.YOH PHYSICAL WING 


SHAPE 

— 3 =PT 

CALL S0RF21 X*Y»1»1*1»EH»0N»EH¥» Sf NX. SFMY»SFMH*KSFN»3> 


SHAPE 

360 

S1«YY(J)-Y*EH 


SHAPE 

361 

IFIABSISl J-Lfc.ERRRI GO TO 735 


SHAPE 

362 

1FIK3.GT.160J Cfi TO 836 


SHAPE 

JbT 

DEL- Sl/(£N*Y4£HY+l.E-20) 


SHAPE 

364 

GO TO 730 


SHAPE 

365 

735 IFIK.NE.K2J AKA( J*I.2>-EH 


SHAPE 

366 

i F 1 K. EQ.K2 ) CTEIJ )-EH 


~SHA>E 

367 

740 1-1*1 


SHAPE 

368 

745 CONTINUE 


SHAPE — 

364 

750 CONTINUE 


SHAPE 

370 

c 


SHAPE 

77X 

790 CONTINUE 


SHAPE 

372 

C 


SHAPE 

* '373 

RETURN 


SHAPE 

374 

■y 

“ ‘ ’ 

'SHAPE - 

' 375 

800 IERROK=0OO 


SHAPE 

37 fc 

T ( 1 * 1 ) =X£0G ( K ) 


SHAPE 

377 

T(2,l >=YEDGIK) 


SHAPE 

378 

IEk = .? - -- - 

' 

SHAPE 

379 

GO TO «40 


SHAPE 

380 

805 ITRRijk"= 8j5 - - - . .. . . . 

' ' 

SHAPE 

'381 

GO TO 840 


SHAPE 

382 

610 IERRl)R = olO 


SHAPE 

383 

T (1 , 1 )=0A<27> 


SHAPE 

384 

T ( 2 , : ) =N 5 P 


SHAPE 

385 

IER=2 


SHAPE 

386 

GO TO '340' - - -- 


SHAPE 

387 

c* 1 5 1 f RROR = o 15 


SHAP.t 

388 

GO 13 3-.0 


SHAPE 

389 

820 1E-R2R-320 


SHAPE 

390 

T < 1,11=aE0GTn5F> 


SHAPE 

391 

T (7 ,1 ) = YEGG IMP ) 


SHAPE 

392 

t 1 3, i 


SHAPE 

393 

; f l = , 


SHAPE 

394 

GC T ' 04y 


SHAPE 

395 

c 1 1 - - 1 " = e 2 *» 


Sr* AP t 

39 fa 

Til.: )=.N L r 


5 nA P t 

397 

i (2 


SHAPE 

39 8 

TI3.:i=.-. u (S t .,rl 


SHAPE 

399 

T(*ill*lL(>.ir ,K-1 1 


SHAPE 

400 

i t 1 -4 


SHAPE 

4 Q 1 

oO TO 8 4C 


SHAPE 

402 

827 I F 4RCR = a2 7 


SHAPE 

403 

T 1 * i i ) - YE 05 IMP ) 


SHAPE 

404 

T < 2 * i )*Yti,G(iV>r) 


SHAPE 

405 

TIj..J=-jP 


SHAPE 

406 

t ( 4 i 


SHAPE 

407 

I IF --4 


SHAPE 

408 

ol T. v - 4 J 


SHAPE 

409 

8 3 J I E. K RGis ^ 73G 


SHAPE 

410 

T ( 1 , 1 3 -- 1 


SHAPE 

411 

T I ? ,1 )- J 


SHAPE 

412 

T ( 3 . 1 3 = a 


SHAPE 

413 

T I *♦ . A J * Y 


SHAPE 

414 



Shape 

415 

T <e,_ )»i«v 


SHAPE 

416 

T ( 7 1 1 3 = J. t. 



•* i 7 

It 1 *-? 


SHAPE 

418 

5*9 W- IT: <!W,3m) l T ( t . 1 ( , 1 = l , 1 L4 3 


SHAPE 

419 

STCf" 


SHAPE 

420 


92 



ORIGINAL' RAGS If 
OF POOR 


S6C IFR--.9 



SHAPE 

421 

CC Tj ,0$v 

t 


SHAPE 

422 

-s 5 I P c =2"** 

1 


SHAPE 

423 

GC. TO £9U 

• 


SHAPE 

424 

860 IPR=30 



SHAPE 

425 

vO to . 39 c 



SHAPE 

426 

665 K- MI*»GIK,N5I1I J 


" '* 

SHAPE 

427 

IFR=2*K+Z9 



SHAPE 

42* 

GO TO 890 



SHAPE 

' 429 

871 !Pf=55 



SHAPE 

430 

89C NRlTt 



SHAPE 

431 

C Uw,10)!PR 



SHAPE 

432 

1C F0RMAT(1H0,10X,1 7HSHAPt — BAD DATA, 15) 



SHAPE 

433 

20 FORMAT l 1H0, 1GX ,38HbA0 NUMBER IN SHAPE NEAR 

STATEMENT NO, ,15, 

SHAPE 

434 

. i_ . ‘./IN - ,iSX,lF8E.I4.b) 



SHAPE 

435 

AO FaRNAT</lrfO,5X,3HNO.,I2,43H REDISTRIBUTION 

OF MING 

LEADING EDGES, 

SHAPE 

436 

1 K S ( , 1 1 , 4h ) * ,I2,7H» I EGG, 3H = ,11) 

1 


SHAPE 

437 

45 F0RMAT!/1H0,5X,25HWING TRAILING EOGES , NSTI,I1,4H) 

” ,12, 

SHAPE 

430 

J lOh, IEDG = ,11) 



SHAPE 

439 

50 FCRMAT(4t6*,I 3.1P2EU.4) ) 



SHAPE 

440 

60 FCRMsTI IH1,5X, A7HL0CAL MACH NUMBER DISTRIBUTION 'TJN 

PHYSICAL WING/) 

SHAPE — 

- • 441- 

65 FORMAT! 1H1) 

f 


SHAPE 

442 

70 FORMAT! 1H0,5XW-Z,19H-TH SPANWI5E COLUMN) 

/ 


SHAPE 

443 

75 F0RMAT«lHC,5X,I2,17rt-TH CHORDNISE ROW) 



SHAPE 

444 

77 FORMAT (AOX, 13FPHYS I CAL WING///) 



SHAPE 

'445 

78 FORMAT! A0X,l6ETRAMSFDRME0 KING///) 



SHAPE 

446 

60 FGRMATI1H , 5X ,6 ! 2X , 13, 1PE13. 51 ) 



SHAPE 

447 

82 FORMAT! 15X , 68FCR0ER OF F IRST! LEADING) AND 

LAST! TRAILING) WING BOX 

SHAPE 

440 

1IN CHORDVISE RCM//20X, 1H J, 3X , 12HHLT 1 NEH ,1’, 

J) ,3X,12HMLT(NEW,2,J)/) 

SHAPE 

449 


1ANVISE CQLUHS//20X,1HI «3X,12HKLC!NEW,1, lit 3X,12HHLC(H£W,2,I) 
i »VX,llH!,MLh!NEM,I),5X» 10HHL ( NE H » I ) J / )f 
8 5 FORMAT! lfeX,I5,5X,I5<3( IPX, 15)) 


SHAPE 

SHAPE 

SHAPE 

SHAPE 


451 

452 

453 


86 F OR MA T! 2DX » 58F.CAUS SIAN INTEGRATION POINTS IN CHQRDWISE ROW - GX4ME SHAPE 454 
lw,K,JJ//J / SHAPE , 455 


87 FQRM*T!///2CX, 52HGAUSS I AN INTEGRATION POINTS IN SEH1-SPAN - GYIHEN 

1,K>//) •’ - 


SHAPE 

SHAPE 


455 

437 


88 FORMA TI///20X, tOHSPANH I S fc COORDINATE' idYP) AND MACH NUMBER <CMP) A SHAPE 
5T"G'xI2',K,J)//1 ‘ SHAPE 

80 FORMAT (1H ,5X, 3(2X, 13, 1PZE16. 6) I t SHAPE 

"STOP " * . SHAPE" 

END SHAPE 


93 





94 



ORIGINAL' mm » 
OF POOR QUALITY 


G>. T„ 


PLNFH 

71 

2fiP I f t Y-'VL + r *■ S- . C E . D ) GO TU 290 


PLNFH 

72 

AM j,l ,st „ J »AR ( J, I,NtW) ♦ AS NG» « Y L-Y+U > /O 


PLNFH 

73 

A 1 = W 


PLNFH 

74 

Y1 = YL 


PLNFH 

75 

IFUtriC.hL.il CO TO 310 . . 


PLNFH ;■ 

76 

1FIJO 284,285,285 


PLNFH 

77 

■ 234 TFrVT2.i.-t;Y-DJ'''Ci;m _ 321Cr ‘ 


PL3FH. 

78 

GC TO Zttb 


PLNFH 

79 

235 IF(Yl2. iL.Y) tC TO J2J 


PLHFH 

80 

256 CONTINUE 


PLNFH 

81 

IFIaEL.GT.aL) CC TC 400 


PLHFH 

82 

If J l=Z 


PLHFH 

83 

GC TC 460 


PLHFH 

S4 

290 X1*X2 - 


PLHFH- 

85 

"Vl=Y2 


PLNFH 

86 

295 1 1 D2 =0 


PLNFH 

87 

tfnCHC.Ke.n GO TO 310 


PLNFH 

88 

Y *Y€ 2 


PLNFH 

89 

CG TO 515 


PLHFH 

90 

310 YP=Y21,TNGMXP~aEII 


PLNFH 

91 

315 IFfJKr JT7T3TETm ~ " ~ 


PLHFH 

92 

316. JF-fYfi .&E.-Y) GC 70 i20 


PLNEH . 

93. 

GC TO 318 “ “ ' 


PLHFH ‘ 

9* 

317 IFCYR-LE^Y-EL GO TQ 320 


PLHFH 

95 ■ 

• ?18 CONTINUE 


PLHFH 

96 

1FCICH0.E0.11 GO TO 325 


PLHFH . 

97 

IF1H-1 


PLHFH 

98 

X2 = XR 


PLHFH 

99 

YZ = YR 


TLNFN 

100 

GO TO 330 ' 


PLHFH 

101 

320 jFIn=1 


PLHFH 

102 

IF(JK) 322*321,321 


PLHFH 

103 

321 Y2*Y 


"PLHFH' 

104“ 

CO TO 323 


PLHFH 

105 

. 322 Y2=Y-D 


PDtEH 

rosr 

323 X2*XE1*(Y 2— Yt 1 1/TNG 


PLHFH 

107 

IF CABSIX2-XR J .GT.ERSR) 60“W3“5fl 


PLHFH 

T08 

I F I N= 1 


PLNFH 

109 

JFIH»0 


PLHFJf 

116 

GO TO 330 

- 

PLHFH 

111 

3Z5 IEOZ-1 

“ 

TLKFH - 

112 

X2=*E2 


■ SLHFN 

ll3 

Y2-YE2 


PLHFH 

114 

3 30 AR4 4, 1 , HE k 1 »A8 4 J, i,NEW) ♦0.5*ASHG4f2.*XR-X2-XlI*tY2 

-Yi»/D2 

PLNFH 

115 

C 


PLNFH 

m 

c . - - - 


- stjmt. 


I F <A 3S t TNG ) .LT.ERRRJ TNG-l.E-20 

, 

PLHFH 

118 

1331 CONTINUE 

* -f 

PLNFH 

11* 

IftJXl 332,331,331 


PLHFH 

120 

331 IFf J.LT.JLTJ GO TO 334 


PLNFH 

121 

SU*XElt1Y»»JUl-Y61UT«fr 

- 

PLNFH 

122 

IF4XLT.GT.XE2. AND. A5NG.GT. 0.01 GO TO 334 

'* A 

PLNFH 

125 

l F ( YFDG4 1 1 .GT.0.0.AN0.YY4 JLT1.LE.YL. AMD.ASHG.GT.O. 

0) XLT-0.0 

PLHFH 

124 

1FCYY (TLT).GE .YTOGt NSTP 1 J XLT*XTDGCNSTP } 


PLHFH 

12* 

1 F 4 ASHG.GT . O.C 1 XLE 4 HEN, 3LT *«XLT — - - — 

- . , - - -A 

• - PLNFH ■ - 

324k- 

. F4ASKG.CT. O.OJ XTE4NEW, JLTJ-XLT 

- 

PLNFH 

127 

GE TO 333 


PLHFH 

128 

-332 IFIJ.GI.JLT) GO TO 334 


PLNFH . . 

- 129 

XLT=XE 1*1 YY 1 JLT l-YE 11 /TNG 


PLNFH 

130 

TF4YYCJLT1.GE.YTOG4N5TP)) XLT-XTDGINSTP 1 


PLNFH 

131 

XTEINEW,JLT)=XLT 


PLNFH 

132 

• 333. JLT*JLT*4K,_. . 

. 

PLNEH 

.133 

GC-TG 1331 — ' 


PLHFH 

134'- 

334 CONTINUE 


PLHFH 

135 

C 


PLNFH 

136 

IFC IET2.E0.1) GO TO 400 


PLNFH 

or 

SAK=SAK*ARt J,I,NEH1 


PI NCV 

138 

' IFCJRr 344*342,342 - - - 


PLNFH"' * 

X39 

342 CCHTINUE- 


PLNFH 

■ 140 




IFI3FIN.EQ.1) GO TO 250 
€0 TO 350 

PLMFM 

PLMFM 

141 

142 

344 

IF(-tFlN.EQ.l) GO TO 555 

PLKf N 

— lis - 

350 

IF! AS MG I 355,360*360 

PLMFM 

144 

355 

I F ( If IN. E G« 1 ) GO TO 579 

PLMFM 

i4S 


GO TO 530 

PLNFH 

146 

360 

IFUFlN.EQ.il CO TO 565 

PLNFH 

" ivr 


GO TO 500 

PLMFM 

148 

C 


PLMFM 

159 

C 

LEADING EDGE AND TRAILING EDGE COMPUTATION 

PLMFM 

• 150 

c 


PLMFM 

1.51, 

500 

IF(ASNG) 550,505,505 

PLNFN 

152 

505 

IF ( JK ) 530,510,510 

PLMFM 

153 

510 

Kl«Kl+ JK 

PLNFM 

154 


IF C K1 .GE.NSP) " GO TU 520 

PLNFM 

155 


XE1=XE0G(K1) 

PLNFM 

156 


YE1=YEDG(K1) 

PLNFM 

157 


K*K 1+ JR 

PLNFM 

159 


”xe"2*xedg i k i : 

' 1>IRFH 

" 1T9" 


YE2-YEDGCK1 

PLNFM 

160 


"GO "TO" '550 

PLNFM 

161 

520 

1FCJK.LT. 0) GC TO 250 

PLNFM 

162 


JK — 1 

PLNFM 

163 


JLT=J 

PLNFM 

164 


MLC (NEW, 2, I 1=C - - - - 

PLNFM 

165' 


K2=K2~JK 

PLNFM 

166 

530 

K2=k2+JK 

PLNFM 

167 


XE 1 = XTDG(K2 ) 

PLNFM 

168 


YE 1=YTDG ( K2 ) 

PLNFH 

169 


K *K2 ♦ JK 

PLNFM 

170 


XE2=XTDGCK) 

PLNFM 

171 


YE 2= YTCG { K > 

PLNFM 

172 


GO TO 550 

PLMFM ' 

173 

W 

"TT2>lC5Tr 

PLNFH ~ 

174 


XEl’XTOGC K2 1 

PLNFM 

175 


YE1*=YTDG(K2 J “ _ ‘ 

PLMFM 

176 


XE2-XTDGCK2U1 

PLNFM 

177 


VEZ-YTDGCK2+1J 

PLNFH 

178 

; 550 

C=XE2-XE1 

PLNFH 

179 


F*YE2— YEl - 

~ PLNFN 

180 


TNG-F/IGU.E— 20) 

PLNFM 

181 


1FCG.LT.0.1 GO TO 650 " “ 

PCNFU 

182” 


IFCNEH.EQ.2) GO TO 555 

PLMFM 

183 


AREA=AREA-,ASNG*G*( YE2+YE11 

PLMFM 

184' 

c 


PLNFM 

185 

555 

Ton tinge - 

PLMFM 

186 


ICHQ=0 

PLNFM 

187 


IF CXR.GT.XE2J ICHQ=1 

PlmFN 

' 188 


1FUED7.EQ.1) GO TO 290 

PLNFM 

189 


IFC IE02.EQ.2) GO TO 295 

RLMFN 

190 


GO TO 250 

PLNFM 

191 

c 


PLNFH 

192 

C 

number of boxes In spanhise column 

PLMFM 

193 

c 

HAKE TYPE AND ORDER OF HAKE BOX ... 

TlNFW “ 

1'94 

c 


PLMFM 

195 

565 

IFCASNG1 579,567,56? 

~POFH ' ' 

1T6' 

567 

IFCJKJ 576,569,469 

PLMFM 

15/ 

4C V 

J J-0 

"PtNFF 

T9S~ 


IFCAR fJ,I ,NEH).LT.ERR2) JJ-1 

PLMFM 

199 


ML(NEH,l)-j-jl "■ ' 

FLWFH ‘ 

200 

C 


PLMFM 

201 


j'J='MtCTfETfn'S " ' ““ - •• 

PLNFH — 

"202 


CO 471 K=1,JJ 

PLMFM 

203 

Is fs -i\ 

"PLNFH 

" 204 


IF ( A* (K •> I jh£H ) *LT .ArtAK) GO TC 472 

PLMFM 

205 

571 

CGNTTNOF - - - 

PC.NFHT 

_ 206 


MLC CHEW ,2*1 ) = KK 

PLNFH 

•207 


GC "TO 473 7 ‘ ' — • 

PLNFH ' 

208 

5 72 

ITi 

MLC 1 NEH,2 , 1 ) *K K— 1 

7*7*0 TTUi'fC — 

PLNFN 

209 


W3 C'CNITNuF — — fLHKH 210 


96 



original' BAG® is 

OF POOR QUALITY 


c 

ADJUSTMENT OF FIRST BOX AT LEADING EDGE ALONG KK-TH CMORDMISE ROM 

PLNFM 

211 


“ KK-HLCOtEB7?m 

IFIXXII».GE.XLECNEH,KK1) GO TO 474 

"PLKFH ” 
PLNFM 

212 

213 


NLC 1N£H»2»T>”MLCINEM»2, 1 1— I 
GO TO 473 

PLMFH 

PLNFM 

214 

213;. 

c 

474 CONTINUE ■ ' 

PLNFM 

PLNFM 

214': 

217 : 


GO TC 560 
476 MLS— 1 

PLHFM 

PLNFM 

2W 

.219 


«L INCH. J 

K-HLOtEW,Il 

plhFI 

PLNFM' 

>20 

an 



453 AR( JJ,I ,N£H)»C.O 
484 CONTINUE 

I F (ML h(NEW,I)» Nt • U ) 60 TO -490 
I F ( AS NG > 467, 486, 486 
mIG'TIP HAKE 
486 IF <ARIJJ,I,NEH),GE.AHAKI CO TO 490 


HLVINEH,! )=jJ 

60 TO 4 88 

Vi I Hj ROOT HAKE- 


437 I F ( A RJ J J , UN EH )»LT »AHAK) GO TO 490 
NLH(NEW,n«JJ-l 

488 ML V t NE H , I J = NL W <NE rt iIIAHLS 


PIWH 

PLNFM 

PLNFM 

PLNFM 


PLNFM 

PLNFM 


490 CONTINUE 

go t: 500 


TRAIL IMG EOCE aox akea adjustment for king tip hake 


PLNFM 

PLNFM 

PLNFM 

PLNFM 


495 CONTINUE 


PLNFM 

PLNFM 


237 

23®. 

239 

240 


241 

242 

243 

244 
24 5 

246 

247 

248 


j= j-i 

Y=Y-J 
JF IN=0 

SAM=sA 1 — A4(J,l,NcHJ 

GC TJ 290 ‘ 

“Flnfm 

PLNFM 

PLNFM 

PLNFM 

"PLNFM 

PLNFM 

249 

250 

251 

252 

253 

254 

END OF A SPANHISE CuLUrtN 

PLMFM 

255 


PLNFM 

256 

i CONTINUE 

PLNFM 

257 

IF IASNC.GT.O. J GO TO 516 

PLNFM 

258 

ADJUST WING TIP BOX FOR H ING ROOT HAKE 

PLNFM 

259 

J J=Nl (NEH, I ) 

PLNFM 

260 


810 IFCAR(JJ,1,FEH). LT . A3AK ) GO TO 512 

HLC t NE K,2,1) =JJ 

GO TO 516' ” 


PLNFM 

PLNFM 

PLNFM 


261 

262 

263 



512 

JJ=*J J-l 


PLNFM 

264 



IFIJJ) 514,514,510 


PLNFM 

265 


514 

MLCINEH,2,I)«C 


PLNFM 

266 



IFINLWINEH.I) -EQ.O) GO TO 516 


PLNFM 

267 



HP ITE ( IW, 1 55 } I»I,, ML HIKE H,I),1»MLC(NE H-, 2*11 


PLNFM 

268 



STOP 


PLNFM 

269 


516 

CONTINUE 


PLNFM 

270 



SUM-SUH+SAM 


PLNFM 

271 

C 




PLNFM 

272 


5 30 

CONTINUE 


PLMFM 

“273 

C 




PLNFM 

274 



IF(ASNG.LT.O.C) GO TO 700 


PLNFM 

275 



IF{XTOGI1).GE.XTOG(NSTP)1 GO TO 700 


PLMFM 

276 

c 




PLNFM 

277 

c 


INITIALIZATION FOR WING ROOT HAKE 


PLNFN 

278 

T 




PLNFM 

WikM 



K2-0 


PLNFM 

280 


97 






ETflilFII4!4 


2 5*i 2 1HC0MPUTAT 1 ON CONTINUES//) 

1^0 "formak ihIi^x.s^hl&aoiSg’ano trailing edge points at EACTTcRtiRDWl 
r *E R0H/13X»lHJ,tX» 7HLEA0ING»8X-»8HTRA1LINC/) 

145 FORH AT I 10X • I 3 * IP 2E 15 • 5 ) 

150 fORMATI 1H0,4X,29HPLNF« — NEGATIVE VALUE NEAR S H, I4,2Xt A2 .1PE11.4) 
155 FORMAT! IH0,4X»25H0R0FR OF KING TIP BOX In ,l2.43H-TH TpSSWTSf '£01 
1MN IS NOT P ROPERLY D EFI NED/ / 

2 15X,8HHLW!NEW,«I2,3H1 » , 1 3/15X, 8HHLL(H£W, , I2«3H> «»t3 


PLNFM 

PtNFH 

PLKFM 

>UfH 

PLNFM 

PtNfH 
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ORIGINAL' BME 1» 
OF POOR QUALITY 



SUBROUTINE PQT2IH2*H0,H0*CK,D,AI < 


POT2 

2- 

C” 



PDTZ 

3 

c 

THE VELOCITY FIELD OF A ONIFORH DOUBLET DISTRIBUTION 


PDT2 

$ 

c 

OV£Sd' * "eO» Ti 'CO«PuTEO AT' ALL POINTS AT MHICR It Uli-L Bfc T" 


POT2 

T 

c 

NEEDED AND STCRED IN THE ARRAY A IN COMMON 


P0T2 

6 

c 



POT2 

~t 

c 

HO (NO CONTROL THE NUMBER OF VALUES CONFUTED 


POT2 

B 

c 



WITZ 

9 

c 

H2 IS THE RANGE OF THE SECOND SUBSCRIPT IN THE ARRAY*" 


POT2 - 

to 

T 

-DTHENSTOMetf *< 2,M2,N2), BUT TREATED" hEWAT AH ARRAY 


P0T2 

11 

c 

RITH TWO SUBSCRIPTS 


P0T2 

12 

“C~ 



TBT2T 

13 


DIMENSION A12,1) 


P0T2 

14 


K»H0 \ 


PDT2 

19 


H«NO . 


POT2 

16 


DK«€*:*D 


POT2 

17 


D*2»0X**2 


P0T2 

18 


TCCTFT 


FBT2 

19 


DXB-DA2/8.0 


P0T2 

20 


DX4»'2.'ff*DK8 


P0T2 

21 


DK12-DR2/12.0 


P0T2 

22 


CN-if.5 


P0T2 

■ 27 


DH-DK*0„5 


P0T2 

24 


OM»0.6*Dft 


POT2 



00-2. 0*DK 


P0T2 

26 


DDH»DD 


FUTZ 

27 


D1»0.25*DK2 


P0T2 

28 


B 5=0X27 Z4T0 


FDTZ 

29 


DO 3 1=1. H 


POT2 

30 


'BT=TJT0 T 


PUT 2 

.... 31 - 


84*2. O/DH 


P0T2 

32 


B2=83/ B4-UH 


POTZ* 

7T 


B3=-0.5*65 


POT2 

34 


D3=UH*B4.B5' ' ' - - 


P0T2 

35 


D4*DK8*B4 


P0T2 

36 


. - - - - .... 


P0T2 

37 


CN=1.0 


P0.T2 

38 


"K^I 


FU.TZ 

"39 


C3-C.0 


_ P0T2 

40 

- 

C4=~07Tr~ . - 


POT 2 

41 


C7=0.O 


P0T2 

42 

. 

C d^TJ.Tr 


P0T2 

43 


DO 2 J*1,N 


PCT2 

44 

— 

Al=Ort/CN 

1 

P0T2 

45 


C1=CN* CJSIA1) 


- POT2 

46 


C2=-tH< STMA-I) 


P0T2 

47 


C5=Ch*CIN(Al»C61 


POT2 

48 


C6=-CK*C6 


PGT2 

49 


C9=C1-,C3 


PCT-2 

50 

■ 

XTa^m^5 ' “ 


P DT2 

51 


C 1 X =C 5— C 7 


POT 2 

52 


C 1 2 ~C"c-C £ 


P0T2 

53 


A ( 1 .a ) = b3*C9-K44Cl6-35PL3-i>l$Cll-t2*C Id 


P CT2 

54 


At? )-D*i^C9->£ j^Cl ! J- 1 i ~ t 1*C 1 c 


P0T2 

55 

23 

C3=C1 


P0T2 

bb 

" * 

C**C'2 


PCT2 

57 


L 7 = C5 


POT2 

58 


Cb=Co 


POT2 

59 


81=dl-01 


P0T2 

60 


tS3 = S3-D3 


PtfT2 

.61 


84=34-04 


P0T2 

62 


04=04.004 


POT2 

63 


CN=CN+2.D 


POT2 

64 

2 

K = K-*3i2 


P0T2 

65 


( 


PQT2 

66 


i,r-Cn*jOrt 


POT2 

67 

3 

DL.“=.dLfl-.DU 


POT2 

68 


GO K 11=1.2 


POT2 

69 


*1 = 1 


POT2 

70 


99 



00 5 J = I,N 

P0T2 

71 

oo 4 1 = 1 , m 

P0T2 

TT 

K= K 1 ♦ h— 1 

P0T2 

73 

4 A Ul_,K) =A ( IL.Kl— A( 1) 

P0T2 

74 

A 1 IL,K1J=2.C*A (iL.kl) 

P0T2 

75 

5 K 1 = K 1 ,(12 

P0T2 

” 76 

CM*0.0 

P0T2 

77 

0(1=0. 0 

PUY2 

TT 

DDn=QA 

P0T2 

79 

JO 12 1=1, (1 

P0T2 

80 

C c 

P0T2 

*1 

ct=o.o 

PQT2 

82 

C9=0.0 

P0T2 

83 

CiC=O.D 

pcrz 

84 

Pi=0.0 

P0.T2 

85 

P 2 = C . C 

P0T2 

86 

CN=1. 0 

P0T2 

87 

66=0. 5*0K 12 

P0T2 

88 

K=I 

P0T2 

89 

00 10 J = l,fi 


iS^KX*a 

Al=Crl/CN 

P0T2 

91 

A2=DM/CN 

‘ POTT' ' 

92 

IP < A 1— 0 • 2 ) 7,7,6 

P0T2 

93 

7 el=4.G~Al**2/3«0 

P0T2 

94 

b2=-0K/ <6.0*CM 

P0T2 

95 

GC To 9 

“P0T2 

96 

8 b3= S IN ( Al ) 7A1 

P0T2 

97 

5 1=2. 0$ 8 3 

por2 

' 98 

82 = 1 83— CCSCA1 ))/A2-0H/CN$B3 

P0T2 

99 

9 6 j= C0S1A21/CP 

P0T2 

100 

64= SlN(A2)/0 

P0T2 

101 

C3 = Sl$d3* 82*04 


102 

C4=62*63— 81*34 

PDT2 

103 

35=QH*CN 

P0T2 

104 

01=88*04—2 . G* C 3 

P0T2 

105 

C2=-2.G*C4-t5*C3 

" ' P0T2 ‘ 

106 

C 5 = C 1-C7 

P0T2 

107 

Cb=C2-Co 

P0T2 

108 

P3=P2— Bb*CN 

P0T2 

109 

P4=P3»2.0*DK12*(CN— 1.0) 

P0T2 

110 

A( 1,K J = A( 1,K ) ,Cb-P 1*C6*P3*C3-P4*C9 

P0T2 

111 

A{2,K) = AC2,KJ j »Cfe + Pl*C6 + P 3$C4— P4*C10 

"P0T2 

112 

P 1=P 1+DH 

P0T2 

H3 

P2=P2+CN*DK4 

>0T2 

114 

CN=CN«2 .0 

P0T2 

115 

C7=C1 

P0T2 

116 


‘rs-cz * P0T2 nr 


C9=C3 

C10=C4 

B6=Sb+DK12 

10 'k»k+«2" 

CM=CM+DK 

0H=DM*DOrt" 

12 DOM=DOM*DO 

D3=CK/ ( 2 .C* 3. 1 41 592651 
fU=M2-H 
K = 1 

.‘.1 = 0.0 
ottt4 j = 1 1 n 

01=03$ SIMAi) 
C2=-D3 * _ 'CCSTaT) 

DO 1 3 1 =1. h _ 

DF£ =A(1,K)$C'1-»V(2;k>*C2 
A(2,K )=AI2,K)*C1-A( 1,K)*C2 
’ "aTT7KT^T5FE 

13 K=K+1 _ 

R = K*wl " 

14 A 1 = A1 +Orl 
" ift TO? fc 

£ND 


PDT2 11* 

" P0T2 " * 114 
P0T2 120 

TOR 121 

P0T2 122 

P0T2 123 

P0T2 124 

" ‘ 'PUT? 125 

P0T2 126 

POTT 127 

P0T2 128 

FTJT2 ITT 

P0T2 130 

POT 2 13 r 

P0T2 132 

TUT2 113" 

P0T2 134 

FDT2 135" 

P0T2 136 

PUT? ” T3T 
P0T2 138 

pot 2 no 

P0T2 140 
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ORIGINAL' 

OF POOF QUALU i 


V “i-UTlNfc 08 E 0 C SF Ox .SFDY.SFDH, KSFO.SFMX.SFHY.SFHH.KSFM. OA.T OREO 

v , A a ,“Y Y ,T* ; L‘, KL'Ci'>iThC)','H'BVREHiHH»KSrS7IP PINT J 'DRED ' 

OlMFNSION SFOX (N6.ll ,SFOY < NB ,HD , 1 J iSFDHCKN.HO, II , KSFDU. ) OREO 

1 -,5FMX (1) ,Sf«Y(l> fmhi i I »OA 111 ,T«NH,1).HLC<2'2,1> ORHT 

2 ,XX(I),YY(i» OREO 

t BKEIT 

£RRR * 1 . E — 06 OREO 

ia'lT€'-U>RlHt/lC “ " BRED 



JK1TE=IPRINT/100 OREO XO 

Tfinew.eqVi 1 TRl’E t -Tpr I NT- lb>XRTt£ ■ ‘ BRED TIT 

IFCNEW.E0.2) GO TO 400 OREO 12 

KSFD«H)=0Af9<J» ■ ' “OREO I3~ 

Np =K SF 3 { M ) OREO 14 

""If (np> ?3g, I7C,ioo DREO T5" 

C OREO 16 


C A SPlINE-SURFACE FOR THE DEFLECT I OfTTS FITTED - TIT VALUES WEED 


C OF DEFLECTION AT GIVEN POINTS. OREO 18 

C ' ‘ ' OREC ITT 


ICO IF (NB-NP) 73C.120.120 

~~IZH TCRTTSuE 

KP=100 

3L 140' I °= 1 ,N F 
SFDX(JP,M)=0A(KP+1)/0A(24) 

5 FC vTT P ,N,1 >«DA<KP+2l/0A(24j 
SFDHt I P . M , 2 ) * C A{ KP +3 1 

SFCH{ Ip-,M,i) = $F0HUP,rt,2I 

140 XP=KP+4 


OREO 20 

BRED Z1T 

DRED 22 

“ ~ TWEO ~ZJ 

DRED 24 

DRED 25* 

DRED 26 

OREO 27 

OREO 28 


C 

c 


c 


i 5 

spTTne-surf ace nTTJTTS 



DRED 

OREO 

29 

30 

COiTlNOt 

IF CIRITl ) 15C. lbO. 150 



DRED 

OREO 

31 

32 

T5C 

** i t i mao* 



. — 

DRED 
OREO 
‘ TIRED - 

"33“ 
.34 
35 

IfcC 

CCNTlNUt 



DRED 

36 


CaLl SURFUNh.KSFOlrtl , T , SF DX ( 1 , MTTSFD Y ( 1, H, 3. ) .STOW* l.M 

nms 

rrrti 

"BRED 

rr 



GO TO 200 

DRED 

38 

c 

PRESENT LTTm PITCH AND PLUNGE DP FORM Z-AO+AITX 

BRED" 

3 9 - 


170 SFDH(1»H,1)=DA(52J 

DRED 

40 


SFDF*(2»N,1 J = 0 A 1 53 1*DAI 24) 

“DEED - 

~4T 


SFGhl3,N,ll«0.0 

DREO 

42 


Go TO 145 

OREO 

43 


200 ICHECK=0 

OREO 

4.4 


Tit TURN 


— 45 

c 


DRED- 

46 

c 

FOR NOD 1 F I fc D RING 

OREO 

4T~ 

c 

i-CALCULATE DEFLECTION ON THE ORIGINAL WING ' 

OREO. 

48 

c 

2— TRANSFORM X, Y, 0EF TO THE TRANSFDRHEO PLANE 

DRED" 

49 

c 

3-SURFACE FIT DATA 

OREO 

50 

~c 


DRED 

5T" 


40C CONTINUE 

OREO- 

52 


IF(KSF3<M).EQ.O) GO TO 500 

OREO 

53 


xP-KSFCtMI 

DRED 

54 

~c 

CALCULATE LOCAL MACH NUMBER OH PHYSICAL MING t INPUT PQIMTSI 

DRED 

55" 


CALL SURF2 ( Sf D XM.HI .SFOY ( 1 ,M , 1 ) ., 1 ,KP 1 0.TI 1*1) .QUN.QUH.SFHX* SF HY 

OREO 

56 



Kill 

rr 

c 

STQRF TRANSFORNEOURt*' 

OREO 

58 


ic j - ; , r » 

OREO 

54 


45G lF„V( I ,M,2) = SFLY< I ,11,1 mi I.1J 

OREO 

60 

c 

FIT -EFLcCTION data IN TRANSFORMED SPACE 

OREO 

61 


Ct l- 5J2H(Nr,KP,T,SF3Xll.H),SFDY(l,H,2) ,SFDH< 1 ,H» 21 .IRITE » 

DRED 

62 


- i Tc- N 

OREO 

63 


D?J ICHeC^=!Cri£uK *1 

DREO 

64 


I F t ICnLC*. * u T • 1 ) oj TO !?60 

DREO 

65 

c 

eel' T •-•iir.Ts cn transformed wing for deflection surface-fitting 

DRED 

66 


-> = ' 

DRED 

67 


jl.CT<L,f.A,NEW,MLC’,KSFS,IW 

DREO 

~ 68~ 


i ,XX, TY.SF .'XU.M) .SFOYti.M.ZI.JRITE) 

DRED 

69 


IF U- 1 Tl.GT.O) rt.ITt IIW.55I “ ' ‘ 

DRED 

“ ‘ 


: 101 



n o ’ I o 1 n r> o 


COMPUTt Y-CLLPCIHAT,. Jr POINTS OH PHYSICAL WING CORRESPONDING 
"TJ T 11 OS b SlLECTlO On T* ANSFU-vHED MING “ 

00 5Sii K«UK5PS . 

0£L=O.C 

X«ST3XIK,H) ■ _ . 

Y«SFDY( I ,H,2I 

520 Y = Y ♦*>{ L _ 

“FiNO'rtACH NGhcERAT {\”,YI OH PHYSICAL' WING 
CALL SURF2(X,Y,l,l,l,E“,DUrt,£MY,SFMX,SFNY,SFHH,KSFN,3) 

Sl“SFOY(K,n,«.)-Y*_n 
IFlAdSCSlJ.Lt.tRRK) ^0 TO 5J0 
DEL*S1/ (LN + Y*t PY*l.r-.’ J) 

GO TO 520 

STORE Y OF PhYSICAL jPaCl 
530 SFUY(K,H,l)=SFDY(k,rl,2)/<EH-rl.t-20> 

550 CONTINUE 

GO TL 5B0 _ ... 

560 DC 570 r. = l ,KSf S 

SF DX< K* lij *S F OXJ! K , MR ) _ 

SF C Y 1<,H»1)=SF~CY(K t Nr* 1 1 ’ " " 

570 SF0Y(K,N«21=SFLY(<,hA,2) _ _ _ 

5 BO CONTINUE 

COMPUTE DtftbCTlON AT SELECTED POINTS ON TRANSFORMED KING 
BY CALCULATING DcFlECTICN OH Tut CORRESPONDING POINTS ON THE 

PH YSICAL WING __ 

” 'CALL SlFRF:<'5FUX<1,M),SFDy7i,M,1),1,'kSF5»0,T(1,1),DUH»0UM 

i _ ,$F0X(l,M^ f SF0Y(l»M»lJ»5FDH( 1 ) t KSFD (N) *1) 

00 610 i = lT,KSFS 
_6_I0 SFDHI 1 tMf NE H I »T( I , 1_) 

IFURITE) 620,630,620 
620 WRITl(IK,6C> f 

KlUTUIXaO] 

‘630 CON f I HUE* 

FIT OEFLECTION - 0?~Tr ANSFORHED' RING BASED ON SELECTED POINTS 
CALL SURF1(NM,K$FS ,T,SFDX<1,MJ , SFDYl 1 ,M,2 ) , SFDHd ,N,2» 


t , IR I TE } 

RETURN _ _ __ _ 

730 IPR«98 

GO TO 750 _ 

750 'HRITE 

0 1 1H ,20) IP R ■ 

STOP ... - — ■ 

1°_F QRHA T< 1_H 0 » JL 0 X » 7 3 HC OH P U_T E 0__DE FlECTION = A0*A 1»X»A2 »Y* SUN. OF HI 

It 1 * (I)A*2)*( AL CClKljl®*'))) 

20 FCRPATI i n C , 1 C X ,1 6H J j — .AT/-, l?) 

30 »Ur‘*AT(lriU,BY,26MPHTSI CAu. Pi-.v-t — ".6: NC . , I 3 ) 

AC LC-mA T t LnC , 6> , 2 yhJ R A - , SF J- * •_ ,■ — “jOc. NC.,13) 

56 FORMAT ( ini ) 
i nE 


DRED ’ 
DRED 
DRED 
’ DRED 
DREO 
ORED 
DRED 
ORED 
DRED 
DRED 
DRED 
OREO 
DRED 
DRED 
DRED 
DRED 
DRED 
DRED 
DRED 
DRED' 
DRED 
“DRED 
DRED' 
DRED 
ORED 
DRED 
MED 
DRED “ 
DRED 
~0rETT" 
DRED 
OREO 
ORED 
"DRED 
ORED 


II *_ DRED 

D RE D’ 

DREO 

DREO 

DRED 

ORED 

DRED 


102 


original 
OF POOR 


SUBROUTINE SLC T <L .NB.HE W«HLC< ISLCTiIW.XX.YV.XI* YI « 3 RITE I 
(NUMBER OF SpAtWISE ROWS! 



c 

_c 

c 

c 

c 

c 


L 

NB 


INPUT 

INPUT 


(HA XI HUM MUH&ER Of POINTS ALLOWED *I<- OIWEMSlOUf- ■ 
(ORDER OF FIRST ANO EAST Sit® BtW 1 -. 1 M CCttMNl 
(SELECT PROW 1 -PHYSiCAL* 2 -TmtSF 0 RME 0 # 


NLC 

NEW 


INPUT 

INPUT 


4 RITE- INPUT (INTERMEDIATE RESULTS PRINT-OUT CONTROL! 
XX«Y.Y— INPUT (COORDINATES OF SOX 6 R 10 I 



DIMENSION MLC(2*2»l»*XXtl*»YY(i»*Xmi,Vmi 
HTOT-O 

SLCT 

SLCT 

♦ 

» 


DO 20 1*1*1. 

stci 

KU I 


IF(HLC(NEW.l*I).EQ.OI 00 TO 20 

SLCT 

« 


NTOT«NTOT+HLC(NEH*2»II-HLC(NEW. 1*1 1*1 

SLCT 

13 

20 

CONTINUE 

SLCT 

14 


NP-0 

ste? 

vi 


NPT*0 

SLCT 

16 


1-0 

SLCt 

17 


K-l 

SLCT 

1« 


DJ-0.0 

SLCt 

19 


DN“FLOAT(NTCT ) /FLOAT (NB) 

SLCT 

20 

1 

IFIDN.LT. 1.0) ON “1.0 

SECT 

21 

! 

IF( JRITE.EO.O) GO TO 25 

SLCT 

22 


WRITE t IH.200) 

SLCT 

23 

L 

IF(NEH.EQ.l) HRI TE ( IH,210) 

SLCT 

24 

r 

I F (NEK .EQ. 2 ) HR ITE ( IW, 220 ) 

SLCT 

25 

L 

HRI TE t IH.23G) 

SLOT 

26 

1 25 

NP“NP+1 

SLCT 

2 1 

! 30 

IF(K.LE.NPT) GO TO 50 

SLCT 

28 

F 

1 = 1 + 1 

SLCT 

29 


NPTP=NPT 

SLCT 

30 


IF (HLC(NEH«i«I I.EO.O) GO TO 30 

SLCT 

31 


NPT=NPT+MLC(NEH,2.I )-HLC(NEW,i, I)*l 

SLCT 

32 


GC TO 30 

SLCT 

33 

50 

xi(NPi=xxm 

SLCT 

34 


J=K— NPTP 

SLCT 

35 


Y I < NP ) =YY ( J ) 

SLCT 

36 


IF(JPITfc.CT.O) WRI TE C 1W , 100) I* J*K*XI ( NP ) tYI(NP) 

SLCT 

37 


IFCNP.GT.N&) CC TO 60 

SLCT 

38 


DJ“DJ*DN 

SLCT 

39 


NJ“DJ+0.6 

SLCT 

40 


CJ-0 J-FLOAT(NJ > 

SLCT 

41 


n = K+fU 

SLCT 

42 


IF (K.GT.NTOT ) GO TO 80 

SLCT 

43 


GC TO 25 

SLCT 

44 

ro 

WRITE ( IW,25C) 

SLCT 

45 


CC 70 1 = 1, NP 

SLCT 

46 

7C 

WRITE ( I M., 11C) IiXlllliYMU 

SLCT 

47 


S TCP 

SLCT 

48 

eo 

CONTIhWfc 

SLCT 

49 


I SLOT =NP 

SLCT 

50 

* ' “ 

Kt TU, S 

SLCT 

51 

100 

f-CR*AT(5X, Ji3,lP4Ll4. ?) 

SLCT 

52 

no 

FCRNaT(1h ,10,14, IP+tl*-.?) 

SLCT 

53 

2 CO 

Fomn ini,*x ,i9MSt<_s.cn priMi >.n i 

SLCT 

54 

2 IU 

FC-RNi.1 1 lH+,2 jX , L2<( PmTSICAL ) 

SLCT 

4 6 

220 

FCRh^T(lH”»^23X, I2«7 i^Afijr ) 

SLCT 

56 

230 

FGRrtATllH+t3?X 1 3Ir , »-Irto Fi'K S r Ll>c — SU-sFACc FITTING// 

SLCT 

57 


S9X,Ul I ,4X , ihJ , 4X ,1m< * *\ .5**1 (R) ,UX,&HYI.(\)/ ) 

SLC.T 

58 

25C 

FLR‘*AT<ltiO,5X,59HNUna£R OF POINTS 5ELECTED IN SLCT EXCEEDS 

ALLOW SLCT 

59 


iloLt LTKIT77T4T, 1HI , 3T,3HXmT,'9XV5HY H I »/) 

SLCT 

60 


ENO 

SLCT 

61 


103 . 




SUBROUTINE WVAL(XX,YY,XTDG,YTDC,XTE,SFDX,SFDY«SFDH.KSFO»S,0,CK 

OVAL 

2 



1 ,IH,L,NL,M1,NEH,NST,MB,H0.NB »NN,KSFS,IPRIHT) 

HVAL 

3 

r. 


CALCULATES OOHNHASh VELOCITY DISTRIBUTION I REAL AND IMAGINARY I 

HVAL 

4 



DIMENSION XXIl),YY(lJ,XTDGCl),YT0G*l>,XTEI2,l),SI2,MB,l),NLt2,ll 

HVAL 

5 



1 , NST (T)»SFOX(NB»1)»SFOY INB, HD.'l 1 ,SFDH(NM, MD,1),KSFD(1 1 

HVAL 

6 



IRITE=IPRINT/1C 

HVAL 

T 



IFINEH.EQ.il 1RITE=IPRINT-10*IRITE 

HVAL 

6 



JLP = KSFDC Ml ) 

.HVAL 

9 



1 F t KSFDIM1 I.EC.O.AND.NfcW. EQ. 2 1 JLP-KSFS 

HVAL 

10 



JL*ML (NEH.L ) 

HVAL 

11 



DO. 5 1 = 1. L 

HVAL 

12 



DO 5 j=l.JL 

HVAL 

13 



DO 5 K= 1 ,2 

HVAL 

15 


5 

S(K,J,I)*0.0 

HVAL 

15 



DO 80 1*1, L 

HVAL 

16 



JL = NL INCH, I ) 

HVAL 

17 



IF(JL) 50.8C.5C 

HVAL 

18 


40 

CONTINUE 

HVAL 

19 



00 70 J = 1 , JL 

HVAL 

20 



IF(XX (I J.GT. (XTEINEH, J)+D) > GO TO 70 

HVAL 

21 



CALL SURF2< XX ( I) ,YY ( J I , 1 , 1 , 1 , V ALU, VALUD . DM, SFDXI1 .Ml > 

HVAL 

22 



i » SF~C Y < 1 »M1, NE H ) .SFDHC liNl,N£H»,JLP,21 

HVAL 

23 



S Cl, J.I >=VALUC 

HVAL 

24 



S(2,J,1I = CX* VALU 

HVAL 

25 


70 

CONTINUE . 

HVAL 

26 


"80 

CONTINUE 


27 



IFCIKITE) 1C0.200.100 

HVAL 

28* 


100 

'HP Iff ' IIH.iOi 

HVAL 

29 



IF CNEH.EQ. 1 ) HRITtC l H » 12 1 Ml » 

HVAL 

30 



IFCNEM.EQ.Z) HRITECIH.15) Ml 

HVAL 

31 


10 

FORMAT ( lrtl ,10X .57NUPHASH (REAL, IMAGINARY, ABSOLUTE, PHASE ANGLE ) I 

HVAL 

32 


12 

FORMAT! 1H+ .56X.28H PHYSICAL PLANE MODE NO., 13/) 

HVAL 

33 


14 

F0RMATC1H+.58X.31H TRANSFORMED PLANE MODE NO., 13/) 

HVAL 

34 



DO 170 I = T,L 

UVSL 

35 



JL=ML (NEW,1 > 

HVAL 

36 



IF IJLi 110,170,110 

HVAL 

37 


iio 

HR ITE ( IH , 201 I 

HVAL 

36 



JLP* JL/2 

HVAL 

39 



IFIJL— 2*JLP.NE .0) JLP*JLP*1 

HVAL 

40 



DO it>0 J* 1 , JL P 

HVAL 

41 



S1=SQRT(S(1,J.I )*S(1,J,I)*S(2,J»I)4S{2,J,I) ) 

HVAL 

42 



IF(SliGT.O.C) GO TO 120 

HVAL 

43 



S 2=0.0 

HVAL 




GC TO 130 

HVAL 

45 


120 

S2=57i2 < 35?S*ATAN2iS(2,J,Ilj5{l,J»lJ) 

HVAL 

46 


130 

J1=J* JLP 

HVAL 

47 



IF ( Jl.LE. JL) GC TO 150 

HVAL 

48 



Ji="0 

HVAL 

49 



$3=0.0 

HVAL 

50 



S5*0.0 

I’L'nvi 

51 



GO TO 160 

HVAL 

52 


150 

CONTINUE 

rrnm 

53 



S3=SQRTlS(l,Ji,I)5SCl,Jl,I )*SI2, Jl , I ) *S ( 2 , Jl, 1 ) > 

HVAL 

54 



IFI53.GT.0.G) GO TO 155 

HVAL 

55 



S4=0.0 

HVAL 

56 



CO TO 160 

mi 



155 

S4«5?.29578*ATAN2iSI2»Jl*It»SU„Jl,l)! 

HVAL 

5® 


160 

rf'R’ITfc { I H * ^ f> 3 J $S(1» J»n > $ t 2 * J* I > *5 !?$'?* J'1*S ( )*St2»Jl*'l)*53«S4 

HVAL 

59 


170 

CONTTNur" 

HVAL 

. 60 

c 


THESE A*t ThE CPHASHES 

HVAL 

61 

c 



HVAL 

' 62 


2-60 

CCNTINUE 

HVAL 

6 3. 



RETURN 

HVAL 

65 i 


20 

FORMAT! 1H01-5X, 22, l.H-TH SPANHISE COLUMN)- 

HVAL 

65 


"2 5“ 

FORMAT (1H ,5X ,2(2Y, I 3, 1P4E13.5) ) 

HVAL 

66 



END 

HVAL 

67 
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ORIGINAL' FAG® » 
OF PQ0R QUALITY 


SUBR OUT I HE BQXP(XX,YY.XTDS«YTDG, XTEVXLftA* A R , DAH ,T ,EOO,S,CK»Q,T£ B OX P 

L 1 1 PRINT t HL7 »HLC , AM A *tJ*P 

l , JHAX ,L,HL«I1LM.)I1,HEM,HST,IM,HB,HC.HE ,HQ 1 fcQXP^ 

ICHECK-<T- PQlENfTIAL AS COMPUTED IM ;80XP» ^ -‘-SLUE 

ICHECK-1 -HUH igjgjW E8€E CORRECTION ' »CXP- 

ICHECK=2 - PRESSURE CGEEFICJEllf CORP^TCO FROM CORRECTED POTENTIAL MJXf 
I CHECK-3 - PRESSURE COEFFICIENT COMPUTED FROM TRANSFORMED *&XF 


C 

c_ 

C 

C 


3 

4 


5 

6 

~T 


-c- 


FOR PHYSICAL MING 


DIMENSION XX(n,YYU.>,XTDMl),Y¥DCtH*XTEI2,l>,XtE<2-,l»«A<2»WC,Ha» 

1 , AR < PE ,HB, 2 ! »EDG (2, 1 ( » S (2 ,NB»MB > ,T C HE,t»E»l I *NST l II 

2 ♦ HL(2,H.HLH(2,I).HLTt2,2,l),HLC<2»2,»B),.)HAX«2) 


BOXP 
8 GXP 


BOXP 

BOXf 


9 

10 

11 

12 


,AMA(MB*H8,1) 
1PRIN-IPRINT 


BOXP 

BOXP 


13 

14 




KRITE=JPRIN/1CCOO 

BOXP 

15 



IPRIN-IPRIN— 100Q0*KRITE 

BOXP 

16 



JRITE-IPRIN/100 

BOXP 

17 



IPRIN-IPRIN— lOO^JRITE 

BOXP 

18 



IR ITE-IPRIN 

BOXP 

19 



IFINEW.EQ.l) CO TO 140 

BOXP 

20 



KRITE-KRITE/IC 

BOXP 

21 



JRITE-JRITE/10 

BOXP 

22 



IRITE-IRITE/1C 

BOXP 

23 



GO TO 150 

BOXP 

24 


140 

KRIT£=KRITE-(KRITE/i0)*10 

BOXP 

25 



JRlTE*JRITE-( JRITE/10)*LG 

BOXP 

26 



IRITE-IRITE— (IRITE/10 1*10 

BOXP 

27 


150 

CONTINUE 

BOXP 

28 



NSH0S=0 

BOXP 

29 



DH- 0.5*0 

BOXP 

30 

c 



BOXP 

31 



CALL £QXPaiXX,YY,XTDC,YT&G,XTE,A,AR,DAN,T,S,CK,D 

BOXP 

32 



,IW,L,MB,HC,HL,HLC,KLVf,HL,NEW,NST, JRITE J 

BfiXP 

33 

c 



BOXP 

34 

c 


"box P'0' computes The' “potential - valuTs in each box". 

BOXP 

35 

c 


THEY ARE STORED IN THE ARRAY S. 

BOXP 

36 

c ' 



BOXP 

37 



1 CHE CK = 0 

BOXP 

38 



PEH-NEH 

BOXP 

39 


* 

JKAK= JMAXCNEH I 

BQXP 

40 



I F C UITE.EQ.O) GO TO 270 

BOXP 

41 



WR1TF C IK, 30) 

BOXP 

42 


200 

Du 205 1=1. L 

BOXP 

43 


-- 

3C 205 J=1,JMAX 
DC 205 K=1 .2 

BOXP 

BOXP 

44 

45 


205 

T« J,I,KJ=S(K.J,I) 

BOXP 

•46 


215 

IMMW.EQ.l) PkITEI I~H, j2 ) 11 

0OXP 

47 



IF (MEh.EQ.E) »-RITttIrf,J4) 11 

BOXP 

48 



ECRiATl Hl.lCX.c L-'PCTE'sT i AL C AL LUL AT t D (REAL, IMAGINARY, 

ABSOLUTE, BOXP 

49 

- 

‘“32 

> Phi'.'ic mNGLL* ) ) 

f CRiATl in. , 71X ,2bM P-?-,IC..i_ PlAhE . V| CDL NO. ,137) 

BOXP 

BOXP 

50 

51 


3a 

FORHAT ( 1H« , / 1> , J1H F L A L NODE NO. ,13/)’ 

BOXP 

52 

r 


PRINT-CUT 

BOX P 

53 


21b 

UC 250 1=1, l 

BOXP 

54 



JL = *L IHL I ) 

euxp 

55 



IF (JL) 2CD, 21 L . 22 j 

PGXP 



0 

KRITs. (IR,2L) I 

BOXP 

57 



J LP= JL / 2 

BOX? 

58 



IF I JL-2PJLP.NE.0I JlP*JLP-*I 

BOXP 

59 



LO 24 0 J^r.JLP 

BOXP 




Sl = SORT(T( J,I ,1)*TI J, I ,1)*T< J,I ,2 J *Tt J, I ,2 1 ) 
IFISl.NE .0.0) CO TO’ 224 

BOXP 

BOXP 

61 

62 



S 2 = 0. 0 

BQXP 

63 



GO TO 22b 

BOXP 

64 


224 

CONTINUE 

BOXP 

65 



S2 = 5 7 .29578* AT AN2(T(J,I ,2) ,T (J,i ,1J) 

BOXP 

~6lT 


22b 

r - SIT * *.»« ir 

BOXP 

67 



J1=J«JLP 

BOXP 

6f 

— 

— 

I F ( J1.LL-.Jl ) GO TO 230 

ji=o 

BOXP 

BOXP' 

69 
“ 70 


105 




53=0.0 

S4-0.0 

tc ji«i»u«o.o 

T'rir»'3V2 >=&".£» 

GO 70 240 


BOXP 

ROW 

BOXP 

BOW 

BOXP 

~ ROW 

BOXP 

HOW 

B OXP 

“BOW 

BOXP 

BOW 

BOXP 

BOXP 

BOXP 


71 


T3Q CONTINUE 

S3*SQRT(T<J1*1»1IPTCJ1»I»1I*HJ1*I»2I*TIJ1»I»2I> 

IM53.NE.0.01 GO TO 27* 

S4-0.0 

nnr 


234 


■GB'nre 236 
CONTINUE 


L,TT2777TiT7I ,11 )' 

236 CONTINUE 

"240 - ST5ITE <T W i 2 5 T~3 , T < J » 1 , i ) » 1 1 J, IV2T7S1.S2' 


• J1,T< Jl.I. 1> ,TC Jl.I, 21, 53,54 


tpt 

85 


m 


BOXP 

“W 

07 

ZTi 

BCUNTI Nut 

IF 1 ICHECK .GE * 1 ) GO TO 400 

BOXP 

BOXP 

nr 

89 

c 


BOOT 

90 


DO 38C J=1,JMAK 

BOXP 

91 


‘ ■TF'THLTTWW ,Tm .T0T01 CTTTCT 7 8 5"' ' " 

BOXP 

92 


I1=HLT(N£W,1, J> 

BOXP 

93 

' ■ 

T2=HLT ( SEH -T2 ,’37~ ' “ “ 

BQXP 

94 


53=0.0 

BOXP 

95 


K«I 


ITTaETTTI ,NfcH).GE. 1. OT”GO T0'~2VU~ 

280 CONTINUE 

TBTT^CT — 



97 
" 98 
99 

~nnr 


290 IF(K.FO.l) S3=0 BOXP 101 




$3*XXIK J— Xtt < NEW , J 1— DH*53 

DO 370 1=1. L 


.BOXP 

BOXP 

102 

103 

104 

105 



S1=(XXU}-XLE(NEH,JH/S3 
IFCS1J 300.34C.310 


"BOW 

BOXP 


e£E 

"51=0.0 


HBDXP 




GO TO 340 


BOXP 

107 


310 

HiHWH'i iiiii mmiiimi i i 



100 


320 

S1=S1*(XX(I J*XLECNEH,Jl ) / « S3*2. *XLE CHEW , J > ) 

BOXP 

109 



IF (Sl.GT .1.01 51=1.0 


BOXP 

— nr 



51=E0G(NEH»JMSQRTIS1> 


BOXP 

in 


”340 

DO 350 K»l,2 


"BOXP 

"ITT 


350 

s(K*j.n=5i4S<K»i,n 


BOXP 

113 


370 

CONTINUE 


now — 

— rxv 


380 

CONTINUE 


BOXP 

115 



1CHE CK = I 


BOXP 

116 



1FURITE) 390, 400, 390 


BOXP 

TIT 


39C 

Hf ITE ( 1H.36) 


BOXP 

118 


36 

F0FM4T(lHl,10X,t.lHP r jTcNTI AL 

COkRECTEO (REAL, IMAGINARY , ABSOLUTE, 

BOXP 

119 



PhASE ANGLE >) 


BOXP 

120 

c 


TRANSFER in TH P-tlNT-JuT OF 

CORRECTED VELOCITY POTENTIAL 

BOXP 

121 



Gl' T-J 200 


BOXP 

122 


4 00 

CONTINUE 


BOXP 

123 



IF (KRlTt.EL.O) GO TO 700 


BOXP 

124 



IF* ICHECK. Gt.C) GO TO 620 


BOXP 

125 

C 


CALCULATE AND PRINT PRESSUki 

COEFFICIENT 

BOXP 

126 



ICHtCk=2 


BOXP 

127 



jC -'10 I = 1 , L 


BOXP 

128 



JL = NL (HER, I ) 

. 

BOXP 

129 



uC 4J0 J=1,JL 


BOXP 

130 



T ( J , I ,31=5(1, J. I) 


BOXP ' 

131 


410 

T l J.l .4 » = 5 (2 ,J,I ) 


BOXP 

132 


420 

J NAX = JN ax { MEN > 


BOXP 

133 



CC 43C I =1 ,L 


BOXP 

134 



Jl^L ( rt r Hf l J 


BOXP 

1T5 



CO 43C J=1,JL 


BOXP 

136 



CO 430 8 = 1,2 


BOXP 

137 


430 

T < J,I ,K >=0.0 


BOXP 

138 



00 63C J*l,JMAK 


BOXP' 

139 

• » 


I l = *1LT(MtW, hJ J_ 

_ _ . . . . . „ . _ _ ^ . 

BOXP 

140 


106 









:2=NLTCMEh»2*,J )-ttLH«LiC ,1 ,J)*1 
11=11 
IK’lZ+l 

Till 6,5)=XLfc CPfcW,^ _ 

TCI, 8,5)=XL£'<HtV,J)' 

T C 1,10 ,31 = 0,0 
TC1,12,5)=0.0 

IF CA3S(XLElH6ti»J>-XXCIl) ).GT.l.E-05> CO TO 445 
C AOJUSTNENT — LEADING EDGE AND FIRST- BOX COINCIDES 
11-11*1 _ 

IK=IK-1 
445 CONTINUE 

DG 450 I -=2, IK " 

T(I,_b,5)^XX(JI) _ 

tci, e,s)-xxcn'i 

T C I,tO ,5)=TC J ,. 1I,3) _ 

T*(fl,12.5)=TC J,fl,41 

TCI-1,1,5) -XX C II ) 

"450 II«II*1 

CALL SPLN1 C I 2 *T( 1, 1,5), TCI, 2, 5), TCI, 4,5) 

i ,Ik,tCl, 6,5) ,TC1,I0,5),TC1,14,5I, NSNOS I 

C ALL S PLN 1 ( 12 ,T Cl , 1,5),TC1, 3, 5), TCI , 5,5 > 

S , IK, TCI, 8, 5), TCI, 12, 5), TC l~Tl4,5) , NSNOS") 



BOXP 141 

BOXP 142 

BOXP' ' 143 
BOXP ‘ 144 
"BOXP 145 

BOXP 146 

“ScDfF 147 

BOXP 148 

BOXP - ’ ‘ iV«J 
BOXP 150 

BOXP ‘ 151" 

BOXP 152 


BOXP 

BOXP 

BOXP 


BOXP 

- T5T 

BOXP 

158 

BOXP 

159 

BOXP 

160 

BOXP 

-161 

BOXP 

162 

BOXP 

163 

BOXP 

164 



1 k E t, k i „ .-CltN»lAi_ } t C T I CN 

TO- Ob TAIN 

BOXP 

211 


V,*L'i'-> P - £ 5S«J» i C-itr e ICI ■ NT CCnOjTATICh 

gn physical wing 

BOXP 

212 


GC T' ‘•it 

; 

BOXP 

213 

700 

'CCNtlNUt. 

“ 

BOXP 

214 


kUU«H • 


BOXP 

215 

' 20 

T r r V « T ( X-"C ‘3 * , 1 2 , 1 "y rf- T rt spa.n-ISE COLuHn) 


ffUXF 

21fc~ 

25 

t !■>«! Mri ,'t-A L."tfc 13.5)) 


BOXP 

21? 


END" 


" " BOXP ' ' 

218 


108 



ORIGINAL BAJj-' * 
OP POOR GgiALITV 



subroutine bqxpoixx,yy,xtdg,ytog,xte»a»ar,oah,t-.s»ciud 


soxpo 

2- 


$ ,TW.L«H8,HC,HL»HLC,MLW,N1,NE»I,I»ST«ARITE» 


BOXPO- 


c 

SOLUTION OF SIMULTANEOUS EQUATIONS FOR THE POTENTIAL 


eoxpo 

. 4 


DIMENSION XXtlliYYI i),XTDCUI.YTD60l*XT£<2*lI,A«2,HCVl*, 


eoxpa 

S>. 


1 ,ARIH8,HB,1)*TI2«HB*1)*S42,MB*1) 


BOXfftj 

<6 


2 ,MLI2,lI,HLC<2,2*ll*m.HC2*ll»NST41l 


BOXPO.- 

7 


DH»0.5*D 


BQXPthr 

-a 


IF ( JRITE.EQ.O) GO TO 25 


BOXPO 

9 


IFCNEW.EQ.2.0R.N1.GT.U GO TO 25 


BOXPO 


c 

PRINT INFLUENCE COEFFICIENT 


Boxrd 

u 


WRITE IIW,10C) 


BOXPO 

12 


JL=2«HLINfcW*L) 


BOXPO 

13 


K-JL/4 


BOXPO 

19 


IF( JL-4*K«NE.O) K“K+1 


BOXPO 

15 


00 20 I»i,L 


BOXPO 

It 


I 1*1-1 


BOXPO 

17 


WRITE <IW,110> 11 


BOXPO 

18 


DO 15 J1=*1»K 


BOXPO 

19 


15 WRITE IIW»12fl ) C(J,A(1*J-»I)»AI2,J»IJ)»J*J1,JL»K) 


BOXPO 

20 


20 CONTINUE 


BOXPO 

21 


25 CONTINUE 


BOXPO 

22 


Il-NSTCNEH) 


BOXPO 

23 


ANI-1.0 


BOXPO 

24 


00 26 1*1,11 


BOXPO 

25 


26 AhI*AMINl(XTDG(I) »AMI ) 


BOXPO 

2* 

c 



80XP0 

22 


IF< JRITE.EQ.O) GG TO 28 


BOXPO 

28 


WRITE! IN, 168) 


BOXPO 

29 


28 CONTINUE 


BOXPO 

30 

" 

”'11=0" 


BOXP'iJ 

31 


NFLNS = u An 


BOXPO 

32 


00 90 1*1, L 


BOXPO 

3T 


X=XX<I) 


BOXPO 

34 

T~ 

ADJUST UPSTR t A 1” INFLUENCE 


BOXPO 

35 


KO = 1 


BOXPO 

36 


If <NFLN5 .eG.O) GO TJ 30 


BOXPO 

37 


KO = NAXCll,I-NFLNS*l) 


BOXPO 

38 


30 CCNTlNUc 


BOXPO 

39 


JL=NLInEh,I) 


BOXPG 

40 


IF < JL.fc'O.Cl L-C TJ 90 


BOXPO- 

41 

c 

DEFINE >, I NG A WAKE liJXtS 


BOXPO 

42 

c 

Jw=NU. OF nAKc BOXES IN 

ROW 

BOXPO 

43 

c 

Jc-Nt* Cf" HlN*o 9CXHS 


BOXPO 

44 

' c 

. -- J i =C F L. C » OF F I 6 S T hInG 

BOX 

BOXPO 

45 

c 

JN=C=.Gt' cf LAST WING 

BOX 

BOXPO 

46 

c 

J 5 L* o t *v of - FIRST HAKE 

BOX 

BOXPO 

47 

c 

J NH- i\J c * JF LAST xAKE 

BOX 

BOXPO 

46 


JN = 0 


BOXPO 

49 


If (>.L;.Jri ) CL TO 3t 


BOXPO 

50 


Tl 


BOXPC 

51 


32 I f ( j« jT« JL J ul f'J 


oOXPQ 

52 


IF 1 *.LT. I> T. ll^-. , J )* ,-t) > 0. ’ ‘ 


fcQXPO 

53 


J h x J»* ■» I 


BOXPO 

54 


33 J=J*1 


BOXPU 

65 


L _ TO 32 


BOXPO 

56 


3* Jh = 


BOXPO 

67 


u ) Ja 


BOXPO 

56 


35 JS * JW-»1 


BOXPO 

59 


JSW=1"~ ■ — " " 


BOXPO 

fcC 


GO TO 37 


BOXPO 

61 


36 JS =1 


BOXPO 

62 


JSW*JE +1 


80XP0 

63 


37 JK = JS VjE— ! • 


BOXPO 

64 


JNN* JS « ♦ J W — 1 


BOXPO 

65 


IF iTi.EQ.oT5I! fo jo 


"BIKPO ‘ 

“ "6‘6~ 

c 

SUBTRACTION OF CbNTRIfcUT IONS OF PRECEDING ROWS TO UPWASH 


BOXPO 

67 


00 4? J*JS» JN 


BOXPO 

68 


CC 45 K= K 0 , 1 1 


BOXPO 

69 


KL*MLCNEW,K>' 


BOXPO 

70 


109 



:jOU 


BOXPO 


71 


M=I*1-K 

f'F JKL.ECi; 0) C-5 to *5 

-'C J "* *' 1 > AL ^ 

N 1 = N+ J 

nz= 3 At<s <n- j ) * ; 

1 (f 1 .' t 

Ai=A(2,Nl,Kl)*A(2i 'C.<n 

rfT-1.15 ” 

! F (N.GE .rLC < ft i- > . AnKPU'IHF 4 » 2 «<< } „t!E .0 } HT=AR (N«JK«NEHI 
SI1,J,U*SU, j,Ii-SAi-iU,h,X>-A i »StZ,N 4 K}j*Ht 
*iP i (<r ,J>i )*5 1 «r i J , J >- i-'-/*., i i 5 > A 1*5 1 2 »'•*<) J *WT 
55 CONTINUE 

5 7 CONTINUE 

SEtf JnC""uF Ti7K“:x > j l Puli XNt 3US E^UAf* CHS 
50 Jl‘ >i. J = liJi. 

DC K--i,jl 
Nl=j4> 

M = I AbS { j-r )*1 

KT=1.C 

I F CR. (7E".SLC(r4-_ Kih , I l 4 556 « rlLlH HtVi i > II «NETO 1 ~VT=AmK'il*'HtM'>'' 


BOXPO 73 

BOXPO 74* 

BOXPO 74 

" BtfXFO — 7* 

BOXPO 77 

BOXPO 78 

BOXPO 79 

UBxFB “80 

BOXPO 81 

“BOXPO 52 

BOXPO 83 

— BOXPO" 3'V 

BOXPO 85 

* Boxpir oo 

BOXPO 67 


IFIJ.cC. v ) fi 7 = 1 » 

T = j+. 

T ( Z , J,<) = { A< 2 ,M , I ! » 

bi cuhttnue — * 


i(J, 1 £KT 


SUBTRACTION OF CONTRIBUTION 
S t»I.1G5-T{HAKE>*PHHMAKE> 

IF'( JW.E J“.0) CC TO bu 
DC So J=JS,JN 

ITC " T5 R=TS ff» j K V •' ~ 

55 S (2 » J 1 1 f=S fEtj, [ l-'d.j/K! r2,'^»l i -T < 2 S JiSjSsTi ,N,fI 
3b CONTINUE 

IF < NL V C N E bVf J • L T . 1 7 OuYJ 60“ 

RC--POS IT I ON ELEMENTS OF 7 


NH=N» JN 

' OtT 5”o 7= iTTt 

M H = J ♦ J H 

' Do 55 k=1,2 ' 

53 T(K,J,M=TIK,Pb,NH) 


1 

SCLUTir.N Oh ; CLf T _ _ . 

r b bCNTINbt 

Il(J*ITt.'. -.v.) 0^ T.. ?_ 
xnIJI ( 1 1 * i lb?! I 
.1 > 70 X=;.2 

Nl=l . 

I- £ - " l *~3 


BOXPO 

BOXPO 

BOXPO 

BOXPO" 

BOXPO 

BOXPO 

BOXPO 

BOXPO 

““BOXPO - 

BOXPO 

BOXPO 
6 OXPO 
BOXPO 
BOXPO 
BOXPO 


BOXPO 

BOXPO 

BOXPO 

BOXPO 


BOXPO 

BOXPO 

BOXPO 

BOXPO 

BOXP O 

“BOXPO 

BOXPO 

BOXPO 

BOXPO 

BOXPO 


97 

— W 
99 
100 “ 
101 



. J z7 J-=* * J' 


80XP.0 

124 

67 

■Ka JTt ( 1^4 lc^) l j 1 « N » 

T i - 1 J 1 A 1 t .N = N 1 .hi > 

BOXPO 

125 


1 ( ( '*. • JL J V/C T '• 


BOXPO 

126 


L = i 


BOXPO 

127 


N _ = N I 


BOXPO 

128 


\ 1 T’ fit 


BOXPO 

129 

A>* 

U NT!'*''" 


BOXPO 

130 


*'t 6*5 J - 1 1 J i 


BOXPO 

131 


ITl l N,17J "iJt* 

Oi'.iJtll 

BOXPO 

132 

r y 

c:ntinu^ 


BOXPO 

133 

7 r 

Ct MI <uc 


BOXPO 

134 

167 

h C p V '1 T ( l* 4 -' 4 '-* /■ 4 I A 1 1 7 

-I. SPA-.-! sc COLUMN) 

BOXPO 

135 

1 6b 

f 0-r r ‘± T ( if 1 i * i C> i nC 

-aTFIX (A J .uF 

IA)* IX 1=481 FOR VELO BOXPO 

136 

J 

CITY PuT‘M in u) 

- - I.'.o is. - C U..LINN/ / ) 

BOXPO • 

137 

I b5 

eOhnA T ( 3X»MZ*«Ir(» 

id 1 .Mt ) 1 ’ tl‘! ) 

BOXPO 

138 

i 7C 

f-QRfJlT ( 9Xicn\ l ,iM 

3nii..i,, 12 ,2M} = fE13.b) 

BOXPO 

139 

?; 

v. JmTI . 

. 

BOXPO 

440 


no 


ORIGINAL': BAG® «g 

op «*» wS? 


aoxp.o 

K x ilSJ'JlEC ( f*c t Jt: f if T i j 1 1 1 t i ] i BOXPQ 

IFl^.ht.D' CD If 4 “ ’ "BOXPO 

.COMPUTE HAKE POTENTIALS — BOXPO _ 
PHJ t TE J*EXPI— IK*<X— XTEIl" BOXPO 

IFIX.Lt .tAMI-U) uU TO 36 BOXPO 

3s-i ‘ - - - "esxps 

Y-VYCJS_) _ _ BOXPO 

IFfX.LT.tXTt 1N£K t J 3 )— OH I ) GO TO iZ ' BOXPlff 

IF < x_»GE • (XT E I NeWt JS )»OH ) I GO TO 82 BOXPO 

IMX.GfixTtllitXi JSJ) GO TG 76 • " BOXPO~ 

f TR=S <1 » JS « 1) BOXPO 


P7I~S(2*JS*I) 

kK=J - BOXPO 

XB=XVo " BOXPO 

GO TO _77_ BOXPO 

P«IX-XTFCNEM*JSJI/D " " BOXPO 

ptk=( i„o-p)*s ti,js,n*p*sa,js, i-i> boxpq 

PTI = U.O-P)*SC2,JS,n*PPS(Z,JS,I-l)' " ’ " BOXPir 

K K*0 8 0XP0 

“xF^X “ BOXPO 

CONTINUE BOXPO 

IFIXB.GT.i.C) GO' to' <J2 ' ' " BOXPO 

XH=(XB-XTE(h£h,JS))«CK BOXPO 


IKH=I*KK BOXPO 

S<l,JS*IKH)=PTR*COS<XW) »PTI*SIH«XH) BOXPO 


SI2»JS,IKH)=PTI*COS(XW)-PTR*SIN<XW) BOXPO 

K K^KK *1 _ BOXPO 

“ xe'-X B*D ' " * “ ' BOXPO 

GO TO 80 BOXPO 


JS-JS-»1 BOXPO 

IFC JS^LE. JLI G£ TO 75 _ _BOXPO 

CONTINUE ’ BOXPO 

11=1 L+i BOXPO 

RETURN BOXPO 

WRITE (IW,140) BOXPO 

STOP BOXPO 

FCKNATCIHi»2C>i6?HINFLtJ!:NCL COEFFICIENT (REAL AND IMAGINARY!///) BOXPO 
FORMAT t, IhO, i a i! IX 3 E P / R A T IvJN IN CHORDhISE DIRECTION) BOXPO 

FGRMATUh ,Fa ,M:A,IP2e13.6| ) BOXPO 

FC=MATI 1HCIGXi50-iSC_uTICN OF SIMULTANEOUS LOCATIONS FOR THE POTENT BOXPO 
IIAL FAILED BOXPO 

;NC BOXPO 


154 

155 

156 

157 

158 
~TW 
160 

161 

162 

163 

164 
S' 
6 


167 

168 

169 

170 


111 



c 

T" 


'biwhisicw hlt i 2 ,i . n»> .xUiaati* ,XTg<?,»t>yVXiiiy4 «yy< w »i 

1 «SFOXIH8*«0**Sf O^§8#HO*M*SFeH<IW»to*lti>HSf0<^S5 

2 aHMia.ittahtiit»l«Tiux{a!i ' 


?55K:l---^3 


M 

MO SMOOTHING IH SUtROOflBE 51*30 TH iSFLMit 




'-W 




AM2-1.0 . . 

NZ-KSF0IN21 


H 

12 



IFIKSf0IM2&.£0.0.ANDoN£*.E0.2l HZ*KSFS 


1% 



08 150 J-lWRAK 

. 

14 




fSccT 

ii 




F9KI 

16 



NC«NZ*1 

FSST 

17 



Kl*l 

FORCI 

It 



K2-4 

FWCI 

19 



KC-0 

FORCI 

20 



Ttl,l)-XTEIHEH,J) 

Fft*Cl 

ii 


40 

KK-HLT1 

FORCI 

22 



III, 2)-XLEINEH,JI 

FORCI 

23 



Tl 1,K2)»0.0 

FORCI 

24 



KA«2 

FORCI 

25 



KB-NZ 

FORCI 

26 



IFIK1.EO.3.0R.NZ.EQ.21 GO TO 45 

FORCI 

27 



IFItCC.EO.2l GO TO 45 

FORCI 

28 

c 


EXCLUDE LEAOINC EDGE POINT IF IT IS TOO CLOSE TO FIRST BOX 

FORCI 

29 



IF ( XX IKK J— XLE INEW»J)«GT. 0.2*01 GO TO 45 

FORCI 

30 



KA-1 

FORCI 

31 



KB-NZ-1 

FORCI 

32 



KC-1 

FORCI 

33 


45 

CONTINUE 

FORCI 

34 

’ - 


00 50 K-KA'KB 

"FORCI 

35 



T IK, 2t=XX<KK) 

FORCI 

36 



IF ( Kl.EQ. 3 ) GC TO 50 

FORCI 

37 



T(K,K2)=S(K1, J,KK) 

FORCI 

38 



IF (NEM. EO. 1. OR. Kl.EQ. 11 GO TO 50 

FORCI 

39 



TIX,12>-AHAIJ,KK,NEH) 

FORCI 

40 


50 

KK«KK*1 

FORCI 

41 



IF (KC.EQ.2 ) GC TO 51 

FORCI 

42 



IFIK1.E0.3) GC- TO 60 

FORCI 

43 



CALL SPLN1C 1, Til, It ,T|1 ,B>,T(1, 91 ,KB,T( 1, 2 1 * T 1 1,K2 ) * T t 1, 10) ,N1 

FORCI 

44 



TINQ,K2)-T(l,et 

FORCI 

45 



IFCKC.NE.lt GC TO 51 

FORCI 

46 



aC*2 

FORCI 

47 



GO TO 40 

FORCI 

48 


51 

CONTINUE 

FORCI 

49 

c 



FORCI 

50 

"' r “* ' 


IFCN2.GT.21 GC To" 52 

FORCI 

51 

c 


ADJUSTMENT FOR TRAILING EDGE VELOCITY POTENTIAL 

FORCI 

52 



KK-MLT'(N£W,2,j) 

FORCI 

53 



DUM=SQRTIIXTEINEW,J1-XLEINEV, J1 1/ < XXC KK l-XLE INEW, J 1 1 ) 

FORCI 

54 



TIN0,K21=DUM*T CNZ.K2 ) 

FORCI 

55 


52 

CONTINUE 

FORCI 

56 

c 



FCRCI 

57 



KC-0 

FORCI 

58 



IFIK1.E0.2) GC TO 55 

FORCI 

59 



X 1*2 

' FORCT - 

60 



K2«6 

FORCI 

61 



CO' "TO 40 " ' 

FORCT 

62 


55 

K 1 = 3 

FORCI 

63 


' 

Go To" 40 

FORCI" 

64 


60 

TINQ, 21-XTE (NEH, J1 

FORCI 

65 



T INd, i2)=CTfc< J ) “ " 

FORCI 

66 



Tl 1,1Z1*1.C 

FORCI 

67 



IFIABSIXL£INEh,J t—XXI MLT 11 l.GT. l.E-05 1 GO TO 67 

FORCI 

68 

c 


LEADING EDGE ANO FIRST BOX COINCIDE 

rcRCi 

69 

■ 


NZ=NZ-1 ~ ~ “ 

FORCI 

70 


112 









ORIGINAL* BAGS Ife 
OF POOR QUALITY 




HQ*MO-l 

FORCI 

71 



DO 65 X*7,H0 

Tmu — 

— rr 



TIK, 2)»KK*i, 2 ) 

FORM 

73 



T t K * (Kt^'y ^ J 

""FORCI 

TV 



TtK, 61-TIK+l, b I 

FORM 

75 



IFtNEH.Ea.il efiTff 65 

FORM 

7* 



TCK,12»=*TIK*1,12» 

FORM 

77 


03 CUN T INUE 

FORM 

— nr 


bl 

CONTINUE 

FORM 

79 

c 



FORM 

-- *w 

c 


interpolate deflection at aox center* leading and trailing edges 

FORM 

81 



CALL SURFZtYYtJ),T(i, 2 1 * 1 .NO ,2 *T (1* 1 31 ,1 1 1, 141 *OUH,SfOX< 1, N2 » 

“Wi 

82 



1 ,SF£Y(l,n 2 ,N£W) ,SFDH(1,M2*NEH),NZ,2) 

FORCI 

83 

t 


PERFORM CH0RDRI5E InYeGRAT lON 

FORM 

— “SV 



00 75 K* 1 i NC 

FORM 

85 



I F t REH . Elf. Z7 AN2*TU«lZT« r TTir»l2 _ )' ~ 

FORM 

"86 



T(K,15)~(T(K,14)*T(K,4)+CK*T<K,l3)*T(K,bn/AH2 

FORCI 

87 


lb 

TtR.l5)-=tTtK, 16}"*TTi< ,67-CK*T (X , 1 TP*T t K , aTT7aT42 

“FORM 

“ 88 “ 



CALL iNTGLtTt l,2),Tll,15),T( J,17),TC J,18>,Ttl,19J,2,Ntl,NZ»NHI 

FORM 

89 




FTfKCT 

90" 



DUH=TlNQ,ni/AM 2 

FORM 

91' 



r t j , i ?7 =ri *c ; ^T^DuM^n r; nr ~ 

' FORM 

— 92 


150 

T(J,ld)«TINO,fc )*OUM~T I J » 1 8 ) 

FORCI 

93 


CON T I NUE 

FDRC 1 

- 7} 4 

_c_ 


PERFORM SPANHlSE integration 

FORCI 

95 




mbs 

96" 



Kl = 17 

FORM 

97 



K2«I5 

FORCI 

“98“ 



Tii,n»o.o 

FORM 

99 


165 DO 170 K-l, JMAK 

FORCI 

■ — loo 


170 

TIK, 2 )=YY CK ) 

FORCI 

101 



“NZ-JMAK 


VfM 



CALL SPLHlt 1 ,T (1,1},TC1,S),T(1,9) ,NZ ,7 1 1 ,2) ,TC 1,K1) ,T( 1,10) ,Nt 

FORCI 

103 



T (NQ,K2 1 =0.0 

FORCI 

TL04 



TC 1,K2)-TC1,8) 

FORCI 

105 



00 175 K*1,JMAK 


— 106 


175 

TtK+l*K2)“T CK,K1) 

FORCI 

107 




FORM 

108“ 



Kl-16 

FORM 

109 



K2-16 

■ BK 1 

■ 110 



60 TO 165 

FORCI 

111 


180 

CONTINUE ~ — ' ' — 

FORCI 

"TTST 



NZ»J«AK*1 

FORM 

113 



T 1 1,21-0.5 ~ ' : - 

FORM 

— nr 



T CNO,2l“YNAXt NENJ, . • 

FORCI 

115 



DC 15 3 K " l » JMA X 

FORCI 

116 


1 ■* 5 

T(K*l»2>=YYIxJ 

FORCI 

" lTT 

c 


■cl No TIP CORRECTION 

FORCI 

118 



»* 1-0 

FORCI 

119 



Xl =T < h"J,. 1-T (XZ-i 

FORCI 

120 



PR=Tt62-l,l5) 

FORCI 

121 



P I=TCN7-1, it) 

FORCI 

122 



IFIT<MZ,2)*.5«L.LI.tlNu,2)) 60 TO 195 

FORCI 

123 


19C 

L'Uf^Gw^TT lTIf-<.,c)-I(NZ,Z) )/XDS 

For.ci 

124 



T (N/ . 15 > si-Ur.+ FR 

FORCI' 

125 



T [Ml, It. )xDUr« F 1 

FORCI 

126 



if :.n ut to iv<> 

FORCt 

127 


1Y3 

IFIYYC JbAa)*: ,G5*J.6t.YMAX(NfcN) ) GO TO 196 • 

FORCI 

128 



<1=1 

FORCI 

129 



n;= 6Z «i 

FORCI 

130 



■* l 

FORCI 

131 



T (NO,: ) =T 1 hu- i ,2 ) 

FORCI 

132 



T IN/ ,2 1 = T (N/- 1 ,:)*0 

FORCI 

133 



T ('* ) = 1 1^1 

FORCI 

134 



tlN5,lG)=fTNW=;,ib) ■ “ “ 

FORM 

T3’5' 



or tj 1*0 

FORCI 

136 


l -<6 

CONTINUE 

FORCI 

137 



IMAiSCTIt-0,/)-T(N/,2)).CL.l.E-05) GO TO 197 

FORCI 

138 



NC»NQ-1 

FORCI" 

139 



NZ=«7-1 

FORCI 

140 


• 113 


TCHJ-, 2 ) =T.( MU* i , 2 1 

FORCI 

141 

T(NC,l?)-T(N0-*l,l2) 

FORCI 

142 

TtM^.lM-UN.-tl.K 1 

forci 

143 

C lkIJ 

FORCI 

-144 


FORCI 

145 

CALL INTGL(T<l,2),T[l,li),T(l,n,T(2,l>,Tll,19> ,2', NO > NZ , Mrt J 

FORCI • . 

146 

KE tU*N 

~F0RCT~ - 

— 147" 

t IhC 

FORCI 

148 



r^nnn j o r* no 


ORIGINAL' RAGjE IS 
OF POOR QUALITY 


SUBROUTINE HRECIDA.T.NM.NB.KSFHtSFHX.SFrtY.SFHH.IH.IPRlNT) 

- 'SIT-THE^Uf<FACE FTT~7TF~7iACFrTiITK!rER “ 3 ' - 

QIHfcNS 1UN DA C 1 J . T < Nrt . 1) , SF NX < X > , SFnY { I ) «SEMH ( 1 ) 

t0NST*0.2S&Tl54V ' ' ' 

KSFH«DAI<77) 

'IFI'KSFH) 80. SC, 10 ' 

FITTING OF GIVEN PkESSOxETn'aCH TCT A SPLINE- SURF ACE ' 

IF(Nb-KSFHJ OC. 15^15 
CONTINUE^ 

KP = 70 1 

013 30 IP=1.KSFP 

SFHX'UP )=OA't'KP I/OAti^l ' 

SFHY(IP ) =» 0 A ( KF*l)/uAI2A) 

S F MH t TP I =0 A < Kf '♦ c > 

JA(96>=1, IisPLT DAT* Axt PRESSURE COEFFICIENT 

DAfS6J«2» InVlT CAT'a A < t lJCAl «A(,h NUflbck 
I F (DA 1 96 ) — 1 . J ) 7 5.20,25 

CCNVE k t P R t SEUR t C J c FrTC"I E N f~ T«T5 ~ L'QC A L NACH NUHBER' 

Sf HH( IP )=S u kT (5.41 1.2/1 { I..0. ?« S FttHf IP I l+RCONST 1-1.) ) 
CONTINUE 
K P =K P. A 

SPLINE-SURFACE FITTING of data 
continue 

TF I IP R 1 NT .~NE . C'J' VkTtE C"f ri > 100’) ’ 

CALL SUkFl (Nn.kSEn, T.SFnX, SF nY.SFNH.I PRINT) 

RE TURN 


PRESENTLY InPLT jf pressure COEFFICIENT in 
A PULYNOMAL fora is not allowed 
THE''F' biLChiNC m=j“^ACTr inPut "Tr** ' ttilyndh ial • 

SO CONTINUE 

kseh=3" - ' 

SFHHI 1 ) =0A ( 71 } 

■ ' S FMH n 1 LTJ aOT) Ada Cl <, ) - • — 

SFHH(3)=0.0 


75~ TPR=Tt> 

GO TO 85 

80 IP'F^97 

35 WRITE 


STOP 

loo' FoRNATriHirrm^HCMrjTFo' 
S(R<I)**2)*CALCGIRUI**2))) 
i TO “FBrnaTT iho. 10X . i-^hhrEo — SaC 
END 


TAT A. T5T 


1 PR I NT ) HREO 

HRED 
HREO- 
HRED 
HRED 
HREO 
HREO 
HRED" 
HREO 
HREO 
HRED 
HREO 
HRED 
HRED 
HRED 
HREO 
HRED 
HREO 
HRED 

ft" HRED 

. )> HREO 

HREO 
HRED 
HRED 
HRED 
HRED 
HRED 
HRED 
HRED 
HRED 
HRED 

=TSB*AI*X HREO" 

HRED 
HRED 
HRED 
HRED 
HRED 
— HREO - 
HRED- 
HRHO 
HRED 
HRED 
HREO 


HREO 

, Y« r- STJH'OF — HITl* 'HEED" 
HREO 
HREr 
HREO 
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SUBROUTINE INT GL< X , Y» V* » V I ,5 *N, NO,NZ , NK 1 

~t IHTFSTaT 10* BASfc'5'TlfrTPTlNE FUNCTION ~ "" * ’ . 

DIHENSUIN Xin,Y{Nrt,2>,S(NN,l) 
t " "OtriNE LIJI " 

DO 20 I =2 . NO 

20 sTT,u»xt i i-xcT-n 

IF INC.ECi.21 CL TO 50 

T~ DEF lNt _ T R f— "2T Too NA L. COEfFTClfNT hajrix”” 

1 9° ijL Is,2 » N ^ 

sTI.Tl* SiI,iJ/o.o 
su.Ji»(su,ii*su*x.i.n/3.o 

25 S 11,41= SflU.U/i.'J 
S( 2,21=0.0 
S<n2, 41=0.0 

C DEFINE RIGHT-FAND-S10E lOLUMN hatrix 
K = 1 

35 DO 40 1=2, NC 

40 $ 1 1 , 5 1 = I Y < l »x 1 — Y ( I — 1 ,K 1 1 /S C I » 1 1 
DO 45 1=2, f»i 

' “ '45' S ( I »fe >* -STr*T,T)^l I ,51 
C SOLVt FOR CCtFF ICIcNTs Of SPLINE FUNCTION FIJI 

CALL TRI DI t 2, N2, St 1,21,:. (1,3), SI 1,41 ,S 1 1 ,6 1 ,S tl , 7 1 , S <1 ,8 > ,S 1 1 , 91 ) 
50 CON TINU E 

S I 1,71 = OVS 
S INa, 71=0.0 

7T _ TJ _ T5 ... --- 

DO 60 1=2, NC 
VI-VIVO.’^SII ,1> 

J *(Y( I ,K1 *Y( I— 1,K 1-5 < I , 11*51 I ,11* <5 Cl ,7 1 *S U-l , 7) 1 /12 . ) 

60 “com NOE 

IF(x;EC,21 RETURN 

IF(N.EG.l) RETURN 
- K - 2 . 

GO TO 35 

TPD ~ " ~ 


INTGL 

2 

INTGL' 

3 

INTGL 

4 

INTGL 

5 

INTGL 

6 

INTGL 

7 

INTGL 

8 

INfGL 

0 

INTGL 

10 

‘INTGL 

11 

INTGL 

12 

INTGL 

13 

INTGL 

14 

INTGL 

1’5 

INTGL 

16 

INTGL 

17 

INTGL 

18 

INTGL 

19 

INTGL 

20 

INTGL 

21 

INTGL 

22 

INTGL 

23 

INTGL 

24 

INTGL 

25 

INTGL 

26 

INTGL 

27 

INTGL 

28 

INTGL 

29 

INTGL 

30 

INTGL 

31 

INTGL 

32 

INTGL’ ‘ 

33 

INTGL 

34 

INTGL 

35 

INTGL 

36 

‘ INTGL 

37 


SUBROUTINE TRI0I(K1,K3,A,8,C,D,V,E,FI 

* 



BlKENSJQH A411VB41),C(1I,D(11,V(11,E(1>,F41) 
IF Ht3.NE.Kll GO TO 5 



.4 

7tX3I*OtX3JS&itC3 1 

Horn 




5 C8|lTINy£ 
EttlMWl) 


f mate* ' 

- 

& 

FIRH*3MHli/£*KU 


-~fmc 

M 

Sf 

t rd W I=*2,K3 

E I ****!«-* m*C 1 1-1 )/£ U-l 1 


■**m r 

TRISI 

it 

u 

io F* II-tOtll-Af 114FU-11 )/E(l) 
VMC3»*FtK3J 


Tjtloi 

HUM. 

13 

14 

K2*K3-K1 
DO 20 3=3*K2 - 

■ 

— TSIBt 

T*«a„. - 

.iS- 

I-K3-M 

20 ¥II»^fM-CUM'VtI*l)/EIII 



w 

RETURN 

END 


TRJ0I 

TRIDl 

s 

» 
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■nr> c* 


ORIGINAL 1 BAGS) TS 



. OF POOR O.UArA n 



SUSP. OUT I nE SPLMlNl,X,Y,DY,N,Xx , YY ,D Y Y , NSHOS ) 

SPLN1 

2 

c 

. X , Y tDY 5 *! NTbRP CL A I ION INDEPENDENT, OfcPENGfcNT VAR I ABLES , . ANO DY/DX 

SPLN1 

3 

c 

Nl -=NC« OF IN IE RPO'_ AT ION PUiNTS 

SPLN1 

A 

c 

XX, YY = I NPUT InuLPcNuEnT ANt, 'IL P EnDcN T VAo 1 AtLbS 

SPLN1 

5 

c 

N =NG. OF INPUT POINTS 

SPLN1 

6 

c 

OY Y = ii’tYY I/l(xa) FDrt Input Data 

SPLN1 

7 

c 

NSNDS =CunTRULS uF jrtOOTHlNG A^C PRC- InTEkPOLATION 

SPLN1 

8 


D I "EnSIGn XU >»Ytl> .DYtll »XX( 11 ,YY 1 11 ,OYY< 11 

SPLN1 

9 


N2=N 

SPLNL 

10 


ipre'=nsngs/ io ' 

SPLN1 

11 


NShC=NSnUS-IPktVi3 

SPLN1 

12 


I F ( I"Kt.Lo.U Gu Tu _>0 

SPLN1 

13 

c 

STDRc INPUT DATA P _K PRt— INT dKPJLATIJN 

SPLN1 

l A 


DC 23 J = 1,L- 

SPLN1 

15 


YtJl=XX( J) 

SPLN1 

lfc 


2C DYtj)=YY(J) 

SPLN1 

17 


CALL 5P 1 St T I N ,aa » Y Y , u Y Y , 0 • 0, •' J 

SPLN1 

18 

c 

PRL—INTckPlLAT Iufc 

SPLN1 

19 


DC cr> 1=2, N 

SPLN1 

20 


11=1* i-2 

SPLN1 

21 


12=11*1 

SPLN1 

22 


aX(I2)=Y(I 1 

SPLN1 

23 


Yr (12 ) = „Yt 1 1 

SPLN1 

2 A 


XA( il )=-j. 3* I Y l I-i J*Y( i 1) 

SPLN1 

25 


2; (.ALL SPLNt I / a ( 1 1 i , i ,il , Y » ,-Y ,o Y Y , Y Y I I-l 1 * J‘JN» 1 ) 

SPLN1 

-26 


’* 2=2- > f-- l 

SPLN1 

27 


id LLNllNur 

SPLN1 

28 


IFIN^c.tU.) CfJ L 'll 

SPLN1 

29 

c 

ShCOTr Input D/T« aa, YY 

SPLN1 

30 


Call. ShCuTh ( It £ f a ^ , Y Y t t' » Nb ^0 } 

SPLN1 

31 

c 

INT=R°uLAU Y AT a -k,-, aa, YY, OYY AND CALCULATE DY=DIY)/DIX) 

SPLN1 

32 


NC CAtL SPlSfcTIN£, , /,YY,jYY,0,0,uJ 

SPLN1 

33 


CALL SPI.N2I X,M ,NA ,Xa, YY ,DYY,Y, UY,2) 

SPLN1 

3A 


RETURN 

SPLNI 

35 


FRG 

SPLN1 

36 


— - - — 

— 

— _ 


SUBROUTINE SPLNZtXP,NP,N»X»Y,D,SPF,SPD*K> ' 

SRLH* 

a 


BTHEHSTUS XU>,YU),<Hl*,Xm).SMPU»,'SI»&ill ““ 5 MJi 2 T 

EVALUATES * NATURAL CUBIC SPLINE AND ITS f IRST DERIVATIVE USIMC Sf Vt& .» 


SLOPE ARRAY D CALCULATED BY SPISET AND USING THE IN>UT DATA 5FOS Y 


C ARRAY 5 X ANE Y 


SPLN2 

* 

DO 10 J=1*NP 

If! XPIJI.LT. XII). QR.N.tQ.l) GO TO 6 


SHUT 

SW.N2 

' ' f 

« 

DO 2 1=2, N 


mm 

♦ 

IMXPtJI.LT.XUll GO TO A 


SPL«2 

10 

2 CCNTINUE 


SPLN2 

il 

SPF(J) = Y(N)*C(NJUXPl JJ-X(N}> 


SPL82 

12 

IF (K.EQ.l) GO TO 10 


i IT YWM& 

li 

SPOt J)=0(N) 


SPLH2 

1.4 

GC TO 10 



I5~ 

a ci*i./< xt'i )-x t i-m 


SPL82 

16 

! C2 = x( l I-XP t J I 


5>Ln2 

17 

: C3«xPt Ji-xti-i) 


SNIH2 

IS 

! " CA=C2*C 1 


SPLN2 

14 

C5=C3*Cl 


5PLH2 

20 

SPFt JI=C5*C5*I 11.+2.PCAJ *Y(I »-CZ*D( II 1 


_ SPCN2 

21 

: i ■♦CAACA’M tl.*2.*C5)*YU— 11 +C3*D< I— 1 J 1 


SPLH2 

22 

IFtK.EC.il GO TO 10 


5102 

23~ 

C6-2.*C2-C3 


SPL82 

2-4 

C?=2.*C3-C2 


SPLH2 

25 

SPDtJ)=Cl*ClAtC3*t2.*t l.+Cl*C6I*YU 1-C6*0tl>» 


SPL82 

26 

A -C2> 1 2 1 1. *C 1$C? 1 *Y < 1 -1 1*C7 *D t I-1III 


5PC82 

27 

GO TU 10 


SPLN2 

28 

6 sPmi=?TTr-cmiA(xt ij-xpcji i " "t 

“ " 

~ SPLN2 ~ 

29 

IFtK_EC.ll GO TO 10 


SPLN2 

30 

TFOtJl^ITii) 

- , - 

‘ SP182 

31* 

io continue 


SPLN2 

32 

' RETURN ■ 


' "5PLN2 

33 

END 

. 

SPi.82 

3.4 


117 



SU8«0UTIHt SPJS£T(N»X»r*0*RHS*IFRMS» 


1 


DIHEHSION All) ,Y(1> ,oui 

sw»t 

4 


DAT* JNAX,Rtj,k/2 0, .16666667,1.07179677/ 

SPISfT 

- $ 


0(11=0. 

SH5S? 

3 


If(N.EQ.l) RETURN- 

SPISfT 

t> 


D 1 N 1 = 0 . 

SPSSST 

7 


AN = N 

SPISET 

8 


DO 1 II«2,N 

SPISET 

9 


l=*+2-ll 

SPISET 

10 


>m >=x('i i-x (i-ii 

SPISET 

11 

1 

.v(ii=mii-T( i-ni/xii > 

SPISET 

12 


1MN.EC..2) Gl) TO 5 

SPISET 

13 


DC 2 I 1 = 3 iN 

SPISET 

14 


I=N+3-I I 

SPISET 

15 


D( 1-1 )*2.« <Y( 2 1— Y (1 — Ul/CXIU+X ( 1-1) ) 

SPISET 

16 


YtIS=l. 5*0(1-11 

SPTSET " ’ 

17 

2 

X(Il=0.5*XtI-l)/(x(I»*X( 1-1 J » 

SPISET 

18 


DC 3 J= 1 * JKAX 

SPISET 

19 


DC 3 1 = 3, N 

SPISET 

20 

3 

L ( I-l)=w*fr (I J-X ( I 1 *0 ( I-2H0.5-XU > KOII j )-(R-1.01*D( 1-11 

SPISET 

21 


DC t 1 = 3, N 

SPISET 

22 


X(I)=X(I-11*(<..5/X(I1-1.) 

SPISET 

23 

<t 

Y ( 1 1=0. 33333333* Y( 1 1 * ( X ( 11 *XU-1 1 1 +Y ( 1-1 1 

SPISET 

24 

c. 

SAYt=Y12»-i'fc*> (2)*(2.'*0(11 *dT2> 1 

SPISET 

25 


L»C ‘j l = ^ i N 

SPISET 

26 


C=Rt:*X( I ) 

SPISET 

27 


Sl=Y(i 1-C* 12. *Ct I-il+QU ) 1 

SPISET 

28 


5; = vm,CM0U-lK2.»3(ill 

SPISET 

29 


Cl I~11 = G. SUSAVE+Oil 

SPISET 

30 

O 

SAVt =52 

SPISET 

"31 


0(M = $A*E 

SPISET 

32 


ifUCR-S.M.i) OH T-j 3 

SPISET 

33 


t-N!> = C-. 

SPI5E.T 

34 


IfCN.ta.E) oC TO 6 

SPISET 

3'5 


C£ 7 1 * j f tt 

SPISET 

36 


C=2.*C( l-Il 

SPISET 

37 

7 

Kl'S = Rhi+((C*o< 1—2 1-3. *Y ( I-m/Xl I-11-UD(I)*C-3.*Y(II l/X(I> 1**2 

SPISET 

38 


R*»S=2.*S-2RTUPS/AN1 

TPTSeT 

“T9" 

6 

DO 5 i=2 ,N 

SPISET 

40” 


Y ( 1 )"= Y ("I 1 *X ( I 1 ♦ V U - ll ' 

' SPISET " 

41 

9 

XI I 1 = X( tl+XI 1-11 

SPISET 

42 


RETURN 

SPISET 

43 


tMC 

SPISET 

44 
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SCbnuOliM irtCCTri (N,X, Y,T,NSH0S> 

SMOOTH 

2 






CUBIC weighted by h 

SMOOTH 

-4 


DIMENSION XI 1) t.Y(l> »T( 1) 

SMOOTH 

5 


IPIN.LT.5J RETURN 

SMOOTH 

b 


DO 10 N5 = l i NSMOS 

"swjoth 

7 


t«u=ns 

SMOOTH 

a 


AN =N 

■rnri 

9 


S=(Tl ll*IXtH]-m) 1 /AN) 4*2 

SMOOTH 

10 


DO 4 L*i>N 

SMOOTH" 

jr 


K = KlNOCN-'nHAXO<l,L-2n 

SMOOTH . 

12 


K4=K»4 

~ SMOOTH 

T3 


DO 1 I = l» 20 

SMOOTH 

14 

"1 

~rm~a. ' 

SMOOTH 

13" 


00 3 h=K»K4 

SMOOTH 

16 


W=1 ./TS ♦ fX <i J - XTHI ) 4'$2 > 

T SMOOTH - 

17 


K-l « 0 

SMOOTH 

18 


DO 3 

— SMOOTH 

IT 


14=1—1, 

SMOOTH 

ZO 


-R'Kvr: o 

SMOOTH 

“ 7T 


u-U i J= i * *\ 

SMOOTH 

22 


J4=44j+I4 

SMOOTH 

73 


TC JiJ =1 C J41«R*RR#H 

SMOOTH 

24 

-» 

Rft = R R 4* I N ) 

SMOOTH 

Zb 


T(I*lfc)=T( l*lt)+R4Y(rt)4W 

SMOOTH 

26 

T 

R=ft*XlH> 

SMM ililJi 1 1 WMm 

z / 


CALL ChLi*YtTj4»T(L7)*l*4) 

SMOOTH 

28 


"fi=L-< (L-IT/5)'*i ” - 

5HU0IH 

29 


IFlL.GT.5) Y<L-b) = T<rl«-20) 

SMOOTH 

30 


T l N + 20 1 =0. 

3M00TH 

3T 


R-l.C . 

SMOOTH 

32 


L C 4 a=l»4 

"SMOOTH 

"33" 


T|.l*-0I^T(1*2C)tS«r ( J ♦ 16 ) 

SMOOTH 

34 

'i 

R = tiTTEJ •' 

SMOOTH" ' 

35 


14=6-5 

SMOOTH 

36 


DC 6 L=1,d ' '" 

"“‘SMOOTH' - 

“37 


NL=*-*l-HI-*L-1 » /=> 1^5 

SMOOTH 

38 


J<»=14 + L 


3.9 

5 

Y( J4 >=T(nL*20 1 

SMOOTH 

40 

10 

CONTINUE 

' SMOOTH 

41 


* t TU*v N 

SMOOTH 

42 


ENt’ 

SMOOTH 

43 
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SUBROUTINE ChLS*.Y{i»N,rl,H,NA> 


CHLSKY 

2 



OIHEHSION A(N> ,i) ) 


CHLSKY 

3 

c 


CHULESk-Y CECGPFJS I T ION 13 UStD 13 3Gi.Vt THE MAT R 1 X 

EOUATION A X=d 

chlsky 

A 

c 


tjHERE TnE CCEEflCItMl taTkU, A, U SYMMETRIC. ON 

OUTPUT X IS 

CHLSKY 

5 

c 


STORED IN B 


CHLSKY 

6 



ifltN.fcc.'n goto o 


CHLSKY 

7 



00 2 1*2, N 


CHLSKY 

8 



I 1*1-1 


CHLSKY 

9 



DC 2 J - I , h 


CHLSKY 

10 



DO 2 »_ = i , 1 1 


CHLSKY 

11 


2 

AiI,J > = A( 1, J)-Ml, I»*4(L, J1 /A I L»i-) 


CHLSKY 

12 



00 5 K*l,n 


CHLSKY 

13 



00 J 1 = A»N 


CHLSKY 

16 



11 = 1-1 


CHLSKY 

15 



00 3 l =1,11 


CHLSKY 

16 


3 

B (I »K ) =8'( I,K)-A(L,i)*l(L,K)/AILfLI 


CHLSKY 

17 



CG A 1=2, N 


CHLSKY 

18 



I X = I— 1 


CHLSKY 

19 



DO 6 L=l,li 


CHLSKY 

2 C 



N I = N— 1 1 


CHLSKY 

21 



NL = N+i- l . 


CHLSKY 

22 


S 

"B«NI,K)=B (Nl ,M-tl ’■I ,’L. !='<■'<->■ ) / A ( <m_ , r,L ) 


chlsky 

23 



CC 3 I = l , N 


CHLSKY 

24 


5 

'silt *)=«>( [,Kl/AUni 


CHLSKY 

25 



r t TuRn 


CHLSKY 

26 


6' 

’a<1»1) = iJ/A('a , l’l 


CHLSKY 

27 



00 7 L=i,h 


CHLSKY 

28 


7 

e’lliU='A t i, j. 1 *>( 1 1 


CHLSKY 

29 



kETUhn 


CHLSKY 

30 



tNu 


CHLSKY 

31 
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ORIGINAL' RAGE IS 
OF POOR QUALITY 


T 

C 

F1TIUTA Ik 

T — TEMPORARY ARRAY' FOR SW. INE~ SURFACE FITTING . 

.MM: 


T~ 

C 

*0XU) - iNsesipfeMT VA«c2'Al^ ; |t'’ 



C 

■■ ^~z 

A0HU1 - COMESTH RS DtirfXBfKT- «RS«8I.* 0F *A#0*- 




’fififs wr *s" W&zmmr s & sKiwE^sanras 

DIMENSION TntH,l),AflXIlJ,A6T<lI,AEHTll 

; . smtFX 


c 

" IV-6 

StMPF»> 

y&m 


HPXSH— fi 

MP2*NP1-H 

an 

"tx 

- *-m 


— srjsa^r?! 

HP4-NP3*! 

3«8Ei- 

SU*F1- 

16 


IFtM.EG.Ol Gfl TO 13 
SO 2 I-1,N 

son fi 

susFi - 

"'" 'TLJf 

---*-£ 


— rmw*2T=*5rrn 

TU.NP3I-A8YII1 

SORFl 

SURFl 

20 


ZTtl ,WMr»ASHl 1 1 
00 4 l-l.H 

SURFl 

SURFl 

21 

22 


"ninvo. 

TU,NP1)-1.0 

SURFl 

SURFl 

"73“ 

24 


nNPiyn^T.G 

TtNP2,Il-TTI,HP2I 

SURFl 

SURFl 

25 

2* 


NM1»N-1 


00 6~~ I = 1*"NH1~ 
IP1*I*1 


“ D(T"6~J=IP1,K 

XX-TII ,NP21-T< J.NP2J 

— YY=TU.NP3)-HJ,NP31 

H=XX$XX-*YY*YY 

' T(TTJT^}mX C CT>f) 

6 T<J,IJ = TU,U 

"DO TT=1 >3 

IPN=I+N 

— tnn~j^nr? : 

JPN= J+N 

s rrTPKTjrtrr='c." 

K^HSIN^IUNJitNPS^liJt lilt ,T<1,NP4>) 

YMK^EG.lT GJ 1c 9 

WRITE* IW,22C> 


9 CONJINUfc 

STOK f ' (N+3)-CCtFF ICl'fcNT 
L'C 12 I = 1»NP3 

12 AtM IJ = m ,Ni>4 j“ 

13 If (K1H) le, 18,14 

~14 «IU'U».V0C) - 

*R ITt l !W, 12C1_ Ubttl ! J . 1 = 
IF(N.EC.O) uu TiJ Is 
DC In I=1,N 
16 WSITfcUW.llCJ 
Id C CNT l iHUg 
•It TURN 

110 FCWMAT{iOX,l6,lPjc l9, 

120 EURH4TI62X, 1P3E14.7) 


IN ARRAY A5H 


•JP1,NP3 ) 


liAoHl I ) » Au a t I ) ■> A b Y ( I ) 


7i 


SURFl 
' SOW! — 
SURFl 
SURFl 
SURFl 
— SURFT - 
SURFl 
SURFl 
SURFl 
SDRFI 
SURFl 


200 FORMAT! 1H0,20X,37HHHE«E R( H**2“ IX-XM >>**2» IY-YII ) )**2/ 
_U.H0,10X,bAHUN DIMENSIONLESS COORDINATES - DISTANCE /CHORD 
Z1H0, lJA,lni . IX ,‘mrlll J,10X,4HXII ) , 10X,4HY<I 1* 

3 15X,2hA0, 1 2Xi2HAl«12X,2HA2/) ; 

220 f GRHATI 3X ,26HCiU NUT CONVERGE IN* SURF1//1 
ENO 


2< 

"2t 

3C 

"71 

32 

'3T 

34 

35 

36 

37 

38 


' ' ' SURFl 

"39 

SURFl 

40 

SURFl 

41 

SURFl 

42 

SURFl 

43 

SURFl 

44 

SURFl 

45 

SURFl 

46 

SURFl 

47 

SURFl 

48 

SURFl 

49 

SURFl 

50 

SURFl ' 

51 

SURFl 

52 

SURFl 

53 

SURFl 

54 

SURFl 

55 

SURFl 

56 

SURFl 

57 

SURF-1 

58 

SURFl 

59 

SURFl 

60 

11/ SURFl 

61 

SURFl 

■62 

SURFl 

63 

SURFl 

64 

SURFl 

65 
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noon jr> r> o <*■> o r> 


SUBROUTINE SURF2IZl,Z2,Jl,J2,HXY T VALU»VLUX,VtUY,XI,rl*HI*H,K» SURF 2 .2 

- 'COMPUTE VAtUc Cf SraNr-3uRFAC€' FITTED 0ATA AT A POINT <K,V1 SUKP2 3 

21,22- COORDINATES OF THE POINT WHERE THE FITTED VALUE IS SOUCHT SURF2 A 

VALU - F ITT E0~VaLU1 SOUGHT" “ " " ' ‘ "“SURF? " 5 


VLUX - GRADIENT OF FITTED VALUE IN X 
VLUY GfiADffcTT OF 'fTTTEO VALUE IN Y 

Xlt-YI'HI - ARRAYS FOR KNOWN PROPERTIES IN SPLINE-SURFACE FORM 
' ~N — hUNBEK OF POINTS IN XI, Y1 ARRAYS 
MXY-J X=211J)» Y=22 I J ) WHERE 4«J1,J2 
rtXY» 1 Y*=Z2<4) WHERE 4*41,42 

NXY = .: Y = Z 1 ( 1 ) » X = Z 2 I J ) WHERE J~J1,J2 

DIMlNSIUN Xl ( I),YI I 1) ,HII 1) 

Dint NS I ON Z111J»Z2{1),VALU(1) , VLUX ( 1 ) , VLUY II) 

' nPi=n »i 
NPZ=NPl+l 
NP 3 = NP2 ♦ 1 

IF(NXY.EO.i) X=Zl(li 
IF IflXY.EQ . 2 i Y =1 1 (1 ) 

DC AO J=J1,J2 

iFiFixy rn3vrr~r=Trm ; 

IMnxY.cO.Zl X-Z2I4) 

IFtnxY.NE.O) Gt TO’ ID 
X=Zl<4) 

Y=Z2I j) * ~ ' ‘ ‘ 

10 CONTINUE 

GC Tu TT3 ,"iZ ,TITi K 

11 VLUYI JJ=HI (NP3) 

12 VLOXI JJ=HITNP2} “ ' — " ~ 

13 YALU1 J>=HI 1NP1 )»ril (NP2)*X+rlI ( hP3}* V _ 

IFlN.Evl.01 CO TO AO 

UC 30 1 = 1 , N 

TV^r-xrm 

T Y- Y— Y I (I ) 

H = 'Tx*TX,'TY*TY 

hA-O. ^ 

I F ( h . G f .T) .1 HA'ALtlGOF} " " ~ - - - 

H6=2.*I1» + HAJRHIU 1 

- - c cnrrrTTTTZTTZT J K : 

21 VLUYI 4) =VLUY( J >+rtB*TY 

22 VluxI j) =vluxT0T*H5*Tx ” 7 " ~ 

23 VALJI J1=VAL_U( J1+H1I1 J*HRHA 

30 CONTINUE " ‘ - - - - ' “ 

AO CON TINUE 

END 


SURF2^ 6 

"5URF2 “ * 7 

SURF2 8 

SURFz 9 

SURF2 10 

SURF2 11 

SURF 2 12 

SURF-2 13 

SURF2 14 

SURF2 T5 

SURF2 16 

SURF2 “ 17 

SURF2 18 

‘ SURF2 19 

SURF2 20 

5URF2 "21 

SURF2 22 

SURF2 23 

SURF2 24 

SURFZ 25 

SURF2 26 

SURF 2 ~ * 27 

SURF 2 28 

5URF2 29 

SURF2 30 

SURF2 31 

SURF2 32 

— SURF 2 * 37T 

SURFZ 3A 

5URF2 35 

SURF2 36 

SURF2 37 

SURF2 38 

SURFZ ~T9 

SURF2 40 

SURF2 41 

SURF2 42 

SURF2 43 

SURF2 44 

SURFZ --4T5 

SURF2 46 
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nnn ; n o n . f rvo n'o n] 


ORIGINAL PAGE i& 
OP POOR QUALITY 


FUNCTION C IN (X 1, S J CIN 

SINE AND COSINE INTEGRAL SUBROUTINE ■ CIN' 

__ , CIN 

' " IF CALLED AY THE STATEMENT C=CIN(X,SJ “ CIN 

C AND S ARE TEE INTE GRALS OVER T FROM 1 TO INFINITY OF CIN 

COS 4 XT } /I AND SIN I XT J/T ’ CIN 

. _ CIN 

SOX.O " “ CIN 

X»XI CIN 

IF CX) 1,2'.'/“ CIN 

1 SG--SG CIN 

X«-X ' CIN 

2 X2«X»X CIN 

IF (X-1,0) 3,3,4 ' " * " ' " CIN' 

_ _ CIN 

' "FOR ABSIx) LESS "Than 1 A SfckltS EXPANSION IS USED CIN 

_ CIN 

3 V‘=l U X2 /98. C-cVfa )♦. 0S» X 2* 1 . 0 ) *X2/ 18 . 0-1. 0>*X* I. 5 70 79633 CIN 

U=( (X2/45.0-1. Cl 4X2/24. C+1.0 1 4X2/4 . 0- . 577215665-AL0G t X » CIN 

GO to 5 ' ' ... ■ CIN 

CIN 

FOR AdS(X) GFtATcK Ti!AN 1 APPROXIMATIONS OF HASTINGS ARE USED CIN 

CIN 

4 P=(((X2+'i9. 394U9|9XZ+47. 411538)4X2+8. 4933J6I/IIUX2+21. 3610551 CIN 

1 4X2+?0.37t49fc)4*2+30.038227)4x J CIN 

Q=C UX2+2I.Jd3724]^X2»4V.-7i977S)4x2 + 5Voa9i04)/] (I 1X2+27.1779581 CIN 
I 4X2+U9. 918932)4x2+76. 70787b) *X2) CIN 

CO=COS lx) ' ’ CIN 

SI*SIN lx) CIN 

U=0*CC-P4aI ' CIN 

V=P$CC+G4SI CIN 

"5 s - v * sg ‘ ™ ' ‘ " c I n " 

CIN=U CIN 

RETURN CIN 

ENO - CIN 


2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 
31 

34 

35 

36 
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FUNCTION N5IMER(X.H.L»A»BJ 

— ffTHENSIDF AIH,l»,a<H,i» 


DO 30 I 

U7«r 


1 ?N 


— Z 

00 10 J - 1 « N 
ITT - ^ ahaXI tC,A8S(A<I7JT 
IfiC.EQ.O.O) CO TO 1000 


rr 


HSINER 
“H SI HER 
MSINER 
HSTHER — 
MSItgg 
NSIHER 
NSIHER 


r 

4 

-~r 

6 - 


e 


20 

DO 20 J - 1,N 
Ail.JJ * A 1 1 * J I/C 

NSTHER 

HSINER 

9~ 

10 

r 


fTSTHEK - 

11 

30 

£ < I , J ) « 6 l I , J I/O 

NS I HER 

12 


THN.EO.il SCTTo To5 

NSIHER 

13 


NH » N - 1 

NSIHER 

14 


'DO 200 J - 1,nH 


ii 


C « 0.0 

MSINER 

16 

- 

K * 0 


IT 


DO 40 1 s J,N 

NSIHER 

18 • 


D » ABTT Al I , J 1) 

HSIHEft 

19 


IF (C.GE.D) 00 TO AO 

NSIHER 

20 


_______ 

NSIHER 

21 


C = 0 

NSINER 

22 

“ * A 0‘ 

CONTINUE 

' ' NSIUPR 

23 


IF(K.EQ.0.QR.C.LT.l.E-7> GO TO 1000 

NSINER 

HSINER 

24 


IF(K.FO.J) G"D — TO 73 

25 


00 50 JJ - J,h 

NSINER 

26 


t *= A (J iJ J t 

““NSIHER 

27 


A(J,JJ» = ACR.JJJ 

NSINER 

28 

_ - 50"A(K>JJ> Z 

' HSINER 

' 29 


DC 60 JJ * 1,L 

NSINER 

30 


"t 5 — ant, J3) • 

‘ “ HSINER 

31 


BIJ.JJ) = 81K,JJ) 

NSIHER 

32 

“ STT 

"BTKVJTT™' C * 

~ ~ “ HSINER 


70 

C = A(J,J> 

NSIHER 

34 


JP = J '+ 1 ‘ ‘ 

HSINER 

35 


00 SO JJ = JP,f\ 

HSINER 

36 

80 

A(J,JJ) - AT J, J 1 ) / C" 

NSIHER 

37 

90 

DC IOC J J = 1 ,L 

HSINER 

38 

"Ttrcr 

E>IJ,JJ) - ts ( JV3 JT/C 

H5T HER 

39 


DO 200 I = I,N 

NSIHER 

40 


mr.ctf.J) oD TC" 2 o 0 ■ 

NSIHER 

41 


C = A(I,JJ 

MSINER 

42 


DU 110 JJ = JP ..-1 ' '• 

NSINER 

43 

110 

AII,JJ) * All » j J J - C*AU,JJ1 

NSINER 

44 


DG~l2irJ3 = T,l' ~ 

NSIHER 

45 

12 0 

9(1, Jj) = PIIiJjJ - 

HSINER 

46 

20 C 

rCNtlNUt 

NSINER 

47 

208 

C * 4(N,ri) 

MS I HE k 

48 


if ia:s(c .r-7) o: t„ i^uo 

NSINER 

49 


UP 4 10 J = 1,L 

NSIHER 

50 

■ ■ rnr'TTFnr =~rfr, nvz 

HSI'MFR 

51 


i F IN.iC.U GO TO 2 J-J 

MSINER 

52 


l*C il'G I - i 9 

NSINER 

53 


C =■ 4U.N1 

NSINER 

54 


op .'.*o jj = :,u 

NSINER 

55 

220 

:U iJJ) - _ l i i J j j - ( *tjj) 

H S I M L ^ 

56 

230 

mSTniTK 1 

HS IHE R 

57 


fv t 1 U **r N 

"S InEk 

59 

100C 

nSlrtlK = 2 

NSINER 

59 


K£ TJRN 

NSINER 

60 


£ NO 

HSINER 

61 
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ORIGINAL BAG’iS 1* 
OP POOR QUALITY 


FUNCTJXIM H$IM£C(M,N,L,A,8) 


■iifn*:n€»ni 




COHPLtX A,8,G 
DO 30 1 * 1»N 
C - 0.0 

on 10 j - i , n ' ‘ " ’ 

10 C>AMAX1 1C ,ABS (REAL I At I,J))),ABS ( AIHAG ( A ( l , J ) ) } ) 


IFIC.EQ.O.O) GC TO 10 O 0 
00 20 J_<= 1 ,N 
20 Afl.J) = Alfjl/t 
00 30 J * 1 ,L 

30 BUiJ) = Bir.J>/C ‘ 

IFCM.fcQ.l) GO TO 205 


Nr - n — i 

DO 200 J “_l,hN _ __ _ __ ___ 

K «_0 

D0**0 1 - 37n “ ----- - * - - 

CABS (REAL! A I I,J) ) ) ♦ABS (AIHAG (A(1,J)) J 


K * I 

•”-r»TT ■ -■ 

AO CONTINUE __ __ 

VflK.EQ.O. JR.C.LT . 1 . E -?T 'GO TO' TWO' 
IF(K.EQ.J) GO TO 70 


HSIMEC 


HSIMEC 
HSIMEC 
HSIMEC 
' HSIMEC 
HStHEC 


MSIHEC 

HSIMEC 

MSIHEC' 

HStHEC 

BSIMTC 

MSIHEC 


HSIMEC 

HSTHEC - 

MSIHEC 

HSTHEC 

HSIMEC 

HSIMEC - 

MSIHEC 

H5THFC 

MSIHEC 

HSIMEC " 

MSIHEC 



16 
— 17 
18 
19* 
20 

ZT 

22 

23 

2A 

25 

26 


G =■ A(JtJJ) 

AtJ'iJJ) '= AlK.jJ) ~ - “ 

50 AIK, JO) = G 

"DC 65 “j"3 = ITT - 

G = SIJ.JJ) 

5TT ,T3T" — STkTjTT 

60 b(K,Jj) = G 
70 G = l‘.0/A(J,J) 

OP = J ♦ 1 

OG 80 JO - J'P,K * * 

60 A I J , J J ) = A(J,JJ)$G 

9C Tra“IC'0"Tr^TTC 

100 810,00) = 3lJ,JJ)*C_ 

oo 200 r = ivf “ ' ' ■ ' *" '■ 

IF I I.fcC.J) GO TO 200 

G - A { I , J ) '* ' ' " 

00 110 JO = OP,N 

a'o'Tn , j j i *= Ai 1 , 00 ) - g*ai o,Tn 

DC 120 JJ - -1 ,L 

120 bUi'JJ) ’* 6(1, JO) - G>3'(J,JJ) 

200 CONTINUE 
205 G'V A(N,N) 

IF IA6SIREALIG)) ♦ A3S ( AIMAGI G ) ) .LT . 1.E-7J GO TO 1000 
J - 

210 B IN , J ) = BIN, J)/G 

iFIN.EO.U GO TO 230 
DO 220 I = 1,NM 

OC 220 JJ ^ 1,L " 

2 20 B;i,JJ) = BIT, 00) - A { i , N) *B IN, J J) 


MSIHEC 

HSIMEC 

HSIMEC 

HSTHEC* 

HSIMEC 

HSIMEC 

HSIMEC 

HSIMEC 

H5IHEC 

HSIHEC 

HSIMEC 

HSTHEC* 

MSIHEC 

MSIHEC 

HSIMEC 

HSTHEC 

MSIHEC 

HSIMEC 

MSIHEC 

HSTHEC ' 

HSIMEC . 


HAinti. 

HSIMEC 
MStWC* 
HSIMEC 
' hTTSeT 
HSIMEC 
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